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field.  Individual  microchannels  are  formed  in  a  substrate  and  filled  with  a  1  pm  scale 
doped  porous  polymer  matrix  material  that  provides  the  required  ^-potential  when  wetted 
by  the  corresponding  electrolyte.  The  resulting  micro  electrokinetie  actuator  (MEKA) 
arrays  have  characteristics  that  make  them  potentially  suited  for  practical  sublayer  control 
on  full-scale  aeronautical  and  hydronautical  vehicles.  Several  such  arrays  have  been  fabri¬ 
cated  from  a  basic  three-layer  design.  Essentially  loss-less  frequency  response  has  been 
demonstrated  up  to  10  kHz;  theoretical  bandwidth  is  in  the  MHz  range.  The  final  MEKA- 
5  is  a  full-scale  hydronautical  array  with  25,600  individual  electrokinetie  microactuators 
on  350  pm  spacing  in  a  7  x  7  cm  mylar  tile  with  a  novel  unit-cell  architecture.  The  array 
has  a  MEMS-fabricated  top  layer  with  leadouts  to  unit-cell  processing.  In  addition  to  the 
electrokinetie  microactuators  themselves,  the  MEKA-5  design  is  based  on  an  array  archi¬ 
tecture  that  provides  dramatic  reductions  in  both  the  sensor  and  processing  requirements 
needed  to  achieve  practical  sublayer  control  on  real  vehicles. 


CHAPTER  1 


INTRODUCTION 


1.1  Practical  Applications 

Active  control  for  drag  reduction  on  aeronautical  and  hydronautical  vehicles  is  among  the 
highest-impact  applications  of  microsystems  technology  for  military  and  commercial  pur¬ 
poses,  and  represents  one  of  the  longest-standing  objectives  in  the  field  of  fluid  dynamics. 
A  reduction  in  the  drag  on  an  air  vehicle  of  just  a  few  percent  translates  into  enormous 
system-wide  reductions  in  fuel  weight  and  operations  costs,  and  corresponding  increases 
in  vehicle  range  and  payload  delivery.  Similarly,  flow  control  devices  capable  of  on- 
demand  vortex  generation  over  delta  wings,  strakes,  and  other  control  surfaces  in  fighter 
aircraft  would  permit  dramatic  increases  in  maneuverability,  and  large  reductions  in  radar 
cross-section  through  elimination  of  traditional  control  surfaces.  Such  benefits  of  micro- 
systems-based  flow  control  extend  to  naval  surface  and  undersea  vehicles  as  well,  includ¬ 
ing  ships,  submarines,  and  torpedoes,  as  well  as  to  unmanned  vehicles  used  in  a  variety  of 
applications  such  as  UAV’s,  UCAV’s,  and  MAV’s.  Similar  technologies  could  be  also 
applied  to  supersonic  vehicles  for  aerodynamic  noise  reduction,  and  would  find  use  in  a 
wide  range  of  flow  control  situations  involved  in  propulsion  systems,  spanning  from  con¬ 
trol  of  mixing  in  advanced  gas  turbine  combustors,  to  active  control  of  surge  and  stall  in 
compressors. 

The  prospects  of  actually  achieving  such  revolutionary  advances  in  flow  control  are  far 
more  realistic  than  they  might  at  first  seem.  It  is  precisely  in  the  realm  of  vehicle  control 
that  the  underlying  physics  of  the  problem  at  hand  and  the  inherent  characteristics  of  a 
microsystems-based  approach  are  ideally  matched.  The  thin  boundary  layer  that  exists 
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directly  adjacent  to  the  surface  of  any  practical  vehicle  is  one  of  the  most  powerful  nonlin¬ 
ear  systems  found  in  nature,  capable  of  amplifying  the  effect  of  small  microactuator- 
induced  perturbations  into  very  large  changes  in  the  aerodynamic  forces  and  moments  act¬ 
ing  on  the  vehicle.  Moreover,  the  inherent  problem  of  matching  the  length  and  time  scales 
between  microactuators  and  the  physical  system  being  controlled  makes  the  viscous  sub¬ 
layer  of  a  turbulent  boundary  layer  a  natural  choice  for  microsystems-based  control.  The 
exceedingly  small  length  scales  associated  with  flow  structures  that  are  naturally  present 
in  the  viscous  sublayer  of  turbulent  boundary  layers  are  ideally  matched  to  microscale 
actuators.  However  the  high  performance  required  for  sublayer  control  places  rather  strin¬ 
gent  limits  on  the  types  of  actuators  that  may  be  suitable  for  this  task. 


1.2  Active  Sublayer  Control 

One  approach  for  controlling  the  vehicle  boundary  layer  is  to  exploit  the  structure  and 
dynamics  of  the  streamwise  vortices  that  exist  naturally  in  the  viscous  sublayer.  As  will  be 
seen  in  §2,  these  sublayer  vortices  are  typically  located  at  distances  as  small  as  10  pm 
above  the  vehicle  surface,  are  spaced  about  100  pm  apart,  are  roughly  1mm  long,  and 
advect  past  any  fixed  point  on  the  surface  at  frequencies  up  to  10-100  kHz. 

Early  flow  visualization  studies  performed  in  the  1960s,  and  especially  those  by  Kline  et 
al  (1967),  first  showed  the  presence  of  these  streamwise  streak  structures  in  the  near-wall 
region  inside  the  turbulent  boundary  layer,  and  revealed  their  role  in  the  production  of  tur¬ 
bulence  within  the  boundary  layer.  Those  studies  originally  focused  on  alternating  arrays 
of  low-  and  high-speed  streak-like  structures  at  the  upper  edge  of  the  viscous  sublayer. 
Experiments  by  Corino  &  Brodkey  (1969)  suggested  that  the  presence  of  an  inrush  of  high 
speed  fluid  towards  the  wall,  termed  “sweeping”,  was  involved  in  the  production  of  turbu¬ 
lence.  Subsequent  work  indicated  that  most  of  the  turbulence  production  occurs  within 
the  first  100  viscous  length  scales  above  the  wall  due  to  the  dynamics  of  these  streak-like 
structures  {e.g.,  Kim  et  al  1971).  These  early  discoveries  also  revealed  the  potential  impli¬ 
cations  of  the  streak-like  near-wall  structures,  suggesting  that  control  of  turbulence  pro- 
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duction  in  boundary  layers  could  be  achieved  by  proper  manipulation  of  these  coherent 
structures. 

Subsequent  experiments  (e.g.,  Blackwelder  &  Eckelmann  1979)  found  that  alternating 
pairs  of  counter-rotating  streamwise  vortical  structures  at  the  upper  edge  of  the  viscous 
sublayer  were  the  probable  source  of  the  elongated  low-  and  high-speed  streamwise 
streaks  in  the  near-wall  region.  These  structures  were  found  to  be  responsible  for  low- 
speed  fluid  being  “ejected”  from  the  immediate  vicinity  of  the  wall,  resulting  in  the  “burst¬ 
ing”  process  by  which  relatively  large  amounts  of  such  low-speed  fluid  are  suddenly 
ejected  from  the  near-wall  region.  The  term  “bursting”  was  first  coined  by  Kline  et  al 
(1967)  to  describe  the  overall  process  by  which  the  near-wall  streaks  break  up.  Since  then, 
additional  details  about  the  origin  and  dynamics  involved  in  the  bursting  process  and  its 
role  in  momentum  transport  within  the  boundary  layer  have  been  revealed  by  numerous 
studies.  While  the  key  role  of  the  counter-rotating  streamwise  vortical  structures  in  the 
bursting  process  has  been  firmly  established,  the  precise  dynamics  involved  remain  con¬ 
troversial  even  today,  and  much  research  still  continues  to  consider  the  precise  physics  of 
bursting.  Further  observations  of  the  streamwise  sublayer  vortices  in  the  near-wall  region 
of  turbulent  boundary  layers  and  their  influence  on  turbulent  transport  can  be  found  in 
Cantwell  (1981)  and  Robinson  (1991). 

A  number  of  methods  have  been  examined  to  achieve  skin  friction  reduction  by  either 
directly  or  indirectly  controlling  the  streamwise  sublayer  structures  to  interrupt  or  other¬ 
wise  alter  the  bursting  process  in  wall-bounded  turbulent  flows.  These  have  included  both 
active  as  well  as  passive  means.  Among  the  earliest  passive  methods,  preceeding  the  dis¬ 
covery  of  the  sublayer  vortical  structures  themselves,  was  the  use  of  exceedingly  long- 
chain  polymer  molecules  seeded  in  very  low  concentrations  (on  the  order  of  a  few  parts 
per  million)  in  the  fluid.  The  interaction  of  these  long  molecules  with  the  sublayer  struc¬ 
tures  has  been  shown  reduce  the  skin  friction  as  much  as  80  percent  {e.g.,  Lumley  &  Kubo 
1985;  Hoyt  1990).  This  phenomenon,  first  observed  by  Toms  (1949),  is  used  commer¬ 
cially  in  long  oil  pumping  systems  such  as  the  Trans  Alaska  Pipeline  to  reduce  the  wall 
shear  stress  and  thereby  lower  the  required  pumping  power.  These  polymers,  however. 
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can  only  be  used  in  liquid  flows,  and  their  elfeet  is  largely  degraded  after  the  fluid  has 
been  subjected  to  very  high  shear  rates,  as  oecurs  in  pumps.  More  recent  work  {e.g., 
Chara  et  al  1993;  Miska  et  al  1995;  Warholic  et  al  1999;  Zakin  2001)  has  used  very  low 
concentrations  of  surfaetant  ions  which  self-assemble  into  long-chain  micelle  structures, 
and  in  so  doing  produce  a  similar  effect  as  long-chain  polymer  moleeules.  Unlike  poly¬ 
mers,  however,  after  passing  through  regions  of  high  shear  rates  these  mieelles  rapidly 
reassemble  to  recover  their  original  effeetiveness.  A  related  method  uses  microbubbles, 
introduced  in  the  near-wall  region  of  turbulent  boundary  layers,  to  aehieve  wall  shear 
stress  reduetion.  While  the  preeise  meehanism  by  which  the  drag  reduction  occurs  is  not 
yet  well  understood,  experiments  have  shown  that  wall  shear  stress  reductions  of  80%  or 
more  are  achievable  under  laboratory  eonditions.  As  in  many  drag  reduction  methods, 
however,  a  key  coneem  is  that  the  pumping  power  required  to  produce  sueh  mierobubbles 
on  the  subsurfaee  hull  of  practical  hydronautical  vehicles  can  more  than  offset  any  power 
benefits  from  the  drag  reduction. 

Another  class  of  passive  drag  reduction  methods  has  used  fine-scale  “riblets”  -  eontinu- 
ous  grooves  aligned  in  the  mean  streamwise  direction  -  on  the  surface  to  alter  the  structure 
and  dynamics  of  the  near-wall  region  and  thereby  achieve  wall  shear  stress  reductions. 
Such  riblets  are  typically  located  in  the  viscous  sublayer  of  the  boundary  layer  and  have 
scales  of  the  order  of  the  structures  found  in  that  region.  Walsh  &  Weinstein  (1979)  were 
among  the  first  to  examine  this  teehnique,  and  it  has  been  examined  in  numerous  other 
studies  {e.g,  Walsh  1983;  Johansen  &  Smith  1986;  Lazos  &  Wilkinson  1988;  Choi  1989, 
1993;  Pollard  1997).  Drag  reductions  on  the  order  of  8%  have  been  demonstrated  in  labo¬ 
ratory  experiments  and  in  direct  numerical  simulations,  however  the  precise  meehanism 
by  whieh  the  wall  shear  stress  reduetion  occurs  remains  a  subjeet  of  some  controversy. 

A  number  of  other  passive  control  methods  have  also  been  used  to  demonstrate  manipula¬ 
tion  of  vortieal  struetures  in  the  near-wall  region  of  boundary  layers.  Noteworthy  among 
these,  for  the  present  study,  are  experiments  by  Gad-el-Hak  &  Blackwelder  (1987,  1989), 
which  showed  that  slot  suction  could  prevent  the  bursting  proeess  of  artificially  generated 
structures  in  a  laminar  flow.  In  principle,  this  approaeh  could  be  used  in  turbulent  bound- 
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ary  layers  to  reduce  turbulence  production  by  the  bursting  of  naturally-occurring  stream- 
wise  vortical  structures. 

All  of  these  passive  drag  reduction  methods  attempt  to  reduce  the  skin  friction  in  turbulent 
boundary  layers  by  modifying  the  coherent  structures  present  in  the  viscous  sublayer  with¬ 
out  actively  sensing  the  pattern  of  these  structures  or  actively  reacting  to  changes  in  these 
coherent  structures.  A  number  of  other  studies  have  attempted  to  couple  sensors  and  actu¬ 
ators  to  provide  active  control  of  wall-bounded  turbulent  flows.  Among  the  earliest  of 
these  were  Alshamani  et  al  (1982),  Wilkinson  &  Balasubramanian  (1985),  Nosenchuck  & 
Lynch  (1985),  Nosenchuck  et  al  (1987),  Breuer  et  al  (1989),  and  Nosenchuck  &  Brown 
(1993).  Wilkinson  (1990)  has  reviewed  a  number  of  these  early  attempts  at  active  control 
of  the  near- wall  structures. 

The  extensive  development  of  micro-electro-mechanical-systems  (MEMS)  technologies 
over  the  past  decade  has  opened  a  new  avenue  for  such  active  control  approaches  for  the 
near-wall  region  of  turbulent  boundary  layers  to  achieve  drag  reduction.  MEMS  fabrica¬ 
tion  processes  allow  comparatively  inexpensive  production  of  large,  dense  arrays  of 
microscale  wall  shear  stress  sensors  and  pressure  sensors,  typically  having  length  scales  of 
the  order  of  a  few  hundred  pm.  Such  sensor  arrays  are,  in  principle,  capable  of  detecting 
the  wall  signature  of  the  instantaneous  coherent  structure  pattern  in  the  near-wall  region  of 
turbulent  boundary  layers.  Information  of  this  type  could,  in  principle,  be  used  with  large 
dense  arrays  of  MEMS-fabricated  wall  actuators  and  an  appropriate  control  system  to 
manipulate  the  near-wall  coherent  structures  and  their  dynamics  in  order  to  interfere  with 
the  bursting  process.  Various  attempts  to  develop  sensors,  actuators,  and  control  systems 
to  accomplish  this  goal  are  discussed  by  Bushnell  &  McGinley  (1989),  Eiedler  &  Eemholz 
(1990),  Wilkinson  (1990),  Gad-el-Hak  (1989,  1994,  1996,  2001),  Choi  etal  (1993,  1994), 
Jacobson  &  Reynolds  (1993,  1998),  Moin  &  Bewley  (1994),  McMichael  (1996),  Eumley 
(1996),  Ho  &  Tai  (1996,  1998),  Pollard  (1997),  Eumley  &  Blossey  (1998),  and  Eofdahl  & 
Gad-el-Hak  (1999). 

At  present,  MEMS  sensor  technologies  exist  that  are  capable  of  detecting  these  vortical 


5 


structures  in  the  viscous  sublayer  via  their  wall  shear  stress  signature.  Deteetion  of  high 
speed  sublayer  streaks  in  turbulent  boundary  layer  has  been  demonstrated  in  relatively 
low-speed  wind  tunnel  tests.  A  variety  of  sensors  are  available  today  that  can  accomplish 
this  task,  including  micro  hot-wire  anemometers,  micro  surface  stress  balances,  and  micro 
pressure  sensors;  these  are  discussed  in  detail  in  §6.  Moreover,  the  eontrol  algorithms  and 
assoeiated  processing  electronics  necessary  to  eouple  such  sensors  to  an  appropriate  set  of 
microaetuators  are  also  relatively  well-developed  {e.g.,  Choi  et  al  1993,  1994;  Rebbeck  & 
Choi  2001).  The  required  algorithms  and  processing  capabilities  are  not  very  far  beyond 
what  has  already  been  demonstrated  in  a  number  of  less  demanding  low-speed  wind  tun¬ 
nel  tests.  In  fact,  reduction  of  surface  shear  stress  has  been  successfully  demonstrated 
with  remarkably  simple  control  algorithms  in  several  laboratory-scale  experiments  and 
numerical  simulations  (e.g.,  Breuer  et  al  1989;  Moin  &  Bewley  1994;  Jacobsen  &  Rey¬ 
nolds  1998).  The  prineipal  remaining  teehnical  hurdles  required  to  achieve  microsystem- 
based  aetive  sublayer  control  for  drag  reduction  are  in  the  area  of  microactuators  for 
manipulation  of  the  near-wall  coherent  structures. 


1.3  Previous  Microactuator  Approaches 

Reeent  experimental  and  numerical  studies  specifically  related  to  the  development  of 
mieroaetuators  for  control  of  near-wall  turbulenee  have  been  reviewed  by  Moin  &  Bewley 
(1994),  Gad-el-Hak  (1996,  2001),  MeMichael  (1996),  Ho  &  Tai  (1996,  1998),  Pollard 
(1997),  Lumley  &  Blossey  (1998),  Lofdahl  &  Gad-el-Hak  (1999)  and  Bewley  (2001), 
among  others.  Among  the  earliest  actuation  approaches  for  sublayer  control  was  the  use 
of  eleetromagnetie  fields  to  induce  a  transversally-oscillating  Lorentz  foree  on  the  near¬ 
wall  region  of  turbulent  boundary  layers.  Experiments  by  Nosenchuck  &  Brown  (1993) 
showed  a  reduetion  in  wall  shear  stresses  through  use  of  such  actuation  on  a  fiat  plate  tur¬ 
bulent  boundary  layer.  Similar  results  have  been  obtained  in  experiments  by  Henoeh  & 
Staee  (1995)  and  in  numerical  simulations  by  Crawford  &  Karniadakis  (1997).  Numerical 
simulations  by  Lim  et  al  (1998)  showed  that  sueh  a  transversally-oscillating  applied  mag¬ 
netic  flux  can  inhibit  the  motion  of  the  streamwise  vortiees  and  reduee  their  strength. 
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Simulation  results  by  Berger  et  al  (1999)  showed  similar  wall  shear  stress  reduetions  by 
eleetromagnetie  Lorentz  foree  actuation  in  the  near-wall  region  of  turbulent  channel  flow. 
However,  the  power  required  for  such  actuation  has  been  shown  to  be  at  least  an  order  of 
magnitude  larger  than  the  power  saved  due  to  the  wall  shear  stress  reduction  (Lim  et  al 
1998). 

Carlson  &  Lumley  {\996a,b)  examined  the  effect  of  a  local  wall-normal  deformation 
applied  at  the  wall  -  specifically  a  gaussian  “bump”  -  on  the  near-wall  streamwise  vorti¬ 
ces.  Their  direct  numerical  simulations  showed  that  an  appropriately  moving  bump  could 
produce  a  local  reduction  in  the  wall  shear  stress.  Related  numerical  simulations  of  wall 
deflection,  using  a  combination  of  sensors  and  deformable  wall  actuators  to  manipulate 
the  streamwise  vortices  in  wall-bounded  turbulent  flows,  have  been  performed  by  Mito  & 
Kasagi  (1996)  and  showed  almost  no  drag  reduction.  However,  similar  numerical  experi¬ 
ments  by  Endo  et  al  (1999)  and  Kang  &  Choi  (2000)  exhibited  wall  shear  stress  reductions 
of  10%  and  13-17%,  respectively. 

Ho  &  Tai  (1996,  1998)  and  Tsao  et  al  (1994)  demonstrated  an  electromagnetically  acti¬ 
vated  microfiap  actuator  that  could  be  deflected  to  produce  an  upward  velocity  in  opposi¬ 
tion  to  the  wall-normal  velocity  produced  by  the  streamwise  sublayer  vortices.  This  type 
of  actuator  showed  a  limited  frequency  response  of  about  100  Hz,  and  required  a  relatively 
high  power  consumption  that  would  likely  offset  any  drag  reduction  benefits.  More  recent 
work  by  Sherman  et  al  (1999)  has  attempted  to  solve  some  the  problems  of  this  earlier 
microflap  actuator  by  using  electrostatically  driven,  high  aspect  ratio,  large  deflection,  in¬ 
plane  microactuators.  That  study  has  shown  the  capabilities  of  this  type  of  actuator  to 
induce  a  Stokes’  flow  pattern  in  the  vicinity  of  the  actuator,  but  to  date  this  approach  has 
not  yet  been  applied  to  demonstrate  drag  reduction  in  turbulent  boundary  layers. 

Another  set  of  investigations  has  examined  the  use  of  piezoelectrically-driven  resonant 
cantilever  actuators  to  control  streamwise  vortices  in  turbulent  boundary  layers.  This  type 
of  actuator,  first  studied  by  Jacobson  &  Reynolds  (1993),  produces  a  pair  of  counter-rotat¬ 
ing  streamwise  vortices  moving  away  from  the  wall.  Experiments  by  Jacobson  &  Rey- 
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nolds  (1998)  in  a  laminar  boundary  layer,  using  two  actuators  and  six  sensors,  have  shown 
the  ability  of  this  type  of  actuator  to  displace  synthetically  generated  disturbances. 
Related  to  this,  a  different  type  of  piezoelectric  actuator  has  been  used  by  Rathnasingham 
&  Breuer  (1997)  to  produce  a  control  jet  that  discharges  from  a  slit  into  the  boundary 
layer.  Using  such  an  actuator  with  a  linear  feed-forward  control  algorithm  and  three  sen¬ 
sors  to  detect  synthetically-generated  near-wall  structures  in  a  low  Reynolds  number  tur¬ 
bulent  boundary  layer,  they  achieved  a  maximum  of  31%  reduction  in 

Other  attempts  to  control  streamwise  vortices  include  the  selective  suction  actuators  of 
Myose  &  Blackwelder  (1995).  However,  in  their  experiment  the  suction  was  applied 
between  counter-rotating  pairs  of  Gortler  vortices  in  an  otherwise  steady  laminar  flow. 
Direct  numerical  simulations  of  turbulent  channel  flows  using  a  similar  actuation  method¬ 
ology,  based  on  local  suction  and  blowing  at  the  wall,  was  performed  by  Choi  et  al  (1993, 
1994).  Their  most  successful  result  used  a  control  approach  in  which  suction  or  blowing 
at  the  wall  was  applied  in  such  a  way  as  to  be  exactly  opposite  to  the  wall-normal  compo¬ 
nent  of  the  velocity  at  «  10,  yielding  20-30%  reduction  in  skin  friction.  That  control 
strategy  is  now  commonly  termed  “opposition  control”.  More  recent  numerical  experi¬ 
ments  by  Hammond  et  al  (1998)  showed  that  the  largest  wall  shear  stress  reduction  is 
obtained  by  the  use  of  such  an  opposition-control  strategy  based  on  the  wall-normal  veloc¬ 
ity  aty+  «  15  above  the  wall,  with  the  resulting  average  drag  reduction  being  about  25%. 
They  suggest  that  the  mechanism  for  this  reduction  is  the  creation  of  a  “virtual  wall”  half¬ 
way  between  the  detection  plane  (y^  «  15)  and  the  actuation  plane  (y^  =  0)  that  reduces  the 
sweep  and  ejection  events  thereby  interferes  with  the  transport  of  high  momentum  fluid 
towards  the  wall  and  low  momentum  fluid  away  from  the  wall.  Despite  the  success  of  this 
opposition-control  approach  in  such  numerical  experiments,  it  would  be  difficult  in  prac¬ 
tice  to  measure  the  wall-normal  velocity  at  that  distance  above  the  wall.  Closely  related  to 
these  earlier  studies,  Rebbeck  &  Choi  (2001)  applied  opposition-control  using  wall  jets 
created  with  piston- type  actuators.  Their  experimental  investigation,  which  used  a  sensor 
and  an  actuator  in  a  wind  tunnel  turbulent  boundary  layer,  exhibited  changes  in  the  inten¬ 
sity  of  near-wall  bursting  events  when  the  jet  was  issued  into  the  boundary  layer.  They 
suggest  that  this  mechanism  will  reduce  the  skin-friction  drag  because  the  wall  jet  pre- 


vents  high  momentum  fluid  from  reaehing  the  wall. 


The  studies  noted  above  have  shown  that  active  sublayer  control  of  turbulent  boundary 
layers  is  feasible  under  laboratory  conditions,  and  have  demonstrated  several  types  of 
actuators  suitable  for  active  sublayer  control  under  laboratory  conditions.  However  unlike 
such  laboratory-scale  demonstration  experiments,  the  turbulent  boundary  layers  on  real 
aeronautical  and  hydronautical  vehicles  are  typically  at  much  higher  Reynolds  numbers, 
and  thus  involve  much  smaller  length  scales  and  much  shorter  time  scales  than  are  gener¬ 
ally  achievable  in  laboratory-scale  boundary  layers.  Furthermore,  actuator  arrays  on  real 
vehicles  will  be  subjected  to  a  variety  of  non-ideal  environmental  and  operating  condi¬ 
tions  that  are  not  encountered  under  laboratory  conditions.  Practical  microactuators  for 
use  on  vehicles  will  thus  need  to  be  suflhciently  robust  and  infragile  to  permit  operation 
under  such  conditions.  Finally,  the  power  required  for  large,  dense  arrays  of  such  micro¬ 
actuators  to  achieve  the  desired  level  of  drag  reduction  must  be  sufficiently  small  to  pro¬ 
vide  a  significant  net  reduction  in  the  power  requirements  of  the  vehicle.  All  of  these 
factors  complicate  the  practical  application  of  existing  microactuator  technologies  to  per¬ 
mit  active  sublayer  control  for  drag  reduction  on  full-scale  vehicles  under  realistic  operat¬ 
ing  conditions. 


1.4  Present  Work 

To  date  there  has  been  no  microactuator  technology  demonstrated  that  can  meet  all  of  the 
requirements  noted  above  for  a  practical  microsystem-based  approach  to  active  sublayer 
control  for  drag  reduction  on  aeronautical  or  hydronautical  vehicles  at  full-scale  operating 
conditions.  By  comparison,  the  requisite  wall  shear  stress  and  pressure  sensor  technolo¬ 
gies  are  relatively  mature,  and  the  underlying  control  theory  required  to  link  the  sensors 
and  actuators  is  also  relatively  well  developed.  The  central  technology  required  to  achieve 
practical  vehicle  control  rests  largely  in  the  development  of  robust,  high-performance, 
highly-survivable,  low-power  microactuators  for  integration  in  large,  dense  arrays  that  can 
meet  the  realistic  operating  requirements  noted  earlier. 
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The  present  study  (Diez-Garias  et  al  2000,  Dahm  &  Diez-Garias  2001)  is  the  first  to  exam¬ 
ine  the  eleetrokinetie  principle  as  the  basis  for  a  new  class  of  microactuator  arrays  that  are 
potentially  suitable  for  active  sublayer  control  on  real  aeronautical  and  hydronautical 
vehicles  under  practical  conditions.  The  eleetrokinetie  principle,  discovered  by  Reuss 
(1809)  and  now  used  extensively  in  a  variety  of  applications,  involves  the  interaction 
between  an  electric  field  and  mobile  ions  in  the  “double-layer”  that  forms  naturally  at  the 
phase  interface  between  a  suitable  electrolyte  and  a  suitable  solid  material.  Under  an 
applied  electric  field,  the  mobile  ions  accelerate  as  they  extract  energy  from  the  electric 
field,  and  then  collisionally  transfer  this  kinetic  energy  to  the  molecules  that  make  up  the 
bulk  of  the  electrolyte.  In  this  manner,  an  applied  electric  field  can  induce  bulk  motion  in 
the  electrolyte.  The  resulting  induced  eleetrokinetie  flow  is  typically  too  weak  for  most 
macroscale  applications,  since  the  double-layer  is  extremely  thin.  However,  the  Helm- 
holtz-Smoluchowski  scaling  laws  that  govern  such  eleetrokinetie  flow  indicate  significant 
performance  advantages  in  microscale  applications.  As  a  result,  this  eleetrokinetie  flow 
potentially  provides  a  highly  efficient  means  for  moving  comparatively  small  amounts  of 
fluid  on  comparatively  fast  time  scales  in  comparatively  small  geometries.  In  this  study 
these  scaling  principles  are  used  together  with  microfabrication  technologies  and  funda¬ 
mental  considerations  of  the  sublayer  structure  and  dynamics  to  develop,  fabricate,  and 
demonstrate  microactuators  based  on  such  eleetrokinetie  flow  that  can  meet  many  of  the 
requirements  for  active  sublayer  control  in  turbulent  boundary  layers  on  real  vehicles. 

The  eleetrokinetie  microactuator  arrays  that  are  the  subject  of  this  study  seek  to  interrupt 
the  bursting  process  associated  with  naturally-occurring  streamwise  sublayer  vortices  in 
the  turbulent  boundary  layer  under  conditions  that  are  typical  of  real  vehicles.  The  sub¬ 
layer  bursting  process  sets  the  local  rate  of  momentum  transport  from  the  vehicle  to  the 
fluid,  and  thus  determines  the  drag  that  acts  on  the  vehicle.  As  noted  earlier,  the  bursting 
process  can  be  interrupted  by  manipulating  the  streamwise  vortices  in  any  of  a  number  of 
ways.  In  the  present  microactuator  arrays,  each  individual  actuator  serves  as  a  local  vol¬ 
ume  source,  applied  over  a  brief  time  interval,  at  the  wall.  The  actuator  thus  impulsively 
displaces  a  fixed  volume  of  fluid  -  either  positive  or  negative,  corresponding  to  local  suc¬ 
tion  or  blowing  -  between  the  wall  and  the  viscous  sublayer  to  induce  lateral  displace- 
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ments  in  the  streamwise  sublayer  vortiees.  The  volume  displacement  produced  by  each 
actuator  occurs  via  electrokinetic  pumping  in  response  to  an  impulsively-applied  electric 
field.  Dense  arrays  consisting  of  large  numbers  of  such  closely-spaced  electrokinetic 
microactuators  have  been  designed  and  fabricated  in  this  study  at  scales  suitable  for  active 
sublayer  control  at  real-vehicle  conditions. 

Furthermore,  a  complete  system  architecture  for  such  microactuator  arrays  has  also  devel¬ 
oped  as  part  of  this  study.  This  novel  system  architecture  is  based  on  groups  of  actuators 
within  a  local  “unit  cell”  operating  together  in  a  “response  vector”  that  depends  on  the  sig¬ 
nal  from  a  group  of  collocated  wall  shear  stress  sensors.  These  local  wall  shear  stress  sen¬ 
sors  determine  the  “state  vector”,  characterizing  the  wall  shear  stress  signature  produced 
by  the  local  streamwise  sublayer  vortex  structures  in  the  unit  cell.  This  unit-cell  architec¬ 
ture  further  permits  a  greatly  simplified  control  logic  approach,  also  developed  as  part  of 
this  study,  that  is  based  on  a  simple  look-up  table  resident  in  firmware  within  each  unit 
cell.  This  lookup  table  provides  the  proper  unit-cell  response  vector  for  the  actuators  in 
response  to  certain  state  vectors  from  the  unit-cell  sensors  that  correspond  to  pairs  of 
streamwise  sublayer  vortices  in  configurations  that  correspond  to  incipient  bursting. 

This  dissertation  is  organized  as  follows.  Chapter  2  first  reviews  the  essential  aspects  of 
the  structure  of  turbulent  boundary  layers,  and  then  from  these  determines  the  required 
spacing  between  individual  electrokinetic  microactuators,  the  frequency  response  required 
of  such  microactuators,  and  the  equivalent  steady  flow  rate  required  of  such  actuators  for 
active  sublayer  control  on  a  wide  variety  of  full-scale  aeronautical  and  hydronautical  vehi¬ 
cle  types.  Chapter  3  then  presents  the  system  architecture  developed  in  this  study  based 
on  the  unit-cell  concept  for  large,  dense  arrays  of  such  electrokinetic  microactuators  with 
collocated  wall  shear  stress  sensors  and  control  logic.  Chapter  4  describes  essential  back¬ 
ground  material  related  to  various  electrokinetic  phenomena,  and  then  uses  this  to  obtain 
results  for  the  theoretical  steady  and  unsteady  performance  of  electrokinetic  microactua¬ 
tors.  Chapter  5  presents  results  from  microfabrication  and  testing  of  several  generations 
of  such  micro-electro-kinetic-actuator  (MEKA)  arrays,  ranging  from  simple  individual 
capillary  actuators  to  full-scale  hydronautical  actuator  arrays  with  25,600  individual  elec- 
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trokinetic  microactuators  on  a  7  x  7  cm  conformally  applicable  tile  (MEKA-5).  This 
ehapter  also  presents  results  from  wind  tunnel  experiments  demonstrating  lateral  displaee- 
ment  of  synthetically-generated  streamwise  vortiees  via  volume  displaeements  indueed  at 
a  wall.  Chapter  6  presents  design  results  for  wall  shear  stress  sensors  suitable  for  integra¬ 
tion  in  sueh  eleetrokinetie  mieroaetuator  arrays,  together  with  a  mierofabrieation  proeess 
sequenee  for  produeing  large,  dense  arrays  of  sueh  eolloeated  sensors  in  an  eleetrokinetie 
mieroaetuator  array.  Chapter  7  summarizes  major  findings  from  this  study  and  presents 
eonelusions  as  to  the  applieability  of  the  eleetrokinetie  prineiple  as  the  basis  for  a  practieal 
elass  of  mieroaetuators  suitable  for  drag  reduetion  via  aetive  sublayer  eontrol  on  full-seale 
aeronautical  and  hydronautical  vehicles.  An  Appendix  provides  eertain  results  relevant  to 
the  performanee  of  slot-type  eleetrokinetie  mieroaetuator  arrays. 


12 


CHAPTER  2 


MICRO  ACTUATOR  PERFORMANCE  REQUIREMENTS 


As  noted  briefly  in  §1,  the  bursting  proeess  assoeiated  with  streamwise  vortiees  at  the 
outer  edge  of  the  viseous  sublayer  sets  the  rate  of  momentum  transport  from  the  wall  to 
the  fluid,  and  hence  the  drag  that  acts  on  the  vehicle.  Drag  reduction  can  be  accomplished 
by  acting  on  the  streamwise  vortices  to  interrupt  or  otherwise  interfere  with  this  sublayer 
bursting  process.  For  the  present  microactuator  arrays,  the  individual  actuators  serve  as 
point  volume  sources  that  impulsively  displace  a  fixed  volume  of  fluid  in  the  near-wall 
region  of  the  flow.  Each  actuator  thus  produces  locally  positive  or  negative  volume  dis¬ 
placements  over  a  brief  time  interval  to  displace  the  streamwise  sublayer  vortices  along 
the  spanwise  direction.  In  order  to  determine  the  performance  requirements  of  individual 
microactuators  needed  to  achieve  such  lateral  displacement  of  the  sublayer  vortices,  this 
section  first  reviews  essential  elements  of  the  structure  and  dynamics  of  the  near-wall 
region  of  turbulent  boundary  layers.  It  then  uses  this  information  to  derive  the  actuator 
spacing,  frequency,  and  flow  rate  requirements  necessary  to  achieve  adequate  displace¬ 
ment  of  individual  sublayer  vortical  structures. 


2.1  Structure  of  Turbulent  Boundary  Layers 

Turbulent  boundary  layers  obey  a  universal  structure  and  scaling  under  so-called  “equilib¬ 
rium”  conditions,  which  typically  apply  over  the  vast  majority  of  a  vehicle  surface.  This 
makes  it  possible  to  approach  the  problem  of  vehicle  control  with  a  general  set  of  microac¬ 
tuator  arrays,  rather  than  having  to  resort  to  a  separate  approach  for  each  vehicle  type.  It  is 
necessary  to  be  aware  of  at  least  the  basics  of  this  general  structure  and  its  scaling  in  order 
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to  understand  the  approaeh  being  taken  here,  and  to  understand  the  performanee  require¬ 
ments  that  they  dictate  for  microactuator  arrays.  This  section  therefore  reviews  certain 
key  aspects  of  the  structure  and  dynamics  of  equilibrium  turbulent  boundary  layers. 

The  central  issue  concerns  the  form  of  the  mean  velocity  profile  u{y)  with  increasing  dis¬ 
tance  y  from  the  vehicle  surface  (referred  to  generically  as  the  “wall”,  located  at  y  =  0). 
There  are  two  distinct  mechanisms  that  contribute  to  momentum  transport  from  the  wall  to 
the  fluid  and  then  within  the  fluid  above  the  wall.  The  first  is  classical  molecular  diffu¬ 
sion,  with  the  momentum  diffusivity  being  the  “kinematic  viscosity”  v,  and  the  second  is 
turbulent  transport  by  velocity  fluctuations  in  the  flow  itself.  Directly  at  the  wall,  the  fluid 
cannot  move  and  thus  the  turbulent  transport  there  must  vanish.  Consequently,  for  a  small 
distance  above  the  wall  the  turbulent  transport  must  therefore  remain  essentially  negligi¬ 
ble.  For  this  reason,  there  is  a  layer  of  fluid  adjacent  to  the  wall  within  which  molecular 
transport  far  exceeds  turbulent  transport.  With  increasing  distance  from  the  wall,  the 
momentum  flux  density  due  to  molecular  transport  must  remain  roughly  constant,  while 
that  due  to  turbulent  transport  increases,  and  thus  turbulent  transport  eventually  becomes 
the  dominant  mechanism.  The  boundary  layer  therefore  naturally  has  a  two-layer  struc¬ 
ture  -  an  “inner  layer”  directly  adjacent  to  the  wall  within  which  molecular  transport  dom¬ 
inates,  and  an  “outer  layer”  within  which  turbulent  transport  dominates.  The  crossover 
between  these  two  transport  mechanisms  occurs  in  the  outer  part  of  the  inner  layer  and  the 
inner  part  of  the  outer  layer  -  called  the  “log  layer”  for  reasons  that  will  become  evident 
subsequently. 

Within  the  inner  layer,  owing  to  the  dominance  of  diffusion  the  momentum  flux  density 
must  be  roughly  constant  at  a  value  called  the  “wall  shear  stress”  literally  the  value  of 
the  shear  stress  at  the  wall.  All  processes  within  this  layer  are  assumed  to  be  properly 
scaled  with  this  wall  shear  stress,  with  the  viscosity  v,  with  the  fluid  density  p,  and  with 
the  distance  y  above  the  wall.  It  is  convenient  to  equivalently  express  such  near-wall  scal¬ 
ings  in  terms  of  the  characteristic  velocity  defined  from  these  quantities,  namely  u^.  = 
often  called  the  “friction  velocity”,  and  the  characteristic  length  l^  =  v  /u^  termed 
the  “wall  unit”  or  the  “viscous  length  scale”.  The  appropriately  scaled  velocity  profile  is 
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then  u^iy^),  where  =  ulu^  and  =ylU-  The  velocity  profile  throughout  the  inner  layer, 
when  scaled  in  this  manner,  is  universal  for  all  equilibrium  turbulent  boundary  layers. 

Within  the  outer  layer,  however,  the  velocity  profile  is  not  entirely  universal,  and  varies 
somewhat  depending  on  the  local  boundary  layer  thickness  5  and  the  local  streamwise 
pressure  gradient  dp/dx.  The  latter  is  determined  via  the  Bernoulli  equation  by  the  local 
Lagrangian  acceleration  of  the  fluid  outside  the  boundary  layer  as  it  negotiates  around  the 
vehicle  shape.  Thus  irrespective  of  the  particular  vehicle,  the  outer  layer  structure  (under 
equilibrium  conditions)  can  be  cataloged  in  terms  of  an  appropriately  scaled  pressure  gra¬ 
dient  parameter  TL  Various  ways  of  expressing  this  scaling  have  been  proposed,  with  the 
most  successful  being  the  “law  of  the  wake”  (Coles  1956).  In  the  overlap  region,  namely 
the  inner  part  of  the  outer  layer  and  the  outer  part  of  the  inner  layer,  the  velocity  profile 
must  become  independent  of  the  pressure  gradient,  and  this  turns  out  to  demand  that  the 
velocity  profile  in  this  overlap  region  must  have  a  logarithmic  form  -  hence  the  term  “log 
layer”.  The  remainder  of  the  outer  layer  is  referred  to  as  the  “defect  layer”,  because  its 
shape  is  described  in  terms  of  the  departure  (or  defect)  it  demonstrates  in  relation  to  the 
shape  that  would  result  if  the  logarithmic  scaling  continued  there. 

Excluding  the  logarithmic  layer,  the  inner  layer  extends  to  about  «  30,  with  the  region 
from  0  <  <  10  termed  the  “viscous  sublayer”,  and  the  remainder  called  the  “buffer 

layer”.  Thus,  the  viscous  sublayer  is  about  10  wall  units  thick,  and  the  inner  layer  up  to 
the  log  layer  is  about  30  wall  units  thick.  The  physical  thicknesses  that  these  correspond 
to  depend  on  the  value  of  ,  and  hence  on  the  wall  shear  stress  (and  thus  on  the  boundary 
layer  thickness  5  and  pressure  gradient  parameter  Tt),  the  fluid  density,  and  the  viscosity. 
These  vary  depending  on  the  location  on  the  vehicle,  its  size,  its  speed,  and  the  properties 
of  the  fluid  through  which  it  moves. 

Due  to  the  dominance  of  diffusive  transport  very  close  to  the  wall  the  vorticity  there  must 
be  uniform,  and  hence  the  velocity  must  increase  linearly  with  distance  from  the  wall. 
Within  the  viscous  sublayer  {i.e.,  below «  10),  the  form  of  M^(y^)  must  thus  simply  be 
=  y^.  This  linear  velocity  profile  remains  steady  until  the  resulting  uniform  vorticity  layer 
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adjacent  to  the  wall  reaehes  a  thiekness  beyond  whieh  the  no-slip  requirement  at  the  wall 
is  no  longer  able  to  maintain  stability  of  the  layer.  This  instability  is  governed  by  the  Rey¬ 
nolds  number  RCy  =  uyN.  Reeognizing  that  in  wall  variables  RCy  =  then  gives  RCy  = 
y^  within  the  sublayer,  and  since  sueh  shear-driven  instabilities  typically  begin  at  Re  « 

'y 

0(10  ),  this  would  indicate  that  this  uniform  layer  of  vorticity  at  the  wall  will  remain  sta¬ 
ble  up  to  about  «  0(10).  This  is  in  good  agreement  with  the  lower  edge  of  the  buffer 
layer. 

This  two-layer  structure  of  turbulent  boundary  layers  is  illustrated  in  Fig.  2.1.  In  particu¬ 
lar,  Fig.  2.1a  shows  typical  boundary  layer  velocity  profiles  in  outer-layer  variables  (the 
local  boundary  layer  thiekness  5  and  the  loeal  veloeity  U  at  the  “edge”  of  the  boundary 
layer)  for  pressure  gradient  parameters  Tl  ranging  from  mildly  adverse  to  strongly  favor¬ 
able.  Figure  2.\b  shows  the  same  three  velocity  profiles  in  inner-layer  variables  W"  andy^, 
where  the  linear  and  logarithmie  scalings  noted  above  ean  be  readily  identified  in  the 
semi-logarithmie  axes. 

Above  the  upper  edge  of  the  viscous  sublayer  aty^  «  10,  naturally-oceurring  infinitesimal 
disturbanees  eause  the  otherwise  uniform  vorticity  layer  to  develop  sinusoidal  instabilities 
along  the  spanwise  (cross-stream)  direction,  denoted  z.  Amplifieation  and  nonlinear 
growth  of  this  spanwise  instability  leads  to  the  formation  of  elongated  concentrations  of 
vorticity  oriented  along  the  streamwise  direetion,  denoted  x.  As  shown  in  Fig.  2.2,  the 
resulting  streamwise  vortieal  struetures  at  the  outer  edge  of  the  viscous  sublayer  have  a 
spaeing,  in  the  sealed  spanwise  coordinate  z^  =  z/^ ,  of  about  z^  «  100.  The  length  of  these 
“sublayer  vortices”  is  typically  about  «  1000,  but  varies  from  about  400  -  1500.  In 
other  words,  the  sublayer  vortiees  are  about  1000  wall  units  (4)  long,  and  are  spaeed  about 
100  wall  units  apart.  Sinee  the  vortiees  are  loeated  roughly  10  wall  units  above  the  wall, 
they  will  adveet  at  a  speed  «  10,  and  thus  move  over  any  fixed  point  on  the  wall  with  a 

'y 

frequency  of  about  /+  «  10 

Figure  2.2  is  an  illustration  of  the  pattern  of  these  streamwise  sublayer  vortieal  struetures 
in  the  near-wall  region  of  turbulent  boundary  layers.  In  partieular.  Fig.  2.2a  shows  a  typi- 
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cal  example  of  the  instantaneous  strueture  of  these  vortiees,  sueh  as  has  been  doeumented 
in  numerous  direct  numerical  simulations  of  turbulent  boundary  layers.  Note  that  the 
struetures  are  not  spaeed  at  preeisely  equal  intervals,  nor  do  they  have  precisely  equal 
lengths,  however  it  is  eommon  to  idealize  this  pattern  as  shown  in  Fig.  22b  with  the  char- 
aeteristie  dimensions  as  noted  above. 


2.2  Bursting  of  Sublayer  Vortices 

It  is  these  streamwise  vortices  at  the  outer  edge  of  the  viseous  sublayer  that  are  the  key  to 
momentum  transport  from  the  inner  layer  to  the  outer  layer  and  thus  from  the  wall  to  the 
inner  layer,  and  henee  determine  the  drag  that  aets  on  the  vehiele.  The  detailed  dynamies 
of  the  near-wall  struetures  are  eomplex  and  remain  the  subjeet  of  eonsiderable  eontro- 
versy,  however,  the  essential  mechanism  can  be  understood  within  the  eontext  of  a  simple 
inviscid  deseription.  Under  the  influenee  of  their  mutual  Biot-Savart  indueed  motion  as 
well  as  that  of  their  images  on  the  opposite  side  of  the  wall,  these  streamwise  vortiees 
must  naturally  group  into  pairs  having  alternating  sense  of  eireulation.  Eaeh  streamwise 
vortex  pair  remains  nominally  steady  until  the  eombined  indueed  motion  of  the  two  elon¬ 
gated  vortiees  abruptly  lifts  them  away  from  the  near-wall  region  and  into  the  outer  layer, 
in  a  proeess  termed  “bursting”.  In  doing  so,  low  momentum  fluid  in  the  wall-fixed  eoordi- 
nate  frame  (high  momentum  fluid  in  an  absolute  frame)  is  transported  from  the  inner  layer 
to  the  outer  layer,  at  a  mueh  greater  rate  than  would  be  the  ease  by  diffusion  alone.  It  is 
due  to  this  bursting  meehanism  and  the  aeeompanying  greatly  inereased  rate  of  momen¬ 
tum  transport  from  the  wall  to  the  fluid  that  turbulent  boundary  layers  produee  far  greater 
drag  on  a  vehiele  than  do  laminar  boundary  layers.  The  bursting  proeess  is  reeognized  as 
the  main  source  of  turbulenee  production  in  turbulent  boundary  layers  (Kim  et  al  1971; 
Blaekwelder  &  Eekelman  1979;  Cantwell  1981;  Robinson  1991). 

While  it  is  aeeepted  that  the  streamwise  vortiees  are  eritieal  in  this  proeess,  the  evolution¬ 
ary  dynamies  behind  the  bursting  proeess  are  not  yet  elear.  Earlier  reviews  of  the  bursting 
proeess  (Robinson  1991  and  referenees  therein)  explained  that  low-speed  streaks  of  fluid. 
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generated  by  streamwise  vortiees,  were  ejeeted  outward  from  the  wall  after  going  through 
some  transitional  oscillations,  generating  most  of  the  turbulence  production  in  the  bound¬ 
ary  layer.  Further  studies  of  this  process  during  the  last  decade  aided  but  the  use  of  DNS 
simulations  has  modified  some  of  those  earlier  concepts.  DNS  simulations  (Hamilton  et 
al  1995)  have  proposed  a  mechanism  for  the  regeneration  of  near-wall  turbulence  struc¬ 
ture.  This  includes  the  formation  of  streaks  by  streamwise  vortices,  and  the  breakdown  of 
streaks  and  regeneration  of  the  streamwise  vortices.  The  formation  of  high  and  low  speed- 
streaks  is  generated  by  the  “sweeping”  and  “ejection”,  respectively,  of  near-wall  fluid  by 
the  coherent  structures  towards  the  wall  or  away  from  the  wall.  This  process,  also 
described  by  Jeong  et  al  (1997)  as  cross-stream  advection  generated  by  the  vortical  struc¬ 
tures,  is  a  major  source  of  negative  and  positive  Reynolds  stress  in  the  boundary  layer. 
Schoppa  &  Hussain  (2000a,  h)  have  indicated  that  the  wallward  motion  enhances  the  drag 
and  heat  transfer  by  steepening  the  wall  gradient  of  streamwise  velocity.  This  has  been 
confirmed  by  Orlandi  &  Jimenez  (1994).  Further  simulations  by  Jeong  et  al  (1997)  have 
shown  that  the  coherent  structures  are  elongated  quasi-streamwise  vortices  inclined  nine 
degrees  in  the  vertical  plane  and  tilted  in  the  horizontal  plane.  Furthermore,  they  indicate 
that  vortices  of  alternating  sign  overlap  in  the  streamwise  direction.  Schoppa  &  Hussain 
(2000a, h)  have  successfully  used  this  later  model  to  show  that  combined  transient  growth 
and  instability  of  lifted,  irrotational  low-speed  streaks  generates  new  streamwise  vortices. 
Thus  completing  the  regeneration  cycle  of  the  turbulent  fluctuations  in  the  boundary  layer. 
These  recent  simulations  have  shown  that  “bursting”  process  does  not  require  a  “violent” 
break  up  of  the  streaks.  Instead,  this  is  rather  a  dynamic  process  by  which  growing  distur¬ 
bances  of  ejected  low-speed  streaks  generate  new  streamwise  vortices,  and  the  mechanism 
by  which  these  new  streamwise  vortices  extract  energy  from  the  mean  flow  to  create  alter¬ 
nating  low-  and  high-speed  streaks  {e.g,  Jimenez  &  Pinelli  1999;  Gottero  &  Onorato  2000; 
Schoppa  &  Hussain  2000a, h). 


2.3  Sublayer  Control  Principles 

Interrupting  this  delicate,  nonlinear,  Biot-Savart- induced  lifting  mechanism  between  adja- 
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cent  streamwise  sublayer  vortices  before  bursting  ean  occur  has  long  been  an  objective  in 
fluid  meehanies  and  aerodynamies.  As  noted  in  §  1 ,  it  has  been  demonstrated  both  in  labo¬ 
ratory  experiments  and  in  numerical  simulations  that  this  can  be  aecomplished  by  acting 
on  the  streamwise  vortices  in  any  of  a  number  of  ways  so  as  to  offset  the  Biot-Savart 
induced  velocities  that  lead  to  the  sudden  bursting.  For  the  present  microactuator  arrays, 
the  individual  actuators  serve  as  point  volume  sourees  that  displace  (either  positively  or 
negatively)  a  fixed  volume  of  fluid  between  the  wall  and  the  inner  layer.  This  is  shown 
sohematieally  in  Fig.  2.3.  The  motion  that  this  volume  source  field  induces  serves  to  dis¬ 
place  the  streamwise  vortiees  along  the  spanwise  direction,  thereby  loeally  increasing  or 
deereasing  the  separation  between  the  two  elongated  vortices  in  a  streamwise  vortex  pair. 
This  counteracts  the  Biot-Savart  indueed  motion  of  the  vortices  on  each  other,  thus  inter¬ 
rupting  the  bursting  process. 

As  for  a  closed-loop  feedbaek  control  approach  that  drives  an  array  of  sueh  individually 
addressable  “point”  microaetuators,  there  are  two  schools  of  thought  regarding  the  best 
way  to  do  this.  The  first,  whieh  may  be  termed  the  “deterministic”  approach,  is  to  use  sen¬ 
sors  to  identify  the  wall  shear  stress  or  wall  pressure  signature  of  a  pair  of  streamwise  vor¬ 
tiees  that  are  about  to  burst,  and  then  use  colloeated  actuators  to  manipulate  the  vortiees  in 
such  a  way  as  to  delay  or  defeat  the  bursting  proeess.  An  example  of  this  is  the  “opposi¬ 
tion  control”  approach  (Choi  et  al  1993,  1994;  Hammond  et  al  1998;  Rebbeck  &  Choi 
2001)  discussed  briefly  in  §1.  The  seeond,  which  might  be  termed  a  “stoehastie” 
approach,  ignores  the  detailed  dynamics  of  the  bursting  proeess  itself  and  instead  uses  an 
artifieial  intelligence  procedure.  In  that  case,  a  neural  network  between  a  group  of  sensors 
and  actuators  is  trained  over  time  in  such  a  way  as  to  minimize  the  wall  shear  stress.  Both 
approaches  have  met  with  considerable  and  roughly  eomparable  suceess  in  laboratory 
experiments  and  in  numerical  simulations,  demonstrating  reductions  in  the  range  of  20- 
30%  in  the  wall  shear  stress. 

The  seenario  outlined  above  involves  closed-loop  control  of  the  equilibrium  structure  and 
dynamics  of  the  boundary  layer.  An  alternative  uses  open-loop  control  of  the  nonequilib¬ 
rium  dynamics  in  the  near-wall  region  of  the  boundary  layer.  This  is  motivated  by  experi- 
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mental  and  numerical  observations  showing  that,  when  a  boundary  layer  is  subjected  to  a 
sudden  transverse  pressure  gradient  or  crossflow,  the  flow  experiences  temporary  reduc¬ 
tions  in  its  turbulent  transport.  This  suggests  that  it  may  be  possible  to  sustain  this  effect 
by  maintaining  the  boundary  layer  in  a  continual  nonequilibrium  state.  One  way  this 
might  be  done  is  through  spanwise  oscillations  of  the  wall  relative  to  the  streamwise  vorti¬ 
ces  {e.g.,  Jung  et  al  1992,  Akhavan  et  al  1993,  Laadhari  et  al  1994).  Experiments  and 
numerical  simulations  have  shown  that  this  produces  significant  reductions  in  turbulent 
transport  and  thereby  reductions  in  the  wall  shear  stress  of  up  to  40%.  The  effectiveness 
depends  on  the  oscillation  frequency,  with  optimal  drag  reduction  again  occurring  at « 
10  .  Moreover  such  an  approach  does  not  require  any  collocated  array  of  sensors  to 
determine  the  instantaneous  pattern  of  sublayer  vortices,  or  control  logic  to  discern  from 
this  the  appropriate  set  of  actuator  signals.  The  only  requirement  is  for  an  occasional  sen¬ 
sor  to  measure  the  average  wall  shear  stress,  which  is  needed  to  properly  scale  the  oscilla¬ 
tion  frequency  to  achieve  optimal  drag  reductions.  This  “oscillating  wall”  approach  thus 
potentially  provides  a  means  for  open-loop  control  of  the  sublayer  vortices,  and  thus  obvi¬ 
ates  the  difficulties  associated  with  collocating  large  arrays  of  sensors,  actuators,  and  local 
processing.  In  practice,  it  is  the  fluid  above  the  wall  that  would  be  oscillated  by  a  set  of 
microactuators  designed  to  generate  volumetric  displacements  above  the  wall  and  thereby 
induce  oscillatory  spanwise  motions  of  the  fluid  in  the  inner  layer.  However,  while  the 
inherent  simplicity  of  open-loop  control  based  on  this  oscillating  wall  approach  is  appeal¬ 
ing,  Appendix  A  shows  that  it  cannot  be  practically  implemented  with  the  electrokinetic 
microactuators  that  are  the  subject  of  this  study.  This  study  therefore  focuses  on  closed- 
loop  control  using  electrokinetic  microactuator  arrays  with  a  collocated  array  of  sensors 
and  appropriate  control  processors. 


2.4  Actuator  Spacing,  Frequency,  and  Flow  Rate  Requirements 

Irrespective  of  the  control  approach  used  and  the  microactuator  type  it  implies,  it  is  neces¬ 
sary  to  have  an  array  of  actuators  that  allow  the  streamwise  sublayer  vortices  to  be  acted 
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upon  appropriately.  Key  performanee  requirements  involve  the  microaetuator  spaeing, 
frequency,  and  flow  rate  needed  to  achieve  adequate  displacement  of  individual  sublayer 
vortical  structures.  For  closed-loop  control  approach  based  on  individual  point  microactu¬ 
ators,  from  the  forgoing  summary  of  the  sublayer  vortex  structure  and  dynamics  it  is 
apparent  that  such  actuators  must  be  separated  by  typically  100  wall  units,  and  displace  a 
volume  of  fluid  with  an  equivalent  hemispherical  radius  of  the  order  of  10  wall  units  4 , 
with  a  step  response  that  corresponds  to  a  scaled  frequency  of/  =  10  . 


With  this  information,  it  is  possible  to  determine  the  size  and  performance  requirements  to 
which  these  conditions  correspond  for  various  types  of  vehicles.  These  performance 
requirements  will  depend  on  the  fluid  type  (principally  air  or  water,  which  differ  in  their 
density  and  viscosity),  the  vehicle  speed,  and  the  vehicle  length.  More  precisely,  the  uni¬ 
versal  structure  and  scaling  of  equilibrium  turbulent  boundary  layers  shows  that  the  local 
sublayer  properties,  and  hence  the  local  microaetuator  performance  requirements,  depend 
on  the  local  boundary  layer  thickness  5,  the  local  fluid  speed  U  at  the  “edge”  of  the  bound¬ 
ary  layer,  the  local  pressure  gradient  dp/dx,  and  the  fluid  properties  p  and  p.  Of  these,  the 
direct  effect  of  the  pressure  gradient  is  relatively  weak;  it  is  the  indirect  effect  of  the  pres¬ 
sure  gradient  on  the  boundary  layer  thickness  5  that  dominates.  This  permits  an  analysis 
for  each  vehicle  type  without  requiring  a  detailed  boundary  layer  calculation  for  a  specific 
vehicle  geometry.  In  effect,  the  performance  requirements  can  be  obtained  at  any  stream- 
wise  position  on  the  vehicle  for  a  variety  of  pressure  gradient  parameters  Tl  ranging  from 
strongly  favorable  to  moderately  adverse  by  evaluating  as  described  in  §2. 1 .  This  means 
solving  for  ,  which  can  be  written  as 


yK 


(2.1) 


In  (2.1),  A,  can  be  obtained  by  solving  the  equations 
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where  k  =  0.41,  Reg  =  (0.057xURex*  *^)/v  and  the  polynomial  shown  in  (2.3)  is  a  fit 
to  a  large  number  of  turbulenee  boundary  layer  data  done  by  Das  (1987). 

Aeeordingly,  Fig.  2.4  shows  the  results  obtained  for  the  required  point  mieroaetuator  spae- 
ing  and  frequeney  response  at  four  downstream  loeations  (x/Z  =  0.25,  0.5,  0.75,  and  1.0) 
and  for  four  different  pressure  gradients.  Vehieles  for  whieh  results  are  shown  in  this  fig¬ 
ure  are  summarized  in  Table  2.1.  It  is  apparent,  both  from  the  results  shown  and  from  fun¬ 
damental  eonsiderations,  that  vehiele  speed  is  the  prineipal  faetor  that  drives  both  the 
aetuator  spaeing  and  frequeney  requirements.  These  requirements  are  relatively  benign 
for  the  UAV  application,  where  actuator  spacings  of  several  millimeters  and  step  response 
of  100  Hz  are  sufficient  to  act  on  virtually  every  sublayer  streak.  At  the  other  extreme,  the 
supersonic  fighter  and  transport  aircraft  require  actuator  spacings  of  100-200  qm  and  step 
response  of  10-90  kHz.  The  various  hydronautical  applications,  which  differ  from  the 
aeronautical  applications  due  to  the  0(10  )  lower  kinematic  viscosity,  require  microactu¬ 
ator  spacings  roughly  300  qm  but  frequencies  of  only  about  1  kHz.  Figure  2.5  gives  the 
equivalent  DC  volume  flow  rate  which  the  microactuators  must  be  able  to  achieve  for 
each  of  these  vehicle  types  in  order  to  impulsively  displace  the  required  volume  of  fluid 
during  the  available  sublayer  vortex  passage  time.  Of  key  relevance  is  that  the  flow  rate 
requirements  for  hydronautical  applications  are  typically  one  to  two  orders  of  magnitude 
lower  than  for  aeronautical  applications.  For  essentially  all  the  hydronautical  applica¬ 
tions,  the  equivalent  steady  flow  rate  required  is  in  the  range  of  10  qL/min. 

For  electrokinetic  microactuators  that  are  intended  for  use  in  dense  arrays,  the  maximum 
flow  rate  achievable  by  each  actuator  is  limited  by  the  area  per  actuator  available  for  the 
electrokinetic  pumping  channels,  as  will  be  seen  in  §4.  As  a  consequence,  the  flow  rates 
in  Fig.  2.5  are  closely  connected  with  the  mieroaetuator  spacings  in  Fig.  2.4.  This  is 
shown  in  Fig.  2.6,  where  the  maximum  flow  rate  achievable  is  shown  as  a  function  of  the 
spacing  between  individual  microactuators  within  the  array.  The  dashed  line  shows  the 
result  for  simple  geometric  scaling  of  a  given  actuator  design,  for  which  the  flow  rate  Q 
increases  with  actuator  spacing  I  as  Q  ~  The  practical  implication  of  this  is  that  a 
mieroaetuator  array  designed  for  use  in  full-scale  hydronautical  applications  can  simply 
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be  geometrically  scaled-up  in  size  by  a  factor  of  15  to  provide  the  required  microactuator 
array  for  a  UAV-scale  aeronautical  application.  The  other  aeronautical  applications,  how¬ 
ever,  would  require  microactuator  arrays  that  can  achieve  higher  electrokinetic  pumping 
efficiencies,  namely  higher  values  of  the  equivalent  DC  flow  rate  Q  per  unit  microactua- 
tor  area  /  .  For  this  reason  the  MEKA-5  microactuator  array  developed  in  this  study  (see 
§5)  has  been  designed  for  full-scale  hydronautical  applications,  with  25,600  individual 
electrokinetic  microactuators  arranged  on  325  pm  center-to-center  spacings  in  a  7  x  7  cm 
array.  The  same  array  can  be  directly  scaled  up  -  with  no  increase  in  per-actuator  perfor¬ 
mance  required  -  for  the  UAV-scale  application. 

These  frequency  and  flow  rate  requirements  represent  what  is  needed  to  act  on  essentially 
every  streamwise  vortical  sublayer  structure  in  the  turbulent  boundary  layer.  While  the 
actuator  spacing  requirements  are  inflexible,  adequate  control  may  be  possible  with  some¬ 
what  lower  frequency  performance.  This  would  involve  acting  on  some,  but  not  all,  of  the 
sublayer  vortical  structures.  In  fact,  it  is  only  necessary  to  act  on  those  structures  that  are 
at  the  point  of  incipient  bursting,  since  it  is  this  vortex  bursting  phenomenon  that  the 
microactuators  are  intended  to  suppress. 
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Vehicle  class 

Length  (m) 

Speed  (m/s) 

K  (Fm) 

/x(kHz) 

Transport  aircraft  (wing) 

3 

200 

1. 8-3.1 

15-41 

Supersonic  fighter  (wing) 

5 

300 

1.3-2.1 

30-80 

UAV  (wing) 

1 

30 

25-54 

0.05-0.22 

Naval  surface  ship  (hull) 

180 

30 

3.6-5.0 

0.41-0.84 

Submarine  (hull) 

170 

10 

2.9-4.2 

0.5-1.2 

Torpedo 

5 

10 

2.7-4.5 

0.5-1.4 

Automotive  (body) 

2 

30 

10-19 

0.4-1.4 

Table  2.1.  Full-scale  aeronautical  and  hydronautical  vehicles  types  considered  in 
§2.4,  giving  overall  vehicle  length  L  and  speed  U  used  to  obtain  viscous  sublayer  char¬ 
acteristics  for  various  values  of  the  pressure  gradient  parameter  Ft  at  locations  x/T  = 
0.25,  0.50,  0.75,  and  1.00.  The  resulting  ^  values  determine  the  required  performance 
characteristics  shown  in  Figs.  2. 4-2. 6  for  active  sublayer  control  on  each  vehicle  type  at 
full-scale  operating  conditions. 
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(a) 


Figure  2.1.  Typical  mean  velocity  profiles  in  turbulent  boundary  layers  for  dif¬ 
ferent  values  of  the  pressure  gradient  parameter  IT,  showing  the  resulting  profile 
shape  in  (a)  outer-law  variables,  and  (b)  inner-law  variables. 
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(a) 


(b) 


Figure  2.2.  Schematics  indicating  the  structure  of  streamwise  vortices  in  the 
near-wall  region  of  turbulent  boundary  layers,  showing  (a)  typical  pattern  of  high- 
and  low-speed  streaks,  and  (b)  idealized  representation  of  the  near-wall  streamwise 
vortical  structures. 
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Figure  2.3.  Idealized  notional  representation  showing  the  lateral  displaeement  of 
streamwise  sublayer  vortices  induced  by  volumetric  pumping  from  an  actuator.  Note 
that  the  Reynolds  number  Re  associated  with  the  flow  issuing  from  the  actuator  must  be 
0(1)  or  smaller,  so  that  the  actuator  flow  serves  as  a  point  volume  source  to  provide 
effective  lateral  displacement.  At  larger  Re  values,  the  flow  issuing  from  the  actuator 
would  fundamentally  change  to  become  jet-ltke,  and  would  then  no  longer  induce 
effective  lateral  displacements. 
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Figure  2.4.  Microactuator  spacing  and  frequency  requirements  for  active  sublayer 
control  by  manipulation  of  streamwise  vortical  structures  as  indicated  in  Fig.  2.3, 
showing  results  for  aeronautical  and  hydronautical  vehicle  types  listed  in  Table  2.1. 
Note  that  hydronautical  applications  all  require  microactuator  arrays  with  approxi¬ 
mately  300  pm  actuator  spacings  and  1  kHz  frequency  response. 
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Figure  2.5.  Equivalent  steady  volumetric  flow  rates  required  for  active  sublayer 
control  by  manipulation  of  streamwise  vortical  structures  as  indicated  in  Fig.  2.3, 
showing  results  for  aeronautical  and  hydronautical  vehicle  types  listed  in  Table  2.1. 
Hydronautical  applications  only  require  flow  rates  in  the  range  of  10  pL/min,  while 
aeronautical  vehicle  applications  typically  require  one  to  two  orders  of  magnitude 
larger  flow  rates. 
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Figure  2.6.  Combined  results  from  Figs.  2.5  and  2.6,  with  dashed  line  showing  cou¬ 
pling  between  microactuator  spacing  and  flow  rate,  as  noted  in  §2.4,  corresponding  to 
simple  geometric  scaling  of  any  given  microactuator  design.  Increased  actuator  spac- 
ings  provide  increased  area  for  electrokinetic  pumping  channels,  which  in  turn  provide 
for  larger  volume  flow  rates. 
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CHAPTER  3 


SYSTEM  ARCHITECTURE  FOR  MICROACTUATOR  ARRAYS 


The  small  actuator  spacing  required  for  sublayer  control  on  real  vehicles,  as  seen  in  §2, 
implies  that  large  and  dense  arrays  of  microactuators  must  be  used  to  cover  key  parts  of 
the  vehicle  surface,  as  indicated  in  Fig.  3.1.  However,  the  inherently  local  nature  of  the 
sublayer  vortex  dynamics  and  bursting  process  suggests  that  such  large  arrays  can  be  com¬ 
posed  of  much  smaller  independent  “unit  cells”,  each  with  their  own  sensors,  control  pro¬ 
cessing,  and  actuators.  This  highly  local  nature  of  the  problem  thus  permits  a  system 
architecture  for  such  microactuator  arrays  that  greatly  simplifies  the  sublayer  control  of 
turbulent  boundary  layers.  This  section  describes  a  system  architecture,  developed  as  part 
of  this  study,  based  on  this  unit-cell  concept  that  is  sufficiently  simple  to  realistically  per¬ 
mit  integration  of  the  present  electrokinetic  microactuators  into  a  complete  system  for 
active  sublayer  control  on  full-scale  vehicles  under  realistic  operating  conditions. 


3.1  Unit  Cell  Concept 

The  typically  300  pm  spacing  between  microactuators  required  for  hydronautical  vehicles 
of  100  m  scale,  such  as  a  submarine  hull,  implies  that  of  the  order  of  30  billion  individual 
microactuators  would  be  required  to  cover  the  entire  vehicle  surface.  The  massive  diffi¬ 
culty  of  integrating  such  a  large  number  of  actuators,  with  a  comparably  large  number  of 
sensors  and  an  appropriate  control  processing  capability,  can  be  greatly  simplified  by  tak¬ 
ing  advantage  of  the  inherently  local  nature  of  the  sublayer  structure  and  dynamics.  In 
particular,  the  length  of  the  sublayer  vortices  typically  extends  over  only  about  1000  wall 
units  ,  and  the  bursting  process  occurs  between  adjacent  counter-rotating  pairs  of  these 
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vortices  that  are  typieally  separated  along  the  spanwise  direction  by  only  about  100  wall 
units  4  and  are  typically  located  about  10  wall  units  4  above  the  wall  as  noted  in  §2. 
Moreover,  the  bursting  process  itself  is  principally  dependent  on  the  separation  between 
the  pair  of  eounter-rotating  sublayer  vortices,  and  is  largely  independent  of  the  locations 
and  dynamics  of  other  distantly-loeated  sublayer  vortices.  This  inherent  locality  suggests 
that  sublayer  eontrol  of  the  turbulent  boundary  layer  over  the  entire  vehicle  can  be  broken 
down  into  elementary  “unit  cells”,  the  size  of  which  is  set  by  these  length  scales  assoeiated 
with  the  bursting  proeess.  Moreover,  this  also  suggests  that  dynamieal  interactions 
between  sueh  unit  eells  should  be  comparatively  weak  in  relation  to  the  interactions 
between  sublayer  vortiees  within  a  given  unit  cell.  Any  such  interactions  would  extend  at 
most  to  the  next  adjacent  unit  eells,  and  to  a  good  approximation  sueh  eell-to-cell  interae- 
tions  could  be  neglected  altogether.  The  unit  eells  then  beeome  fully  independent,  each 
with  their  own  sensors,  control  processing,  and  actuators.  Since  the  number  of  sensors 
and  actuators  within  such  unit  cells  is  then  relatively  small,  the  problem  of  integrating 
them  into  a  eomplete  control  system  beeomes  eorrespondingly  simpler. 

The  extent  of  this  simplification  depends  on  the  minimum  length  scale  to  whieh  these  unit 
eells  ean  be  reduced.  The  1000  4  average  length  of  the  near-wall  structures  suggests  that 
this  is  the  maximum  dimension  along  the  streamwise  direction  that  could  be  required  for 
sueh  unit  eells.  Further  reduction  is  possible  by  recognizing  that  incipient  bursting  oeeurs 
far  more  locally  than  this,  over  a  smaller  length  of  the  sublayer  vortex  pair.  A  reasonable 
upper  limit  on  the  required  streamwise  dimension  of  the  unit  cells  would  thus  be  about 
500  4  to  allow  detection  of  the  incipient  bursting  configuration  and  subsequent  actuation 
to  interfere  with  this  proeess. 

In  the  spanwise  direction,  the  coherent  structures  are  separated,  on  average,  by  about  100 
wall  units  4.  Since  sublayer  control  involves  acting  on  pairs  of  such  sublayer  vortices  near 
incipient  bursting,  and  since  dynamical  interactions  via  Biot-Savart  induction  are  signifi¬ 
cant  only  between  adjacent  vortices,  it  appears  reasonable  that  the  spanwise  dimension  of 
the  unit  cell  should  be  large  enough  to  aceommodate  a  given  vortex  pair  and  its  next  adja¬ 
cent  spanwise  neighbors.  This  would  suggest  that  a  spanwise  dimension  of  500  /„  equal 
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to  the  streamwise  dimension  of  the  unit  eell,  would  be  suflheient.  Unit  eells  mueh  larger 
than  this  would  not  appear  to  be  neeessary,  and  would  simply  inerease  the  number  of  sen¬ 
sors  and  aetuators  that  must  be  integrated  into  the  loeal  unit-eell  eontrol  system. 

Based  on  this  reasoning,  the  fundamental  size  of  the  unit  eells  in  this  system  arehiteeture  is 
500  4  in  both  the  streamwise  and  spanwise  direetions.  Moreover,  the  500  4  x  500  4  unit 
eells  in  this  arehiteeture  are  taken  to  be  eompletely  independent  of  eaeh  other,  though 
eell-to-eell  interaetions  between  adjaeent  unit  eells  eould  be  readily  ineorporated  by  hav¬ 
ing  eaeh  of  the  eell-based  eontrol  proeessors  share  information  with  adjaeent  cells. 

Figure  3.1  schematically  shows  how  this  concept  can  be  used  to  group  the  sensors,  actua¬ 
tors,  and  processing  over  an  entire  vehicle  into  such  unit  cells,  each  composed  of  an  n  x  n 
arrangement  of  elementary  sensor  and  actuator  pairs  together  with  local  unit-cell  process¬ 
ing  capability.  Given  the  unit-cell  dimensions  noted  above  and  the  actuator  spacings 
noted  in  §2,  it  is  apparent  that  n  will  typically  be  in  the  range  of  4  <  n  <  6. 

Furthermore,  the  relatively  small  number  of  sensor-actuator  pairs  in  each  such  unit  cell 
allows  the  resident  processing  capability  that  “connects”  these  sensors  and  actuators  to  be 
greatly  simplified.  While  most  previous  work  has  examined  comparatively  elaborate  con¬ 
trol  processing  approaches  based  on  fundamental  control  theory,  many  of  which  require 
comparatively  sophisticated  local  processing  capability  to  implement,  the  unit-cell 
approach  permits  investigation  of  potentially  far  simpler  methods.  In  particular,  the  num¬ 
ber  of  sensors  and  actuators  is  sufficiently  small  that  true  “processing”  at  the  unit-cell 
level  may  be  discarded  entirely,  and  replaced  by  a  simple  “look-up  table”  approach.  In 
this  case,  the  actuator  states  are  directly  implied  by  the  sensor  states  via  a  look-up  table 
stored  in  a  simple  programmable  logic  array  (PLA)  in  each  unit  cell.  In  this  manner,  on 
each  clock  cycle  the  new  sensor  states  are  used  to  determine  the  new  actuator  states,  thus 
eliminating  the  need  for  any  true  processing  capability.  Additional  details  of  this  look-up 
table  approach,  including  how  the  look-up  table  in  the  PLA  is  determined,  are  discussed  in 
§3.3. 
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Moreover,  given  the  typical  dimensions  of  these  unit  cells  and  the  dimension  of  the  silicon 
wafers  on  which  microfabrication  technologies,  such  as  photolithography  and  other 
MEMS  processes,  can  be  used  to  parallel-fabricate  large  numbers  of  such  unit  cells,  it  is 
natural  to  group  these  unit  cells  into  “tiles”,  as  also  shown  in  Fig.  3.1.  Each  tile  contains 
the  number  of  unit  cells  that  can  be  fit  into  a  square  array  on  a  single  wafer.  For  typical 
unit-cell  sizes  and  wafer  sizes,  this  dictates  roughly  250  x  250  unit  cells  in  each  tile.  Each 
tile  has  a  single  pair  of  electrical  leadouts,  with  an  internal  bus  structure  that  distributes 
power  to  each  of  the  unit  cells.  These  tiles  form  the  basic  elements  used  to  cover  desired 
parts  of  the  vehicle  surface  on  which  sublayer  control  is  to  be  accomplished. 

Note  that  this  novel  system  architecture  based  on  the  unit-cell  approach  and  the  consider¬ 
able  simplifications  that  it  implies  for  active  sublayer  control  systems  are  essentially  inde¬ 
pendent  of  the  type  of  actuators  or  sensors  used.  A  similar  architecture  could  be  applied  to 
many  different  types  of  active  sublayer  control  approaches,  with  similar  practical  advan¬ 
tages  for  realistically  implementing  active  sublayer  control. 


3.2  Unit  Cell  Components 

As  indicated  in  Fig.  3.2,  each  unit  cell  in  this  system  architecture  is  composed  of  a  collo¬ 
cated  array  of  wall  shear  stress  sensors,  electrokinetic  microactuators,  and  local  digital 
signal  processing  (DSP)  capability.  This  section  briefly  describes  the  key  aspects  of  the 
major  components  that  make  up  each  of  the  fully-independent  unit  cells  on  each  tile,  con¬ 
sisting  of  an  n  X  n  array  (with  typically  4  <  n  <  6)  of  sensors  and  actuators  connected  via 
the  unit-cell  control  logic. 


3.2.1  Microactuators 

Each  unit  cell  contains  an  n  x  n  array  of  electrokinetic  actuators  that  provide  the  volume 
displacements  within  the  near-wall  region  of  the  turbulent  boundary  layer  to  laterally  dis¬ 
place  the  streamwise  sublayer  vortices.  Although  the  fundamental  unit-cell  architecture 
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can  accommodate  essentially  any  type  of  actuator,  the  electrokinetic  microactuators  that 
are  the  principle  focus  of  this  study  offer  several  benefits  that  make  them  potentially  well 
suited  for  active  sublayer  control.  The  electrokinetic  principle  on  which  these  microactu¬ 
ators  function  (see  §4)  requires  no  moving  parts,  with  the  volumetric  pumping  being 
induced  by  an  impulsively  applied  electric  field.  As  a  result,  these  electrokinetic  microac¬ 
tuators  can  be  significantly  more  robust  than  many  other  aetuators  based  on  moving  parts. 
Moreover,  the  sealings  that  govern  the  electrokinetic  effect  provide  signifieant  perfor¬ 
mance  advantageous  for  actuators  fabricated  at  the  microscale.  In  partieular,  such  electro¬ 
kinetic  microactuators  can  achieve  an  extremely  high  impulse  response,  and  equivalent 
steady  flow  rates  sufficient  to  meet  the  requirements  noted  in  §2  for  active  sublayer  con¬ 
trol  under  full-scale  vehicle  operating  conditions.  Such  microactuators  can  also  be  readily 
sized  to  fit  with  the  unit-cell  architecture,  and  do  not  involve  any  parts  that  protrude  into 
the  flow. 

Furthermore,  sueh  eleetrokinetic  mieroaetuators  can  be  fabricated  with  a  comparatively 
simple  three-layer  design,  as  indicated  in  Fig.  3.3.  This  design  involves  a  top  layer  con¬ 
taining  the  individual  actuator  nozzles  and  electrodes  together  with  leadouts  to  the 
unit-cell  DSP,  a  middle  layer  containing  the  electrokinetic  driver  channels  filled  with  a 
porous  polymer  matrix  in  which  pumping  of  the  electrolyte  oecurs,  and  a  bottom  layer  that 
serves  as  a  common  electrolyte  reservoir  for  the  unit  cell  (or,  alternatively,  for  the  entire 
tile). 

Power  for  the  electrokinetie  mieroaetuators  is  distributed  throughout  the  tile  via  a  simple 
bus  structure  composed  of  parallel  horizontal  and  vertical  line  conductors  in  the  top  layer 
across  the  entire  tile,  as  indicated  in  Fig.  3.2.  The  common  electrolyte  reservoir  serves  as 
the  other  eleetrode.  The  high  level  of  redundancy  provided  by  the  top-layer  bus  structure 
ensures  continuity  of  a  conductive  path  to  all  actuators  on  the  tile  even  in  the  event  of  a 
local  damage  site  on  the  tile  or  a  local  failure  along  any  of  the  line  conductors.  A  set  of 
leadouts  run  within  the  top  layer  from  each  individual  actuator  electrode  to  a  terminal 
located  near  one  of  the  common  power  bus  lines.  On  each  clock  cycle,  the  DSP  selec¬ 
tively  eonnects  eaeh  electrode  terminal  to  the  power  bus,  with  either  positive  or  negative 
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polarity,  to  provide  for  positive  or  negative  volume  displacement  by  each  individual  actu¬ 
ator,  or  else  leaves  this  connection  open  to  provide  for  zero  displacement  by  the  actuator. 


3.2.2  Sensors 

Detection  of  streamwise  vortical  structures  at  the  upper  edge  of  the  viscous  sublayer  is 
accomplished  by  an  n  x  n  array  of  wall  shear  stress  sensors  collocated  with  each  microac¬ 
tuator  in  the  unit  cell.  The  detailed  design  and  fabrication  of  the  wall  shear  stress  sensors 
is  discussed  in  §6.  The  unit-cell  architecture  in  Fig.  3.2  uses  one  wall  shear  stress  sensor 
for  each  actuator,  located  between  adjacent  actuators. 

Many  of  the  usual  concerns  about  shear  stress  sensor  calibration  and  accuracy  can  be 
relaxed  in  this  system  approach,  since  the  role  of  each  sensor  is  not  to  measure  the  wall 
shear  stress  distribution  below  the  vortical  structures,  but  rather  simply  to  identify  the 
presence  of  a  vortical  structure  just  prior  to  bursting  with  reasonably  high  probability  of 
detection.  Thus  the  instantaneous  wall  shear  stress  sensor  outputs  can  be  compared  to  a 
local  running  average  value  obtained  from  a  simple  (e.g.,  RC)  low-pass  filter  located  in 
the  resident  DSP  on  each  unit-cell.  When  the  sensor  output  exceeds  a  preset  multiple  of 
this  running  average,  then  the  state  of  the  /-th  sensor  in  the  processing  electronics  is  set  to 
Si  =  1  (or  to  5*,  =  ±1  if  the  sensors  have  directional  capability).  Otherwise  the  input  from 
the  sensor  is  set  to  =  0.  A  simple  circuit  design  that  generates  the  sensor  state  S^  corre¬ 
sponding  to  each  wall  shear  stress  sensor  is  shown  in  Fig.  3.4.  The  sensor  states  are  shown 
schematically  in  Fig.  3.5. 

The  complete  set  of  sensor  states  {Si}  for  i  =  (1,  2,  ...,  ng),  where  is  the  number  of  sen¬ 
sors  in  each  unit  cell,  provides  the  input  to  the  unit-cell  processing  electronics.  The  run¬ 
ning  average  essentially  eliminates  the  need  for  calibrating  each  sensor,  and  eliminates 
difficulties  caused  by  changes  in  the  vehicle  speed  and  attitude  or  by  the  particular  loca¬ 
tion  of  the  sensor-actuator  pair  on  the  vehicle.  Sensor  drift  over  time  scales  significantly 
longer  than  the  averaging  time  thus  becomes  irrelevant,  and  sensor  accuracy  does  not  need 
to  be  high  since  the  sensor  output  is  thresholded  in  the  manner  described  above. 
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3.2.3  Processing 


Since  all  interactions  between  the  sensors  and  actuators  occur  at  the  unit-cell  level,  eaeh 
unit  eell  eontains  its  own  independent  proeessing  eapability.  Beeause  of  the  relatively 
small  unit-cell  size  that  the  dynamies  of  the  sublayer  bursting  proeess  allows,  it  may  be 
possible  to  greatly  reduee  the  required  sensor  eapabilities  and  unit-eell  proeessing  eapabil- 
ities.  The  role  of  the  loeal  unit-eell  proeessing  is  to  use  the  sensor  states  {S',}  on  eaeh 
eloek  eycle  to  determine  the  Ha  aetuator  states  {^4,},  where  ^4,  =  +1,  0,  or  -1  corresponds 
respeetively  to  positive  volume  flux  (blowing),  zero  volume  flux,  or  negative  volume  flux 
(suetion).  Note  that  the  aetuators  are  not  modulated;  they  are  either  on  (±1)  or  ojf  (0). 
Thus  the  voltage  of  the  top-layer  eleetrode  for  eaeh  actuator  is  set  to  where  T{g/is  the 
voltage  of  the  common  power  bus  that  runs  between  the  unit-eells.  The  proeessing  eireuit 
thus  eflfeetively  aets  as  a  three-state  bridge  between  this  power  bus  and  the  eleetrode  eon- 
taet  for  eaeh  of  the  mieroaetuators  in  the  unit-eell.  The  relationship  between  the  sensor 
state  veetor  and  the  actuator  state  veetor  is  via  a  simple  lookup  table,  deseribed  in  the  fol¬ 
lowing  seetion. 


3.3  Look-Up  Table  Control  Logic 

The  set  of  sensor  states  {S,}  implies  a  set  of  aetuator  states  {A,}.  Funetional  approaehes  of 
the  type  A  =/ (S)  may  be  appealing  in  their  generality  and  rigorous  eonneetion  to  eontrol 
theory,  but  are  likely  to  require  signifieantly  more  proeessing  eapability  than  ean  be 
readily  aeeommodated  within  each  unit-cell.  Moreover,  sueh  an  approaeh  would  overlook 
the  faet  that  the  physies  of  the  sublayer  struetures  greatly  restriet  this  funetional  depen- 
denee  and  thus  provide  opportunities  for  eonsiderable  simplifioation.  For  example,  sinee 
the  proeessing  is  purely  loeal  at  the  unit-eell  level  and  the  unit-eell  size  ean  be  made  rela¬ 
tively  small,  it  may  be  possible  to  use  a  simple  look-up  table  approaeh  {e.g. ,  a  programma¬ 
ble  logie  array)  that  eliminates  the  need  for  eomplex  proeessing,  as  indieated 
sehematieally  in  Fig.  3.6. 
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Such  a  look-up  table  approach  can  be  implemented  at  various  levels  of  approximation. 
The  simplest  is  to  make  use  of  the  nominally  streamwise  structure  of  the  sublayer  vortices, 
and  thereby  reduce  the  number  of  spatial  dimensions  in  both  the  sensor  and  actuator  arrays 
to  just  one.  This  can  be  done  by  logically  grouping  the  sensors  in  the  unit  cell  into  nC^ 
columns  along  the  streamwise  direction,  and  then  obtaining  a  single  sensor  column  state 
CSi  for  each  column  {e.g.,  by  a  majority  rule  among  the  sensor  states  {Sj}  within  the  col¬ 
umn).  Similarly,  the  actuators  are  grouped  into  streamwise  columns,  and  a  single  actu¬ 
ator  column  state  is  used  to  drive  all  the  actuators  in  that  column.  This  reduces  the 
logic  circuit  requirements  to  setting  the  nCA-element  column  state  vector  {Cf,}  from  the 
nCs-element  column  state  vector  {C^,},  and  can  be  done  by  the  programmable  logic  array 
shown  schematically  in  Fig.  3.7.  In  that  case,  the  logic  array  may  contain  a  set  of  matrix 
coefficients  that  provide  the  actuator  state  vector  for  any  given  sensor  state  vector,  as  indi¬ 
cated  in  Fig.  3.8. 

For  the  relatively  small  number  of  columns  provided  by  the  unit-cell  architecture,  the 
number  of  possible  sensor  column  vector  combinations  may  be  sufficiently  small  that  each 
of  the  appropriate  actuator  column  vectors  can  be  explicitly  determined  from  simple 
model-based  considerations.  Alternatively,  irrespective  of  whether  the  unit-cell  process¬ 
ing  is  to  be  done  in  one  or  two  spatial  dimensions,  since  the  lookup  table  remains  static 
and  is  the  same  for  all  unit-cells  it  can  be  generated  from  a  detailed  simulation  study.  This 
might  be  done  by  a  neural  net  approach  that  evolves  to  determine  a  look-up  table  which 
minimizes  the  wall  shear  stress  within  a  unit-cell  from  discrete  thresholded  sensor  inputs 
Si.  Such  an  approach  to  generating  the  look-up  table  may  be  able  to  incorporate  recent 
advances  in  applications  of  control  theory  to  wall  shear  stress  reduction  (e.g.,  Kim  2001). 
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250  X  250  actuators 


Basic  Unit  Cell 

6x6  actuators  w/DSP 


Fundamental  Sensor/Actuator  Element 

1  sensor  +  1  actuator 


Figure  3.1.  Schematic  indicating  fundamental  architecture  for  microactuator  arrays 
based  on  a  unit-cell  structure  developed  in  this  study.  Each  unit  cell  is  composed  of 
small  arrays  of  sensors  and  actuators,  with  local  unit-cell  control  processing  capability 
resident  in  each  cell,  as  indicated  in  Fig.  3.2.  Tiles  containing  large  numbers  of  such 
units  cells  can  be  produced  using  microfabrication  techniques. 


39 


100  t 


100  t 


- ► 


actuator 

/ 


sensor 

/ 


D 


Figure  3.2.  Basic  unit-cell  architecture,  showing  collocated  n  x  n  arrays  of  typically 
4  <  n  <  6  microscale  sensors  and  actuators  spaced  100  viscous  wall  units  ^  apart,  with 
space  for  unit-cell  processing  and  with  common  voltage  bus  lines  separating  adjacent 
unit  cells.  The  relatively  small  number  of  sensors  and  actuators  in  each  unit  cell  per¬ 
mits  simple  look-up  table  approaches  for  the  unit-cell  processing. 
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Figure  3.3.  Fundamental  three-layer  design  of  electrokinetic  microactuator  arrays  at 
the  unit-cell  level,  showing  a  top  layer  containing  individual  actuator  nozzles,  elec¬ 
trodes  and  leadouts,  a  center  layer  containing  the  electrokinetic  driver  channels  in 
which  electrolyte  pumping  occurs  in  response  to  an  applied  electric  field,  and  an  elec¬ 
trolyte  reservoir  and  common  electrode  in  the  bottom  layer. 
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Figure  3.4.  Simple  cireuit  design  that  generates  the  sensor  state  S,  for  eaeh  of  the 
wall  shear  stress  sensors  within  the  unit  cell  as  described  in  §3.2.2.  Sensors  outputs  are 
thresholded  relative  to  a  running  average  to  eliminate  the  need  for  individual  sensor 
calibrations,  sensitivity  to  sensor  drift,  sensitivity  to  location  on  the  vehicle  surface,  and 
effects  of  changes  in  vehicle  speed  or  attitude. 
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Figure  3.5.  Schematic  indicating  thresholding  of  wall  shear  stress  measurements 
relative  to  a  running  average  (top)  to  produce  sensor  states  (bottom),  which  in  turn 
imply  unmodulated  actuator  states  via  the  PLA  look-up-table.  When^f,  =  +1  the  cor¬ 
responding  actuator  produces  positive  volume  displacement,  negative  displacement 
when  =  -1 ,  and  no  displacement  when  A,  =  0. 
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Figure  3.6.  Representation  of  a  unit-eell  programmable  logie  array  (PLA)  showing 
sensor  state  veetor  {S',}  input  to  PLA  formed  from  all  individual  sensor  states  within  the 
unit  eell,  and  the  aetuator  state  veetor  {A^}  output  from  PLA  that  determines  eaeh  of  the 
individual  eleetrokinetie  mieroaetuator  states.  The  relatively  small  vector  sizes  permit¬ 
ted  by  the  unit  cell  architecture  may  allow  a  simple  look-up  table  in  the  PLA. 
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Figure  3.7.  Reduction  in  sensor  and  actuator  state  vectors  described  in  §3.3  due  to 
overall  one-dimensional  structure  of  sublayer  vortical  structures,  with  correspondingly 
reduced  requirements  for  the  unit-cell  programmable  logic  array  (PLA). 
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Figure  3.8.  Minimal  programmable  logic  array  (PLA)  content  in  each  unit  cell,  giv¬ 
ing  actuator  state  vector  {A^}  as  a  linear  function  of  the  sensor  state  vector  {8,}  for  a 
restricted  subset  of  possible  sensor  states.  Several  such  response  matrices  Ry,  each 
applicable  for  such  a  restricted  set  of  sensor  states  S',,  would  be  contained  within  the 
PLA  on  each  unit  cell. 
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CHAPTER  4 


ELECTROKINETIC  THEORY 

4.1  Electrokinetic  Flow 

The  electrokinetic  microactuators  that  are  the  subject  of  this  work  function  on  the  basis  of 
the  electrokinetic  effect  (e.g.,  Potter  1961;  Burgreen  &  Nakache  1964),  first  noted  by 
Reuss  (1809).  While  this  effect  inherently  operates  at  the  microscale,  it  is  widely  used  in  a 
variety  of  practical  devices  and  processes  to  produce  macroscale  effects.  Examples 
include  transport  processes  in  emulsions  and  suspensions,  as  well  as  soil  remediation  pro¬ 
cesses,  to  name  just  a  few.  The  present  study  is  the  first  attempt  to  exploit  its  potential 
advantages  as  the  basis  for  a  new  class  of  microscale  actuators  suitable  for  active  sublayer 
control  in  turbulent  boundary  layers  at  full-scale  vehicle  operating  conditions. 

The  electrokinetic  effect  occurs  at  the  interface  between  two  phases,  typically  a  solid  and  a 
liquid,  that  are  in  direct  contact  with  one  another.  In  the  case  of  the  present  electrokinetic 
microactuators,  this  interface  is  between  a  solid  polymer  and  an  electrolytic  fluid.  An 
ionic  “double-layer”  forms  naturally  at  such  an  interface,  as  discussed  in  §4.2.  Four  dif¬ 
ferent  types  of  electrokinetic  phenomena  can  occur  as  a  result  of  interaction  between  this 
ionic  double-layer  and  an  electric  field.  These  can  be  classified  on  the  basis  of  whether  it 
is  the  liquid  phase  or  the  solid  phase  that  moves,  and  whether  the  electric  field  is  imposed 
to  create  the  motion  or  is  a  result  of  relative  motion  induced  between  the  two  phases.  The 
electrokinetic  phenomena  that  result  from  the  four  possible  combinations  are  termed  elec¬ 
troosmosis,  electrophoresis,  streaming  potential,  and  sedimentation  potential.  The  latter 
two  are  produced  by  an  imposed  motion  of,  respectively,  the  liquid  or  the  solid,  each  of 
which  results  in  a  current  due  to  the  flow  of  ions  in  the  double-layer,  and  thereby  produces 
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an  associated  electric  potential.  The  former  two  are  produced  by  an  imposed  electric  field, 
which  produces  an  electromagnetic  force  on  the  ions  in  the  double-layer  and  thereby 
induces  motion  on  them.  The  ions  then  collisionally  transfer  momentum  to  induce 
motion,  respectively,  in  the  bulk  liquid  or  on  the  solid  particles. 

The  term  electroosmosis  is  used  for  the  case  where  an  imposed  electric  field  induces  a 
flow  in  the  liquid  phase,  and  it  is  this  manifestation  of  the  electrokinetic  phenomenon  that 
drives  the  present  electrokinetic  microactuators.  Such  an  electroosmotic  flow,  shown 
schematically  in  Fig.  4.1,  results  when  an  electric  field  is  applied  to  a  liquid  electrolyte  in 
contact  with  a  charged  solid  surface,  producing  a  displacement  of  the  liquid  with  respect 
to  the  stationary  surface.  In  practice,  this  stationary  surface  might  typically  be  a  thin  cap¬ 
illary  tube  or  a  fine  porous  matrix,  with  the  electrolyte  filling  the  open  volume.  In  either 
case,  it  is  the  presence  of  the  ionic  double-layer  and  its  interaction  with  the  imposed  elec¬ 
tric  field  that  is  central  to  the  electrokinetic  flow  principle  on  which  these  microactuators 
are  based,  and  thus  to  the  operating  performance  which  these  microactuator  arrays  can 
achieve. 

For  this  reason,  §4.2  reviews  the  physical  processes  involved  in  the  double-layer  forma¬ 
tion,  including  how  the  solid  surface  acquires  its  charge  and  how  this  interface  charge 
forms  the  double-layer  and  the  associated  zeta  potential.  Following  this,  §4.3  reviews  the 
theoretical  basis  of  steady  electrokinetic  pumping  that  can  be  accomplished  by  subjecting 
this  double-layer  to  a  steady  electric  field.  Finally,  §4.4  then  develops  simple  hydrody¬ 
namic  models,  developed  as  part  of  this  study,  which  describe  the  frequency  response  of 
electrokinetic  pumping  by  an  unsteady  applied  electric  field,  in  order  to  determine  the 
impulse  response  of  electrokinetic  microactuators. 


4.2  Double-Layer  Formation 

There  are  four  fundamental  physico-chemical  mechanisms  by  which  a  solid  surface  can 
acquire  an  “interface  charge”  through  contact  with  a  liquid  electrolyte;  these  have  been 
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widely  studied  (e.g.,  Hunter  1981;  Crow  1988;  Boekris  &  Khan  1993)  and  form  the  basis 
of  the  ionic  double-layer  origin.  This  interface  charge  can  arise  from  fij  the  adsorption  of 
ions,  (ii)  the  loss  or  gain  of  lattice  ions,  (in)  the  absorption  of  ionic  surfactants,  or  (iv)  the 
ionization  of  surface  acid  or  basic  groups.  Perhaps  the  most  important  for  biological  and 
polymer  systems  is  the  latter,  which  occurs  as  a  result  of  chemical  interactions  between 
surface  acid  or  basic  groups  in  contact  with  the  electrolyte  solution.  This  is  the  case,  for 
example,  in  proteins  with  the  ionization  of  carboxyl  and  amino  groups  to  respectively  pro¬ 
duce  COO^  and  NH3'^  ions,  or  in  mineral  substances  such  as  clay  with  the  ionization  of  sil- 
anol  (Si-OH)  groups,  as  shown  in  Fig.  4.2.  In  this  study,  first-generation  electrokinetic 
actuators  were  fabricated  from  silica  capillaries  having  100  pm  inner  diameter,  and  filled 
with  1-5  pm  diameter  glass  beads  to  form  a  porous  solid  structure,  as  shown  schematically 
in  Fig.  4.3a.  The  surfaces  of  the  glass  beads  were  initially  covered  with  Si-OH  groups, 
which  ionizes  when  brought  in  direct  contact  with  a  liquid  electrolyte  to  produce  SiO^  sur¬ 
face  ions  when  the  electrolyte  pH  >2,  or  to  produce  SiOH2^  surface  ions  when  the  electro¬ 
lyte  pH  <  2.  The  current-generation  electrokinetic  microactuators  that  have  resulted  from 
this  study  work  on  a  similar  principle,  but  are  fabricated  from  a  porous  polymer  matrix 
structure  (Peters  et  al  \99%a,b)  with  average  pore  size  around  1  pm,  shown  schematically 
in  Fig.  4.3h.  In  this  case  the  polymer  is  doped  with  sulfonic  acid  groups,  which  have  the 
general  formula  R-SO4H,  where  R  can  represent  any  of  a  number  of  different  organic 
structures.  The  ionization  of  this  acid  that  results  from  contact  with  an  electrolyte,  shown 
in  Fig.  4.2a,  produces  surface  sulfate  ions  with  the  general  formula  R-SO4  . 

Ionization  of  the  surface  acid  or  base  groups  leaves  an  excess  negative  charge  on  the  sur¬ 
face  and  an  excess  positive  charge  in  the  liquid.  The  negative  charges  on  the  wall  are  bal¬ 
anced  by  positive  counter-ions  in  the  electrolyte  buffer  solution,  with  the  counter-ion 
concentration  decreasing  exponentially  with  increasing  distance  from  the  wall  in  accor¬ 
dance  with  a  Boltzman  distribution,  as  will  be  seen  in  §4.3.  At  the  same  time  the  negative 
co-ion  concentration  decreases  as  the  wall  is  approached,  since  these  ions  in  the  electro¬ 
lyte  buffer  are  being  repelled  by  the  negative  wall  charge.  Sufficiently  far  from  the  wall, 
the  counter-ion  and  co-ion  concentrations  are  the  same.  The  interaction  between  the  wall 
charges  and  the  oppositely-charged  counter-ions  in  the  electrolyte  buffer  gives  rise  to  a 
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so-called  “double-layer”  strueture  near  the  wall,  shown  sehematioally  in  Fig.  4.4. 

This  double-layer  structure  of  the  eounter-ions  in  the  vieinity  of  the  wall  is  eomprised  of 
an  immobile  “inner”  layer  and  a  diffuse  “outer”  layer.  The  inner  layer  results  from 
eounter-ions  that  are  eleetrostatieally  held  by  the  oppositely-signed  surfaee  eharge.  These 
are  solvated  ions  -  ions  surrounded  by  a  layer  of  solvent  (in  most  cases  water  molecules)  - 
that  are  in  very  elose  proximity  to  the  charged  surfaee  and  remain  immobile  with  respect 
to  it,  forming  what  is  sometimes  ealled  a  “Stern”  layer.  The  strong  eleetrostatie  attraetion 
between  the  wall  eharges  and  the  eounter-ions  in  this  inner  (Stem)  layer  prevents  these 
ions  from  moving  under  the  influenee  of  an  applied  electrie  field.  However,  further  away 
from  the  wall  the  eleetrostatie  attraction  of  counter-ions  by  the  wall  charge  is  no  longer 
suffieient  to  overeome  the  eleetrostatie  repulsion  between  these  same-signed  eounter-ions. 
As  a  result,  the  remaining  eounter-ions  required  to  balanee  the  wall  eharge  will  be  kept  in 
relatively  close  proximity  to  the  wall  due  to  their  eleetrostatie  attraction  by  the  wall 
eharge,  but  will  form  a  relatively  diffuse  layer  due  to  their  mutual  eleetrostatie  repulsion. 
The  eounter-ions  within  this  diffuse  outer  layer  (sometimes  referred  to  as  the  “Gouy” 
layer)  are  free  to  drift  under  the  influenee  of  an  applied  eleetrie  field,  and  it  is  their  motion 
in  the  presence  of  an  applied  field  that  forms  the  basis  of  the  eleetrokinetieally  induced 
flow.  For  the  eharge  states  considered  here,  the  eounter-ions  will  migrate  toward  the  cath¬ 
ode.  The  boundary  between  these  inner  and  outer  layers  is  often  referred  to  as  the  shear 
plane,  sinee  this  “slip  plane”  separates  the  immobile  inner-layer  ions  and  the  diffuse 
outer-layer  ions. 

The  detailed  stmeture  of  the  charge  distribution  within  the  outer  layer  ean  be  obtained 
from  fundamental  analyses  of  the  physieal  proeesses  involved;  this  will  be  examined  in 
§4.3.  The  resulting  separation  of  eharges  between  the  wall  and  double-layer  gives  rise  to  a 
deereasing  eleetrie  potential  with  inereasing  distanee  from  the  wall.  The  potential  differ- 
enee  between  the  shear  plane  and  the  wall  is  termed  the  “zeta  potential”,  and  provides  a 
measure  of  the  total  driving  potential  (per  unit  applied  field  strength)  available  for  eleetro- 
kinetie  pumping.  In  general,  solid-liquid  pairs  that  produee  larger  zeta  potentials  are  more 
eflfeetive  for  electrokinetie  pumping. 
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4.3  Electrokinetic  Pumping 


The  Navier-Stokes  equation  ean  be  used  to  obtain  the  detailed  outer-layer  structure  and  the 
resulting  electrokinetic  pumping  that  it  induces  under  the  influence  of  an  applied  electric 
fleld  {e.g.,  Burgreen  &  Nakache  1964;  Rice  &  Whitehead  1965).  In  its  general  form  for 
constant  viscosity,  this  equation  can  be  written  as 

P^=  -Vp  +  [iV^V-F.  (4.1) 

Under  steady-state  conditions,  for  the  axisymmetric  flow  that  results  in  an  infinitely  long 
capillary  or  pore  channel  this  can  be  written  in  polar  coordinates  as 


rdr^  dry 


-P.-F, 


(4.2) 


where  is  the  pressure  gradient  -dp/dz  and  F^  represents  the  electric  body  force  generated 
by  an  axial  electric  field  on  the  net  charge  density  p(r)  in  the  double-layer.  This  force 
can  be  written  as 


F^  =  E^p{r).  (4.3) 

The  charge  density  p(r)  at  a  distance  r  from  the  surface  can  be  written  using  the  Boltzman 
distribution,  which  gives  the  number  n,-  of  ions  of  a  particular  kind  i  which  are  at  a  poten¬ 
tial  \|/  above  that  of  a  solution  containing  uniformly  distributed  ions,  namely 

_ _ _  _ _ _  /V  -ez.\i^/{kT) 

p(r)  =  ^ripzi  =  e^UiZiC  .  (4.4) 

i  i 

/V  -eZjy/kT 

where  n,  =  n,e  .  For  a  simple  electrolyte  which  ionizes  into  equally  charged  ions  of 
valence  z,  this  charge  density  becomes 

p(r)  =  (4.5) 

which  can  be  simplified  to 


p(r)  =  -2nezsinh^. 
^  kT 


(4.6) 


This  can  be  further  simplified  for  small  values  of  FP/kT  using  the  approximation 


sinh 


e\\f\  ^  ey 
kPJ  *  kT' 


(4.7) 


This  approximation  is  in  good  agreement  with  experimental  results  up  to  values  of 
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v|/  «  50m  V  as  discussed  by  Overbeek  &  Wiersema  (1967),  Levine  et  al  (1974),  and  Bock- 
ris  &  Khan  (1993). 


Using  (4.7),  (4.6)  can  be  written  for  an  eleetrolyte  with  monovalent  ions  as 

p{r)  =  -2he\/{kT).  (4.8) 

Equation  (4.8)  can  be  solved  with  the  aid  of  the  Poisson  equation,  whieh  relates  the  poten¬ 
tial  \)/  (r)  to  the  eharge  density  p(r)  as 


1  d  _  p(r) 


rdr\  dr 


(4.9) 


where  8  =  8^  8^.  Here  8^  is  the  permittivity  of  vaeuum  and  8^  is  the  relative  permittivity  of 
the  electrolyte.  Substituting  (4.8)  into  (4.9)  gives  the  Poisson-Boltzman  equation  in  cylin¬ 
drical  coordinates  as 


1  d  fd\]f 


kV, 


(4.10) 


rdr\  df' 

2 /V  -1/2 

where  K  =  (ekT/(2en))  is  the  double-layer  thiekness.  In  kinetic  theory,  k  is  known 
as  the  Debye  length.  Typieal  values  of  the  Debye  length  for  aqueous  eleetrolyte  solutions 
can  range  from  0.5  to  1000  nanometers  (van  Oss  1994)  and  deereases  with  increasing  ion 
concentration  n  in  the  electrolyte  as  k  ~  n  ,  suggesting  that  the  double-layer  thickness 
can  be  made  be  exeeedingly  small. 


The  solution  to  (4.10)  is 


V  =  c 


koi^r) 
loiKw)  ’ 


(4.11) 


where  w  is  the  radius  of  the  eapillary,  \|/  =  ^  at  r  =  w  and  is  the  zero-order  modified 
Bessel  function  of  the  first  kind.  Substitution  of  (4. 1 1)  into  (4.8)  then  gives 


p(r)  =  -8k"C 


/o(kw)  ■ 


(4.12) 


with  this  known  charge  density,  the  equation  of  motion  in  (4.2)  can  then  be  rewritten  as 


rdr^  dry 


(4.13) 
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For  the  present  eleetrokinetie  microactuators,  no  pressure  gradient  is  applied,  and  the  flow 
instead  results  entirely  from  the  applied  electric  field  E^.  Solving  (4.13)  under  these  con¬ 
ditions,  the  velocity  is  obtained  as 


v,(r) 


hi^wy  ■ 


(4.14) 


Except  in  the  double-layer  region  r  «  w,  (4.14)  can  be  simplified  to 

v,(r)  =  (4.15) 

The  ionic  mobility  Q  is  defined  from  (4.15)  as 

Q  =  ^  .  (4.16) 

The  volume  flow  rate  Q  can  be  obtained  by  integrating  (4.14)  over  the  cross-sectional  area 
of  the  capillary  or  pore  channel  to  give 


Q  =  nw  1  - 


hi^r) 


(4.17) 


7o(kw) 

which  simplifies  for  capillaries  with  radius  much  larger  than  the  double-layer  thickness  to 

Q  =  nw^Q.E^.  (4.18) 

The  pressure  achieved  in  a  plugged  capillary  can  be  obtained  by  first  solving  (4.13)  as 


7o(Kr)A 


(4.19) 


and  then  integrating  to  obtain  the  flow  rate  Q  similarly  to  (4.17).  For  such  a  plugged  cap¬ 
illary,  the  net  flow  rate  is  zero.  Using  this  condition,  the  solution  can  be  obtained  for  the 
pressure  gradient  as 


P.  =  8pQ  1 


y{Kr)'\E, 

/„(kw)7  ■ 


(4.20) 


Again,  for  kw  >  1,  (4.20)  can  be  simplified  to 

P. 


w 


and  hence  the  force  E  achieved  by  the  actuator  channel  when  plugged  is 

E  =  8kij.QE^L  , 


(4.21) 


(4.22) 
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where  L  is  the  length  of  the  ehannel. 


'y 

These  sealings  are  for  a  single  channel.  For  a  channel  of  radius  R  consisting  of  {R/w) 
individual  pores  each  of  radius  w,  the  flow  rate  in  (4.22)  becomes 

Q  =  nR^QE,  (4.23) 

'y 

with  the  resulting  flow  speed  U=  Q/nR  ,  and  the  force  F  achieved  in  a  plugged  actuator 
becomes 

F  =  E,L.  (4.24) 

Benchtop  testing  conducted  with  electrokinetic  microactuators  formed  from  capillaries 
packed  with  glass  beads  of  size  w  and  operated  under  steady  applied  fields  (§5.1)  verified 
that  the  pressure  increased  linearly  with  applied  voltage  AV and  varied  with  bead  diameter 
(and  hence  effective  channel  radius)  as  w'^,  as  suggested  by  (4.21).  The  volume  flow  rate 
Q  also  was  observed  to  increase  linearly  with  applied  voltage  and  the  ionic  mobility  fl  as 
in  (4.23),  verifying  the  applicability  of  these  scalings  to  the  present  electrokinetic  micro¬ 
actuators. 

Note  in  (4.23)  that  the  flow  rate  achieved  is  independent  of  the  pore  radius  w,  but  the  force 
in  (4.24)  increases  as  the  pores  are  made  smaller.  This  suggests  that  by  fabricating  elec¬ 
trokinetic  channels  with  sufficiently  small  pores,  it  may  be  possible  to  achieve  flow  rates 
adequate  to  meet  the  requirements  in  §2.3  while  at  the  same  time  achieving  sufficiently 
high  pressures  in  any  plugged  actuators  to  allow  these  to  unplug  themselves.  Moreover,  as 
will  be  seen  in  the  following  section,  such  electrokinetic  actuators  fabricated  at  the 
microscale  can  have  extremely  high  frequency  response  owing  to  the  small  time  scales 
required  for  collisional  (diffusion)  transport  of  the  induced  flow  across  the  channel. 


4.4  Hydrodynamic  Models  of  Frequency  Response 

As  noted  in  Fig.  2.5,  practical  microactuators  for  sublayer  control  must  operate  at  frequen- 
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cies  around  1  kHz  for  hydronautical  applications,  and  at  frequencies  as  high  as  90  kHz  for 
the  most  demanding  aeronautieal  applieations.  For  eleetrokinetie  driver  sections  based  on 
porous  channels  with  extremely  fine  pores,  the  resulting  pore  Reynolds  number  UwN  will 
be  very  low,  and  thus  inertial  damping  of  the  unsteady  response  should  be  negligible.  Fur¬ 
thermore,  the  extremely  short  time  seale  on  whieh  the  Coulomb  foree  equilibrates  with  the 
ion  drag  suggests  that,  in  the  presenee  of  a  time-varying  applied  field,  the  frequency 
response  of  mieroaetuators  based  on  the  eleetrokinetie  prineiple  will  be  very  high.  The 
ultimate  frequency  response  limits  of  sueh  unsteady  eleetrokinetie  mieroaetuators  have 
not  been  fundamentally  examined  in  the  literature,  however  assuming  that  eapaeitanees 
are  properly  managed  there  are  two  distinctly  different  physieal  proeesses  that  can  limit 
their  performance.  The  first  is  frequeney  response  limitations  in  the  eoupling  between  the 
eleetrie  field  and  the  ion  drift  in  the  outer  layer,  however  owing  to  the  physies  assoeiated 
with  this  proeess  it  may  be  antieipated  that  this  response  limit  will  be  eonsiderably  higher 
than  the  requirements  in  Fig.  2.5.  The  second  limit  is  due  to  inertial  damping  by  the  flow 
within  the  eleetrokinetie  driver  matrix.  In  order  to  estimate  the  theoretieal  frequeney 
response  attainable  with  sueh  eleetrokinetie  mieroaetuators,  this  seetion  develops  simple 
models  based  on  the  hydrodynamies  of  eleotrokinetieally  driven  flow  to  assess  the  poten¬ 
tial  frequeney  response  and  the  key  parameters  that  set  the  response  limit. 

The  response  of  the  eleetrokinetie  flow  in  a  typieal  pore  of  radius  w  to  an  unsteady 
applied  field  ean  be  approximated  by  simple  hydrodynamie  models  in  two  elosely  related 
limits.  One  limit  applies  when  the  eleetrokinetie  double-layer  thiekness  is  mueh  smaller 
than  the  pore  radius  w,  as  is  typieally  the  ease  in  praetical  applieations.  The  motion 
indueed  in  the  double-layer  by  either  an  impulsively  applied  eleetrie  field,  or  by  a  sinusoi¬ 
dally  oseillating  applied  field,  is  then  equivalent,  respeetively,  to  Stokes’  first  or  seeond 
problem  for  the  flow  indueed  by  viscous  diffusion  above  a  moving  wall.  The  “wall”  in 
this  ease  is  the  thin  double-layer,  shown  sohematieally  in  Fig.  4.5,  and  viseosity  aets  to  dif¬ 
fuse  the  indueed  motion  within  it  throughout  the  pore.  The  solutions  to  both  of  Stokes’ 
problems  are  eompletely  elassieal,  and  thus  attention  here  will  be  foeused  on  their  impli- 
eations  for  the  frequeney  response  of  eleetrokinetie  flows.  In  both  oases,  the  motion  is 
oonfined  to  a  “Stokes  layer”  with  thiekness  5  «  adjaoent  to  the  wall  within  which  vis- 
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cous  dilfusion  is  able  to  induce  flow. 


Stokes’  second  problem  corresponds  to  the  motion  induced  in  the  fluid  by  a  sinusoidally 
oscillating  wall,  and  is  most  directly  associated  with  the  concept  of  frequency  response. 
This  case  applies  to  the  motion  induced  in  the  bulk  electrolyte  outside  the  double-layer  by 
a  sinusoidally  time -varying  applied  electric  field,  which  produces  a  sinusoidal  motion  in 
the  thin  double-layer  directly  adjacent  to  the  wall.  The  resulting  velocity  field  in  this  case 
is  self-similar,  and  given  by 

u{y,t)  =  Ue^e~^cos{(x)t- r\)  ,  (4.25) 

where  r[=y/  J2\/ co  is  the  wall-normal  similarity  coordinate. 


Figures  4.6a,b  respectively  define  the  actual  flow  rate  Q'ac  that  results  from  such  a  sinu¬ 
soidally-induced  motion  during  each  volumetric  displacement  “stroke”  of  the  actuator, 
and  the  flow  rate  Q'dc  that  would  result  during  the  same  time  by  steady  electrokinetic 
pumping.  Inertial  effects  in  the  induced  flow  are  negligible  when  “quasi-steady”  perfor¬ 
mance  is  achieved,  namely  when  Q' ac^ Q' dc  ~  1  •  The  flow  rate  Q' ac  for  a  circular  chan¬ 
nel  or  pore  is 


+  774 


u{r,  t)2nr(dr)dt , 


which  can  be  approximated  by 


(4.26) 


Q'Ac^^2Kw8Ut^,  (4.27) 

where  is  the  average  velocity  in  the  Stokes’  layer  directly  adjacent  to  the  wall,  and  T 
is  the  period  corresponding  to  each  stroke  of  the  actuator.  Substituting  for  T  and  5  in 
(4.27)  results  in 


K  1/2  -3/2 

Qac-^w'^  ra  UeqQ. 

J2 


(4.28) 


An  electrokinetic  microactuator  of  radius  R  filled  with  a  porous  matrix  with  pore  size  w 
contains  {R/w)  pores.  From  (4.31)  the  total  flow  rate  Q'ac  produced  by  these  N 
pores  in  the  entire  actuator  becomes 
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When  a  single  eapillary  or  pore  is  driven  by  a  steady  (DC)  field,  as  shown  in  Fig.  4.6h, 
then  the  resulting  flow  rate  Q'dc^^ 

Q\c  =  (4.30) 

Similarly,  for  the  entire  aetuator  the  total  flow  rate  Q' dc  can  be  written  as 

Q  DC  =  (4.31) 

The  ratio  of  these  two  flow  rates  then  allows  the  influenee  of  inertial  effeets  to  be  evalu¬ 
ated  as 


Q'  AC 
Q'  DC 


1/2  T 

V2— 

W 

^  eq 


(4.32) 


2 

From  (4.32)  a  eut-off  frequeney  cOc  =  v/ (2iv  )  ean  be  dedueed  below  whieh  the  aetuators 
will  show  “quasi-steady”  performanee,  and  above  whieh  a  roll-off  will  oeeur  as  shown  in 
Fig.  4.7.  Assuming  an  aqueous  eleetrolyte  with  kinematie  viseosity  eomparable  to  that  of 
water  (v  «  10  ^  m/s),  this  indieates  a  eut-off  frequeney  of  the  order  of  500  kHz  for  1  qm 
pores.  For  the  porous  polymer  matrix  used  to  fabrieate  the  arrays  in  this  study,  the  average 
pore  size  is  roughly  1  qm,  and  this  would  suggest  essentially  loss-less  frequeney  response 
up  to  around  500  kHz. 


Stokes’  first  problem  eorresponds  to  the  ease  where  the  wall  is  impulsively  set  into 
motion,  as  would  apply  for  an  impulsively  applied  eleetrie  field.  The  flow  veloeity  with 
inereasing  distanee  from  the  wall  and  with  inereasing  time  again  has  a  similarity  solution 
in  terms  of  the  eomplementary  error  funetion.  The  resulting  thickness  5  of  the  layer  adja¬ 
cent  to  the  wall  within  which  viscous  diffusion  is  able  to  induce  flow  is  given  by  5  «  JVt . 
The  volume  flow  rate,  defined  as 

r’" 

Q{t)=\  u{r,t)2nrdr,  (4.33) 

•’  0 

can  be  readily  integrated  to  give 

Q{t)^TiwJvtU,^ 


(4.34) 


where  is  the  equivalent  wall  veloeity,  in  this  ease  the  veloeity  in  the  outer-layer 
indueed  by  the  electrokinetic  flow.  Defining  the  steady  flow  rate  Qj^c  =  aehieved 

by  a  steady  (DC)  applied  voltage  aeross  the  eleetrokinetie  driver  channel,  the  ratio  of  the 
unsteady  to  steady  flow  rates  is  given  by 

(4.35) 

Qdc 

Thus  a  characteristic  startup  time  x  =  {w  /v)  can  be  deduced  from  (4.35).  Assuming  an 
aqueous  electrolyte  with  kinematic  viscosity  comparable  to  that  of  water  (v  «  10 mVs) 
and  w  «  1  pm  pore  sizes,  this  would  indicate  that  the  unsteady  flow  rate  in  (4.35)  will 
remain  essentially  the  same  as  the  steady  flow  rate  up  to  applied  frequencies  of  the  order 
of  10'’  Hz.  In  other  words,  essentially  no  loss  of  electrokinetic  pumping  performance  will 
be  seen  until  the  frequency  of  the  applied  electric  field  approaches  1  MHz.  The  volumet¬ 
ric  pumping  capabilities  of  electrokinetic  microactuators  operated  at  frequencies  below 
this  can  be  directly  inferred  from  equivalent  steady  DC  flow  rates. 

The  second  limit  applies  when  the  double-layer  thickness  is  greater  than  the  pore  radius,  in 
which  case  the  electrolyte  throughout  the  pore  contains  a  largely  uniform  concentration  of 
ions.  The  electrokinetic  forcing  then  induces  motion  directly  throughout  the  fluid,  and 
thereby  acts  like  a  body  force.  The  response  is  thus  equivalent  to  the  oscillatory  flow 
induced  in  a  pipe  by  a  sinusoidally  varying  pressure  gradient.  This  limit  also  has  an  exact 
solution  (Sexl  1930)  that  can  be  used  to  assess  frequency  response  limits  applicable  to  that 
case.  However,  since  the  double-layer  thicknesses  in  the  present  electrokinetic  microactu¬ 
ators  are  much  smaller  than  the  pore  sizes,  this  limiting  case  does  not  appear  to  be  directly 
relevant.  For  actuators  fabricated  with  much  smaller  pore  sizes,  or  for  liquid-solid  pairs 
with  much  larger  Debye  thicknesses,  this  would  be  the  appropriate  limit  for  determining 
the  frequency  response. 

These  considerations  indicate  that  it  should  be  possible  to  meet  the  frequency  response 
requirements  for  all  of  the  vehicle  types  in  Fig.  2.5  with  electrokinetic  microactuators 
based  on  aqueous  electrolytes  and  porous  matrix  structures  having  typically  1  qm  pore 
radii  in  the  electrokinetic  channels. 
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Figure  4.1.  Schematic  showing  electroosmotic  flow  induced  by  an  applied  electric  field  along  a  pore  or  capillary 
channel.  The  electric  field  induces  drift  in  the  mobile  ions  within  the  diffuse  outer  layer,  which  then  collisionally  trans¬ 
fer  momentum  to  the  electrolyte.  Subsequent  collisional  transport  (diffusion)  transfers  this  motion  throughout  the  chan- 


(a) 

solid  surface 


(b) 

solid  surface 


solid-liquid  interface 


negatively  ctiarged  liquid 

surface 


solid-liquid  interface 


Figure  4.2.  Schematics  showing  charge  formation  mechanism  at  a  solid-liquid 
interface  by  ionization  of  (a)  sulfonic  acid  groups,  and  (b)  silanol  groups. 
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diffuse  layer 


bulk  fluid 


Xy  Water  molecule 


double  layer 


Counter-ions 


Figure  4.4.  Schematic  representation  of  the  double-layer  structure  involved  in 
electrokinetic  flows,  showing  wall  charges  and  immobile  inner  (Stem)  layer,  outer 
diffuse  (Gouy)  layer  composed  of  hydrated  counter-ions,  and  bulk  electrolyte  out¬ 
side  double-layer. 
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Figure  4.5.  Schematic  showing  Stokes  layer  of  thickness  5  within  which  motion 
induced  by  a  sinusoidally-applied  electric  field  is  confined.  Oscillating  “wall”  in  this 
case  results  from  oscillating  motion  in  the  exceedingly  thin  double-layer  immediately 
adjacent  to  the  actual  wall. 
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(a) 


Figure  4.6.  Schematics  defining  (a)  the  volume  Q' ac  displaced  by  electrokinetic 
pumping  per  actuator  stroke  time  774  under  a  sinusoidally-applied  electric  field,  and  (b) 
the  volume  Q' dc  displaced  during  the  same  time  774  under  a  steady  applied  electric 
field. 
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Figure  4.7.  The  ratio  of  Q' ac  to  Q' dc  ,  defined  in  Fig.  4.6,  from  (4.32)  for  any  sinuso¬ 
idal  eomponent  of  frequency  co  in  an  unsteady  applied  electric  field.  For 
O'  AC  /  Q'  Dc~  volumetric  displacement  achieved  by  electrokinetic  pumping  is  in  the 
quasi-steady  limit  with  no  losses  due  to  inertia  effects.  Note  that  for  1  pm  pore  sizes  as 
in  §5,  roll-off  produced  by  inertial  effects  does  not  begin  until  frequencies  above  500 
kHz. 
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CHAPTER  5 


ELECTROKINETIC  MICROACTUATOR  ARRAYS 


5.1  Elementary  Packed-Capillary  Actuators 

Initial  experiments  with  electrokinetic  flow  as  the  fundamental  basis  for  a  new  class  of 
microactuator  arrays  were  performed  using  simple  packed  capillaries  to  provide  the  dou¬ 
ble-layer  surface  that  drives  the  flow.  Testing  of  these  actuators  was  done  in  collaboration 
with  researchers  in  the  Microtechnology  Division  at  Sandia  National  Laboratories  (e.g., 
Dahm  et  al  1997)  and  continued  in  the  initial  phases  of  the  present  study.  These  initial 
“actuators”  consisted  of  glass  capillaries  with  100  pm  inner  diameter,  packed  with  glass 
beads  of  various  diameters  ranging  from  0.5  pm  to  5  pm,  and  fdled  with  various  aqueous 
electrolytes  that  provided  differing  zeta  potentials.  The  glass  beads  created  effective  inter¬ 
stitial  channels  at  submicron  scales.  The  purpose  of  these  experiments  was  to  verily  the 
applicability  of  the  scaling  results  in  §4  for  such  packed-capillary  electrokinetic  pumps  at 
these  scales  under  steady  (DC)  applied  fields.  The  configuration  used  for  these  tests  is 
shown  schematically  in  Fig.  5.1. 

Results  verified  the  Q  ~  scaling  in  (4.18)  of  the  flow  rate  Q  with  pore  size  w  (i.e., 
bead  diameter)  and  applied  field  strength  E^,  and  the  Ap  ~  E^LIvP  scaling  (4.21),  to  pore 
sizes  w  at  least  as  small  as  those  produced  by  the  0.5  pm  glass  beads.  In  view  of  the 
exceedingly  small  Debye  length  scales  associated  with  the  double-layer  thickness,  which 
typically  is  in  the  range  of  10-100  nm,  the  observed  applicability  of  these  scalings  even  at 
the  smallest  bead  diameter  (500  nm)  is  not  entirely  surprising.  Nevertheless,  the  results 
confirmed  that  fabrication  of  electrokinetic  microactuators  with  pore  sizes  down  to  at  least 
these  scales  does  not  invalidate  these  fundamental  performance  scaling  laws.  This  indi- 
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cates  that  electrokinetic  microactuators  based  on  porous  matrix  structures  having  average 
pore  sizes  of  1  pm,  for  which  the  entire  range  of  pore  sizes  should  be  within  the  range 
applicable  to  the  tests  above  using  0.5  pm  beads,  should  allow  large  values  of  pressure  rise 
ISp  to  be  achieved  in  plugged  actuators,  while  causing  essentially  no  reduction  in  volume 
displacement  rate  Q  for  open  actuators.  All  subsequent  micro-electro-kinetic-actuator 
(MEKA)  arrays  developed  in  this  study  were  based  on  electrokinetic  flow  produced  in 
such  micron-scale  porous  matrix  structures. 


5.2  Three  Layer  Design 

The  technology  for  creating  porous  polymer  matrix  structures,  with  1  pm  scale  pore  sizes, 
that  can  be  suitably  doped  to  provide  electrokinetic  flow  (Peters  et  al  \99%a,b)  provided 
the  basis  for  a  fundamental  three-layer  design  of  electrokinetic  microactuator  arrays.  This 
three-layer  design,  shown  schematically  in  Fig.  5.2,  uses  a  center  layer  containing  individ¬ 
ual  electrokinetic  driver  channels  formed  in  an  appropriate  substrate  material,  as  discussed 
in  the  following  sections,  and  filled  with  this  porous  polymer  matrix  structure.  The  porous 
matrix  structure  can  be  readily  created  throughout  these  individual  channels  by  introduc¬ 
ing  the  polymer  in  liquid  form  immediately  after  mixing  its  two  principal  constituents. 
Dopants  are  previously  added  to  these  principal  polymer  constituents  as  needed  to  provide 
the  desired  acidic  or  basic  groups  on  the  solid  surface;  the  particular  choice  and  concentra¬ 
tion  of  dopants  depends  on  the  intended  electrolyte  and  the  desired  zeta  potential.  Dis¬ 
solved  oxygen  is  removed  from  the  liquid  polymer  mixture  by  repeated  vacuum  pumping 
before  curing.  Subsequent  crosslinking  of  the  polymer  in  the  curing  step,  during  which 
the  liquid-filled  center  layer  is  placed  between  glass  plates  and  “baked”  at  40-60  °C  for 
several  hours,  produces  the  porous  structure  within  the  center  layer.  Subtle  changes  in  the 
polymer  constituents  and  in  the  curing  process  can  be  used  to  control  the  resulting  nomi¬ 
nal  pore  size.  The  center  layer  in  this  three-layer  design  consists  only  of  the  individual 
electrokinetic  driver  channels  filled  with  such  porous  polymer  matrix  structures,  as  shown 
in  Fig.  5.3. 
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Once  the  euring  proeess  in  the  eenter  layer  is  eompleted  to  produee  the  porous  polymer, 
the  resulting  strueture  ean  be  readily  wetted  with  eleetrolyte  by  simple  eapillary  wieking 
action.  Continual  contact  with  the  eleetrolyte,  and  thereby  eontinual  wetting  of  the  entire 
porous  matrix  throughout  eaeh  individual  aetuator  ehannel  in  the  eenter  layer,  is  provided 
by  a  eommon  eleetrolyte  reservoir  in  the  bottom  layer,  as  shown  in  Fig.  5.2.  Sinee  eaeh 
unit  eell  acts  as  a  zero-net-displacement  entity,  a  fixed  amount  of  electrolyte  is  maintained 
within  each  unit  cell  (or,  alternatively,  within  eaeh  tile,  as  diseussed  in  §3). 

Note  that  no  eleetrolyte  is  diseharged  from  the  unit  eell;  the  eleetrokinetie  flow  indueed  in 
the  electrolyte  when  an  eleetric  field  is  applied  aeross  any  individual  eenter-layer  ehannel 
serves  only  to  displaee  the  working  fluid  (ambient  fluid)  from  the  top  layer.  This  top 
layer,  shown  in  Fig.  5.4,  ineludes  a  nozzle  directly  above  eaeh  aetuator  channel,  as  seen  in 
Fig.  5Ab,  in  whieh  working  fluid  is  displaeed  by  pumping  of  eleetrolyte  in  the  eenter  layer 
and  issues  through  the  top  surfaee  shown  in  Fig.  5.4c.  In  oases  where  the  eleetrolyte  must 
be  separated  from  the  working  fluid,  as  would  typically  be  the  ease  for  hydronautioal 
applioations,  a  thin  layer  of  flexible  membrane  material  would  be  sandwiohed  between  the 
eenter  and  top  layers.  The  top  layer  also  oontains  a  oiroular  disk  eleotrode  for  eaeh  aotua- 
tor  channel,  with  leadouts  to  conneotors  near  the  periphery  of  the  unit  eell.  These  are  oon- 
neoted  with  either  positive  or  negative  polarity  by  the  unit-cell  control  circuitry  to  the 
eommon  voltage  bus  that  runs  between  unit  oells.  This  provides  for  a  positive  or  negative 
eleetrie  field  along  eaeh  aetuator  channel  between  the  top-layer  eleetrode  and  the  eommon 
electrode  in  the  electrolyte  reservoir  in  the  bottom  layer,  produeing  positive  or  negative 
pumping  of  the  electrolyte,  and  thereby  produeing  positive  or  negative  displaeement  of  the 
working  fluid  from  the  top  surfaee. 

This  basic  three-layer  design  allows  for  dense  arrays  eomprised  a  large  numbers  of  sueh 
actuators  to  be  fabrieated  in  thin  eonformally-applieable  tiles  that  aeeommodate  the 
unit-cell  system  architecture  deseribed  in  §3,  and  that  may  in  prineiple  be  easily  applied  to 
the  surfaee  of  real  vehieles.  Elements  of  this  basie  three-layer  design  have  been  developed 
in  sueeessive  generations  of  eleetrokinetie  mieroaetuator  arrays  deseribed  in  the  following 
seetions.  The  most  advaneed  of  these  arrays  (MEKA-5)  eonsists  of  a  eomplete  tile  of 
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25,600  such  individual  electrokinetic  microactuators,  spaced  on  325  pm  centers  and  orga¬ 
nized  into  a  40  X  40  array  of  unit  cells,  eaeh  with  a  4  x  4  arrangement  of  aetuators,  and 
fabrieated  from  this  fundamental  three-layer  design. 


5.3  MEKA-0  Microactuator  Array 

Based  on  the  results  from  experiments  with  elementary  packed  channel  actuators,  an  ini¬ 
tial  3x3  array  of  electrokinetic  microactuators,  shown  in  Figs.  5.5a-e,  was  fabricated  and 
used  for  performanee  testing.  This  array  is  somewhat  smaller  than  the  typieal  5x5 
unit-eell  indieated  in  Fig.  3.2,  however  its  purpose  was  principally  for  an  initial  explora¬ 
tion  of  the  major  ehallenges  involved  in  unit  eell  fabrication,  and  for  initial  performance 
testing  of  unsteady  (AC)  eleetrokinetie  flows.  For  added  simplicity,  this  initial  array  was 
sized  to  the  requirements  of  the  UAV-scale  applieation  in  Figs.  2.4  and  2.5,  and  thus  eon- 
sisted  of  individual  mieroaetuators  with  1000  pm  diameters  and  2000  pm  eenter-to-eenter 
spaeings.  The  basic  three-layer  design  in  Fig  5.2  was  used,  with  eleetrokinetie  driver 
channels  in  the  center  layer  and  an  electrolyte  reservoir  in  the  bottom  layer,  but  with  elee- 
trodes  and  leadouts  formed  on  the  top  surface  of  the  eenter  layer. 

The  eleetrokinetie  driver  ehannels  were  in  this  ease  meehanically  drilled  into  a  3  mm  thiek 
glass  substrate.  This  was  intended  as  a  preliminary  proeess  used  for  this  initial  fabrication 
only,  and  not  intended  to  represent  a  fabrieation  proeess  suitable  for  eventual  mass  pro- 
duetion.  Note  that  the  center-layer  thiekness  determines  the  length-to-diameter  ratio  of 
the  electrokinetic  channels,  whieh  should  be  kept  large  enough  {LID  >  3)  to  maintain  rea¬ 
sonably  straight  electric  field  lines  within  the  ehannel. 

An  individual  ring  eleetrode  with  a  leadout  to  an  electrical  contact  was  formed  on  the  top 
surfaee  around  the  periphery  of  eaeh  of  the  mieroaetuator  ehannels  using  standard  metalli¬ 
zation  and  photoetehing  teehniques,  as  can  be  seen  in  Figs.  5.5c-e.  Due  to  ehipping  that 
oeeurred  in  the  ehannel  drilling  process,  as  shown  in  Figs.  5.6a-d,  the  resulting  edge  qual¬ 
ity  was  relatively  poor.  This  eaused  some  diflfieulties  in  maintaining  eleetrode  continuity 
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on  the  top  surface  of  the  center  layer,  however  as  Figs.  5.6a,b  show  the  resulting  ring  elec¬ 
trodes  were  still  of  acceptable  quality.  The  bottom  surface  electrode  is  the  common  for  all 
the  channels,  and  can  be  seen  in  Figs.  5.6c, d. 

The  individual  center-layer  electrokinetic  driver  channels  in  Figs.  5.6a-d  were  filled  with 
the  porous  polymer  matrix  structure  described  in  §5.2.  For  the  MEKA-1  array  this  used  a 
polymer  consisting  of  several  constituents,  including  butyl  methacrylate  with  a  solvent 
and  with  a  dopant  to  provide  the  zeta  potential.  The  resulting  porous  polymer  matrix 
structure  is  visible  as  the  white  material  in  the  channels  in  Figs.  5.5c-e.  Filling  of  the 
porous  polymer  with  electrolyte  was  readily  accomplished  by  simply  applying  a  DC  volt¬ 
age  across  the  top  and  bottom  surface  electrodes  for  each  channel.  The  electrolyte  quickly 
wetted  the  entire  channel  and  could  be  readily  seen  to  flow  from  the  top  of  each  actuator. 

DC  performance  testing  was  conducted  using  the  apparatus  shown  in  Fig.  5.7.  This  con¬ 
sisted  of  an  8.5  mm  in  diameter  by  2.5  mm  in  height  electrolyte  reservoir  which  had  the 
MEKA-1  array  wafer  located  above  it.  The  reservoir  was  filled  with  a  2mM  sodium  borate 
buffer  solution.  When  a  steady  voltage  difference  AV  was  applied  across  the  electrodes, 
the  electrokinetic  actuators  steadily  pumped  the  electrolyte  into  the  graduated  reservoir 
shown  in  Eig.  5.7.  The  volume  flow  rate  was  measured  over  an  extended  time  and  the 
results  are  shown  in  Eig.  5.8.  Note  that  the  total  volume  displaced  increased  linearly  with 
time,  with  the  flow  rate  Qj^c  being  given  by  the  slope  of  each  set  of  symbols.  The  flow 
rate  achieved  is  seen  to  be  proportional  to  the  applied  voltage,  consistent  with  (4.23).  Eor 
these  conditions,  the  porous  polymer  matrix  material  provided  a  volume  flow  rate  Qj^c  of 
0. 1  pE/min  at  an  applied  field  strength  of  33  V/mm.  The  ratio  {Qacl^oc)  should  thus 
be  independent  of  the  field  strength,  and  represents  the  fundamental  characterization  of 
the  electrokinetic  pumping  capability  of  the  porous  polymer  matrix  structure  per  unit 
cross-sectional  area  of  the  actuator  channels. 

AC  performance  testing  was  also  conducted  to  determine  frequency  response,  by  applying 
a  sinusoidally  time-varying  voltage  difference.  These  measurements  were  performed  with 
a  1  mm  diameter  by  2  mm  thick  polymer  electrokinetic  pump,  driven  by  a  function  gener- 
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ator  into  a  Burr-Brown  3745  high-voltage  op-amp.  The  applied  voltage  was  varied  from  5 
V  peak-to-peak  up  to  90  V  peak-to-peak.  In  this  ease,  the  mieroactuator  drove  a  small 
elosed  eavity  eontaining  an  air  bubble,  and  the  time -varying  pressure  that  developed  in  the 
eavity  as  the  air  bubble  was  eompressed  by  the  time-varying  volume  displaeement 
aehieved  by  the  aetuator  was  measured.  A  pressure  transdueer  (Omega  PX170-14GV) 
was  used  at  low  frequeneies,  and  a  mierophone  (Briiel  &  Kjser  4138  mierophone  with 
B&K  2618  preamp)  at  high  frequeneies.  Results  are  shown  in  Fig.  5.9,  giving  the  relative 
peak-to-peak  variation  in  the  signal  from  the  transdueer  or  mierophone  normalized  by  the 
manufaeturer’s  response  ealibration.  The/  *  sealing  shown  by  the  solid  line  is  the  theo¬ 
retical  result  corresponding  to  loss-less  AC  performance;  i.e.,  the  reduction  in  volume  dis¬ 
placement  per  cycle  resulting  solely  from  the /  *  reduction  in  cycle  period  with  frequency. 
The  equivalent  DC  flow  rate  Q' oc  defined  in  Fig.  4.6a  is  seen  in  Fig.  5.9  to  be  essentially 
constant  up  to  applied  frequencies  at  least  as  high  as  10  kHz.  The  origins  of  the  roll-off 
above  10  kHz  remain  unclear,  and  may  be  due  to  an  experimental  limitation.  The  mea¬ 
sured  frequency  at  which  this  roll-off  begins  would  be  consistent  with  the  onset  of  inertial 
damping  in  a  porous  matrix  structure  with  10  pm  pores,  based  on  the  considerations  in 
§4.4,  whereas  the  nominal  pore  size  in  the  electrokinetic  driver  channels  is  believed  to  be 
about  1  pm.  However  it  is  also  possible  that,  while  the  pore  size  in  the  present  matrix  is 
nominally  1  pm,  the  largest  pores  may  be  as  large  as  10  pm,  and  that  these  determine  the 
inertial  damping  limit  for  the  actuator.  Irrespective  of  these  considerations,  the  10  kHz 
frequency  response  demonstrated  in  these  tests  meets  the  demands  for  most  of  the  vehicles 
in  Figs.  2.4  and  2.5,  and  far  exceeds  the  requirements  for  all  the  hydronautical  applica¬ 
tions. 


5.4  MEKA-1  Microactuator  Array 

Based  on  the  performance  seen  with  the  MEKA-0  mieroactuator  array  in  §5.3,  a  MEKA-1 
array  was  fabricated  to  (i)  assess  laser  drilling  of  the  electrokinetic  channels  as  an  alterna¬ 
tive  to  mechanical  drilling,  and  (ii)  explore  design  issues  associated  with  larger  unit  cells. 
This  10  X  10  actuator  array  was  fabricated  in  glass,  as  was  the  case  for  the  previous 
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MEKA-0  array,  and  again  sized  for  a  UAV-scale  application.  The  individual  electrokinetic 
microactuator  channels  were  2000  pm  in  diameter  and  located  on  4  mm  center-to-center 
spacings.  The  layout  for  this  MEKA-1  array  is  shown  in  Eig.  5.10. 

Drilling  of  the  electrokinetic  driver  channels  in  glass  was  in  this  case  done  with  a  CO2 
laser  by  a  commercial  vendor.  This  fabrication  process  is  amenable  to  eventual  mass  pro¬ 
duction.  A  comparison  of  the  resulting  edge  quality  produced  by  this  laser  drilling  process 
in  MEKA-1  with  that  of  the  mechanically  drilled  holes  in  MEKA-0  is  shown  in  Eigs. 
5.\\a,b.  Edge  quality  achieved  at  this  scale  with  laser  drilling  was  excellent,  and  more 
than  adequate  for  metallization  and  electrode  etching. 

The  10x10  unit-cell  size  of  this  MEKA-1  array  is  considerably  larger  than  the  6x6 
unit-cell  anticipated  in  Eig.  3.2  for  the  proposed  system  architecture,  yet  Eig.  5.10  shows 
that  the  required  electrode  leadouts  can  be  readily  accommodated.  Eike  the  MEKA-0 
array,  the  bottom  surface  electrode  is  common  for  all  channels.  The  channels  were  filled 
with  the  same  porous  polymer  matrix  as  for  the  MEKA-0  array.  Eigures  5. 12a, 6  show  the 
resulting  filled  microactuator  channels,  where  the  1  pm-scale  porous  polymer  matrix 
structure  again  appears  white.  Note  that  the  individual  electrokinetic  microactuator  driver 
channels  in  the  center  layer  can  be  readily  filled  with  the  polymer  using  the  procedure  out¬ 
lined  in  §5.2. 


5.5  MEKA-2  Microactuator  Array 

In  view  of  the  excellent  edge  quality  achieved  by  the  laser  drilling  process  in  the 
UAV-scale  MEKA-1  array,  the  process  was  applied  to  develop  a  hydronautical-scale 
MEKA-2  array.  The  principal  purpose  of  this  array  was  to  test  the  suitability  of  the  CO2 
laser  drilling  process  in  a  glass  substrate  at  much  smaller  scales.  Accordingly,  this  array 
consisted  of  a  simple  3x3  matrix  of  electrokinetic  driver  channels,  each  with  300  pm 
diameter  and  2  mm  length,  on  500  pm  center-to-center  spacings. 
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The  resulting  hole  quality  is  shown  in  Figs.  5.\3a,b,  where  both  the  top  (entry-side  for  the 
laser  beam)  and  bottom  surfaees  are  shown.  As  can  be  seen,  the  smaller  channel  diame¬ 
ters  in  this  case  produced  much  poorer  results  from  laser  drilling  than  was  the  case  for  the 
2000  pm  diameter  channels  in  the  MEKA-1  array.  These  difficulties  were  driven  by  fab¬ 
rication  glass  substrate  material.  Since  glass  was  not  intended  as  a  final  substrate  material, 
subsequent  arrays  examined  fabrication  with  other  substrate  materials. 


5.6  MEKA-3  Microactuator  Array 

The  CO2  laser  drilling  process  was  adapted  for  a  MEKA-3  array  by  changing  the  substrate 
material  from  glass  to  acrylic,  which  is  also  readily  amenable  to  laser  drilling.  The 
MEKA-3  array,  shown  in  Eigs.  5.\Aa,b,  consisted  of  a  5  x  5  matrix  of  electrokinetic  driver 
channels,  each  300  pm  in  diameter  and  2  mm  long,  and  with  1000  pm  center-to-center 
spacings. 

Eigure  5.15a  shows  close-up  views  of  the  hole  quality  resulting  on  the  exit  side  from  laser 
drilling  process  at  this  scale  in  acrylic.  Edge  quality  was  generally  excellent  and  more 
than  adequate  for  metallization.  There  was  indication  of  melt  flow  around  the  periphery 
of  the  hole,  with  further  indications  of  a  circumferential  instability  in  the  melt  flow,  but 
these  do  not  compromise  the  hole  quality  needed  for  metallization.  Eigure  5.15Z)  shows 
much  poorer  edge  quality  on  the  entry  side  of  the  acrylic.  The  degradation  appears  due  to 
cratering  produced  by  ejection  of  molten  material  from  the  hole,  as  well  as  melt  flow 
along  the  surface.  However,  since  only  one  side  of  the  center  layer  needs  to  be  metallized, 
this  can  be  done  on  the  entry  side,  where  edge  quality  is  higher. 

The  results  obtained  from  this  MEKA-3  fabrication  with  laser  drilling  in  acrylic  indicate 
that  it  is  possible  to  mass  produce  the  center  layer  in  thin  plastic  films  by  such  a  process. 
However,  acrylic  is  unsuitable  as  a  center-layer  substrate  material  since  it  is  incompatible 
with  the  organic  solvents  used  in  the  porous  polymer  process.  There  are  other  plastics  that 
can  be  readily  laser  drilled  and  that  are  compatible  with  these  solvents.  Eor  this  reason. 
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the  center  layer  for  MEKA-5  array  in  the  following  section  was  produced  in  mylar,  which 
can  be  laser-drilled  and  is  suitable  for  use  with  the  porous  polymer  process. 


5.7  MEKA-5  Hydronautical  Array 


5,7.1  MEKA-5  Array  Design 

The  MEKA-5  actuator  array,  shown  in  Eigs.  5.16,  demonstrated  fabrication  of  a  full 
hydronautical-scale  array  of  electrokinetic  microactuators  and  their  integration  with  a 
top-layer  containing  the  basic  unit-cell  structure  and  all  electrical  leadouts  required  for 
actuation.  The  array  was  fabricated  in  a  7  x  7  cm  tile,  containing  25,600  individual  elec¬ 
trokinetic  microactuators  with  250  pm  channel  diameters  arranged  on  350  pm  cen- 
ter-to-center  spacings,  as  indicated  in  Eigs.  2.4  and  2.6  for  full-scale  hydronautical  vehicle 
applications.  These  individual  electrokinetic  microactuators  were  grouped  into  1600  indi¬ 
vidual  unit  cells  arranged  in  a  40  x  40  unit-cell  pattern  on  each  tile,  with  each  unit-cell 
composed  of  a  4  x  4  array  of  microactuators.  As  indicated  in  Eig.  5.17,  every  fifth  row 
and  column  of  microactuators  in  the  tile  was  skipped  to  provide  room  for  the  resident 
unit-cell  processing  electronics.  The  top  layer,  shown  in  Eig.  5.18a,  was  fabricated  using 
conventional  MEMS  processing  steps.  The  center  layer  in  Eig.  5.186  was  fabricated  with 
the  same  porous  polymer  matrix  material  in  thin  flexible  mylar  substrate  material  suitable 
for  conformal  application  on  a  vehicle  surface.  The  following  sections  describe  these  cen¬ 
ter-layer  and  top-layer  fabrication  processes  used  to  develop  this  full-scale  MEKA-5 
hydronautical  array. 


5,7,2  MEKA-5  Center  Layer  Fabrication 

Figures  5. 19a,  6  show  the  center  layer  of  the  three-layer  MEKA-5  hydronautical  array  tile. 
Each  of  the  25,600  individual  electrokinetic  microactuators  has  a  250  pm  microchannel 
diameter  and  350  pm  center-to-center  spacing  between  microactuators  within  a  unit  cell. 
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The  channels  were  fabricated  in  flexible  mylar  substrate  material,  having  250  gm  thick¬ 
ness  to  permit  a  high  electric  field  E  across  the  individual  actuator  channels  at  a  low  volt¬ 
age  difference.  A  potential  difference  of  just  20  V  across  this  250  pm  thick  center  layer 
produces  the  necessary  field  strength  to  achieve  the  10  pL/min  equivalent  DC  flow  rate  in 
Fig.  2.5  needed  for  active  sublayer  control  on  hydronautical  vehicles,  given  the  pumping 
capability  of  the  porous  polymer  demonstrated  in  Fig.  5.8.  The  electrokinetic  porous 
polymer  matrix  structure  was  filled  in  the  liquid  state  in  the  microchannels  by  the  two 
component  polymerization  process  described  in  §5.2.  The  curing  process  produced  a 
porous  matrix  with  pore  sizes  in  the  range  of  1  pm.  Figure  5.20  shows  an  SEM  image  of 
the  typical  resulting  pore  structure.  Note  that  pore  sizes  vary  significantly,  but  are  typi¬ 
cally  1  pm  or  smaller,  indicating  a  roughly  1  MHz  theoretical  frequency  response  limit 
based  on  the  considerations  in  §4.4.  The  10  kHz  loss-less  frequency  response  demon¬ 
strated  with  essentially  the  same  porous  polymer  matrix  in  Fig.  5.9  is  thus  more  than  suffi¬ 
cient  for  the  1  kHz  frequency  response  requirements  in  Figs.  2.4  and  2.5  for  sublayer 
control  on  large  hydronautical  vehicles. 


5.7.3  MEKA-5  Top  Layer 

The  top  layer  of  the  MEKA-5  array  was  fabricated  using  photolithographic  etching  and 
other  MEMS  mass  fabrication  processes.  This  layer  provides  the  25,600  individual 
microactuator  electrodes  and  leadouts,  grouped  into  a  unit-cell  architecture  as  shown  in 
Eig.  5.21.  The  layer  was  fabricated  by  a  three-mask  MEMS  process  consisting  of  an 
“electrodes  mask”,  a  “cavities  mask”,  and  a  “nozzles  mask”,  shown  respectively  in  Eigs. 
5.21,  and  5.22a, b.  The  fabrication  process  sequence  developed  for  this  MEMS  top  layer  is 
shown  in  Table  5.1. 

The  major  steps  in  this  process  are  shown  schematically  in  Eig.  5.23.  A  525  pm  silicon 
wafer  was  first  used  as  a  temporary  substrate  for  the  top  layer.  A  2  pm  thick  layer  of  alu¬ 
minum  was  then  sputtered  on  the  top  surface  of  this  wafer  -  the  entire  MEKA-5  top  layer 
was  fabricated  on  top  of  this  aluminum  layer.  The  silicon  wafer  was  used  solely  to  pro¬ 
vide  support  for  the  top  layer  and  to  permit  handling  during  the  fabrication  process.  Elpon 
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completion  of  the  fabrication  process  described  below,  this  aluminum  layer  was  etehed 
away  as  shown  in  the  last  step  in  Fig.  5.23,  to  release  the  top  layer  from  the  silicon  wafer. 

As  regards  the  fabrication  of  the  top  layer  itself,  a  6  pm  thick  layer  of  polyimide  was  first 
spin  eoated  over  the  aluminum  layer.  Using  the  nozzles  mask  shown  in  Fig.  5.22a,  the 
polyimide  layer  was  then  photolithographieally  patterned,  unmasked,  and  exposed  to 
define  the  50  pm  nozzle  array.  A  layer  of  photoresist  was  then  used  to  coat  the  entire  sur¬ 
face.  Using  the  electrode  mask  in  Fig.  5.21,  corresponding  parts  of  this  photoresist  layer 
were  then  etehed  away.  This  defined  the  areas  where  the  electrodes  were  to  be  loeated.  A 
4000  A  layer  of  chrome  was  then  sputtered  onto  this  surface,  followed  by  another  4000  A 
layer  of  gold.  The  use  of  Cr/Au  for  the  electrodes  proved  to  be  important  for  the  fabriea- 
tion  proeess,  as  early  attempts  to  pattern  the  electrodes  using  platinum  repeatedly  failed 
beeause  this  material  would  consistently  peel  away  from  the  underlying  polyimide  layer. 
A  similar  result  oecurred  when  gold  alone  was  used.  This  Cr/Au  double-layer  was  then 
patterned  and  unmasked  using  the  electrode  mask,  and  the  photoresist  then  etched  with 
aeetone  to  expose  the  eleetrodes.  Finally,  a  125  pm  thick  layer  of  Su-8  material  was 
deposited  and  selectively  etched  using  the  cavities  mask  shown  in  Fig.  5.22b.  At  this 
point,  the  top  layer  fabrieation  was  complete,  and  the  layer  is  released  from  the  silicon 
support  wafer  by  etching  away  the  aluminum  layer  with  a  phosphorie-nitrie-aeetie  acid 
mixture. 

Fabrieation  of  the  complete  top  layer  with  this  process  sequence  was  done  via  MEMS 
Exchange  in  The  University  of  Michigan’s  College  of  Engineering  MEMS  fabrieation 
facilities.  The  three  eomponent  sublayers  that  make  up  the  top  layer  are  named  after  the 
masks  used  to  fabricate  them.  A  full-tile  view  of  the  “eleetrodes  sublayer”  is  shown  in 
Eig.  5.24a.  This  sublayer  provides  Cr/Au  electrodes  for  each  of  the  25,600  individual 
microactuators  and  the  corresponding  leadout  for  each  actuator  to  a  contact  located  near 
the  edge  of  its  unit  cell.  A  region  of  this  eleetrodes  sublayer  is  shown  in  Eigs.  5.24b, 
where  the  repeated  pattern  in  twenty  adjaeent  units  eells  can  be  seen.  A  closer  view  of  a 
single  unit-cell  in  this  strueture  is  shown  in  Eig.  5.25a,  with  the  sixteen  unit-eell  eontacts 
for  each  of  the  individual  electrodes  in  the  unit  cell,  together  with  the  eommon  voltage  bus 
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lines  that  run  between  adjacent  unit  cells.  These  voltage  bus  lines  can  be  readily  seen  in 
Fig.  5.25b.  Future  versions  of  this  MEKA-5  array  could  contain  a  digital  signal  process¬ 
ing  (DSP)  layer  that  selectively  connects  each  of  these  contacts  to  the  bus  line  with  either 
positive  or  negative  polarity,  depending  on  the  actuation  state  vector  implied  by  the  sensor 
states.  As  described  in  §3,  this  could  be  readily  done  via  a  look-up  table  approach  using  a 
simple  programmable  logic  array  (PLA),  in  place  of  more  traditional  control  logic 
approaches  based  on  comparatively  elaborate  signal  processing. 

The  “nozzles  sublayer”  is  etched  in  polyamide  and  provides  the  50  pm  nozzles  through 
which  the  pumping  induced  by  the  electrokinetic  flow  occurs.  A  full  view  of  the  “cavities 
sublayer”  is  shown  in  Fig.  5.26a.  This  sublayer,  fabricated  in  Su-8  material  as  described 
in  the  process  flow,  serves  two  specific  functions.  It  provides  the  necessary  structural  sup¬ 
port  for  the  top  layer,  and  at  the  same  time  provides  the  necessary  separation  between  the 
top  of  the  porous  polymer  matrix  structure  associated  each  microactuator  channel  in  the 
center  layer  and  the  corresponding  top-layer  electrode.  A  closer  view  of  the  cavities  sub¬ 
layer  in  Fig.  5.26b  shows  the  250  pm  diameter  cavities  on  325  pm  center-to-center  spac- 
ings,  and  reveals  the  comparatively  well-defined  edges  of  these  cavities  that  result  from 
the  fabrication  process  described  above. 


5,7,4  System  Architecture  for  MEKA-5  Microactuator  Array 

The  architecture  for  arrays  like  MEKA-5  is  based  on  a  4  x  4  unit-cell  composed  of  wall 
shear  stress  sensors  and  electrokinetic  microactuators.  The  top-surface  electrode  for  each 
microchannel  has  a  leadout  that  runs  to  a  contact  near  the  edge  of  the  unit  cell.  All  the  unit 
cells  within  the  tile  share  a  common  ground  electrode  in  the  bottom  layer.  A  power  bus 
for  the  entire  tile,  held  at  constant  reference  voltage  runs  along  horizontal  and  vertical 
lines  between  the  active  areas  of  adjacent  unit-cells.  This  redundancy  makes  the  array 
highly  fault  tolerant  to  damage.  The  look-up  table  logic  circuit  provides  a  three-state 
bridge  between  the  electrode  contact  for  each  of  the  microactuators  and  the  corresponding 
closest  power  bus  line.  As  noted  in  §3,  on  each  clock  cycle  the  actuator  state  vector  {A/} 
obtained  via  the  programmable  logic  array  from  the  sensor  state  vector  {S,}  sets  the  volt- 
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age  AiV^efOi  the  top  eleetrode  for  eaeh  aetuator.  This  produees  positive  volume  displace¬ 
ment  on  some  actuators,  negative  displacements  on  some  actuators,  and  no  action  on  the 
remaining  actuators.  The  space  between  adjacent  unit  cells  on  the  MEKA-5  array  appears 
sufficient  to  accommodate  the  relatively  simple  circuitry  needed  to  implement  this  system 
architecture. 


5.8  Laboratory  Demonstration  of  Synthetic  Vortex  Manipulation 

The  various  generations  of  electrokinetic  microactuator  arrays  shown  in  the  previous  sec¬ 
tions  have  demonstrated  key  performance  and  fabrication  requirements  for  use  of  the  elec¬ 
trokinetic  principle  in  actuator  arrays  designed  for  active  sublayer  control  under  full-scale 
vehicle  conditions.  However,  direct  flow  visualizations  showing  the  lateral  displacement 
of  naturally-occurring  sublayer  vortices  in  turbulent  boundary  layers  would  be  exceed¬ 
ingly  difficult  due  to  the  very  small  size  of  these  vortices  and  the  fast  time  scales  associ¬ 
ated  with  them  (see  §2).  Similarly,  direct  measurements  of  drag  reduction  resulting  from 
active  sublayer  control  using  a  set  of  such  electrokinetic  microactuators  would  require  a 
complete  control  system,  including  sensors  and  control  processing,  which  is  beyond  the 
scope  of  the  present  study.  For  these  reasons,  a  set  of  flow  visualization  experiments  was 
instead  conducted  to  demonstrate  lateral  displacement  of  a  synthetically-generated  stream- 
wise  vortex  in  the  vicinity  of  a  wall  by  volumetric  pumping  from  the  wall.  This  is  essen¬ 
tially  the  mechanism  by  which  the  present  electrokinetic  microactuators  work,  but  allows 
the  experiment  to  be  conducted  at  much  larger  scales. 

These  experiments  to  demonstrate  manipulation  of  synthetic  streamwise  near-wall  vorti¬ 
ces  via  wall  pumping  were  conducted  in  a  simple  wind  tunnel.  The  basic  layout  is  shown 
schematically  in  Fig.  5.27.  A  streamwise  vortex  of  circulation  F  was  synthetically  gener¬ 
ated  by  a  small  rectangular  wing  held  at  a  fixed  angle-of-attack  and  located  a  small  dis¬ 
tance  behind  the  leading  edge  of  a  flat  plate  mounted  in  the  wind  tunnel.  The  fiat  plate  is 
shown  in  Fig.  5.28,  and  was  mounted  on  a  computer-controlled  positioning  stage  that  per¬ 
mitted  precise  positioning  of  the  streamwise  vortex  within  the  field-of-view  of  the  PIV 
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measurement  system.  The  flat  plate  had  a  rounded  leading  edge  and  a  flap  to  inerease  eir- 
culation  and  thereby  prevent  boundary  layer  separation.  In  this  two-dimensional  experi¬ 
ment,  pumping  from  the  wall  was  aecomplished  through  a  slot  aligned  with  the  flow,  and 
the  vortex  displaeement  is  measured  in  the  erossplane  flow.  Because  electrokinetic  pump¬ 
ing  is  effective  only  at  the  microscale,  the  volumetric  pumping  in  this  case  was  accom¬ 
plished  by  steady  air  injection  at  various  rates  through  the  slot  from  a  pressurized  plenum. 
The  principle  aim  of  the  experiments  was  to  demonstrate  lateral  displacement  of  the 
streamwise  vortex  by  wall  pumping,  and  to  characterize  the  relative  lateral  displacement 
achieved  as  a  function  of  the  wall  pumping  rate,  the  vortex  strength,  and  the  position  of 
the  vortex  relative  to  the  location  at  which  the  wall  pumping  occurs. 


5,8.1  Wind  Tunnel  Facility 

The  experiment  was  conducted  in  a  low  speed  wind  tunnel  that  was  specifically  designed 
and  fabricated  for  this  study.  This  suction  wind  tunnel,  shown  in  Figs.  5.29-5.31,  operated 
on  an  induced  draft  created  by  the  ventilation  system  in  the  laboratory,  which  allowed  for 
simple  removal  of  the  submicron-scale  seeding  particles.  The  three  main  sections  of  this 
wind  tunnel  are  shown  schematically  in  Figs.  5.29  and  5.30.  A  seeding  section,  shown  in 
Fig.  5.29a,  is  located  at  the  extreme  upstream  end,  where  the  induced  flow  enters  the  tun¬ 
nel  and  mixes  with  seed  particles  introduced  through  four  jets  issuing  from  the  sidewalls. 
Effective  mixing  of  the  seed  particles  is  enhanced  by  turbulence  created  from  a  set  of  ran¬ 
domly  spaced  holes  of  various  diameters  through  which  the  induced  air  enters  the  tunnel. 
The  flow  then  passes  through  a  straightening  section  composed  of  a  series  of  two  screens, 
a  perforated  plate,  and  a  honeycomb  section  (also  evident  in  Fig.  5.29a)  that  reduces  the 
free  stream  turbulence  level  in  the  test  section  to  less  than  3%  at  speeds  as  low  as  0.1  m/s, 
as  measured  by  the  PIV  system.  The  test  section  has  interior  cross-sectional  dimensions 
of  36  cm  X  30  cm  and  is  96  cm  long.  The  test  section  has  four  glass  windows  for  optical 
access  of  the  PIV  cameras  through  the  side  walls,  as  can  be  seen  in  Figs.  5.30a,h,  to  mini¬ 
mize  the  optical  distortions  produced  by  the  acrylic  test  section  walls.  The  locations  of  the 
four  PIV  cameras  as  well  as  their  orientations  with  respect  to  the  wind  tunnel  are  shown  in 
Figs.  532a-d. 
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5,8,2  Stereo  Particle  Image  Velocimetry  SPIV  System 


A  unique  dual-plane  two-eolor  stereo  PIV  system  was  assembled  for  these  experiments, 
which  in  principle  can  permit  highly  resolved  measurements  of  all  nine  components  in  the 
velocity  gradient  tensor  field  Vu(x,t).  The  basic  layout  of  the  DSPIV  system  is  shown  in 
Fig.  5.33.  The  system  uses  four  frequency-doubled  Nd:YAG  lasers  (Spectra-Physics)  to 
provide  beams  at  532  nm  (green)  with  roughly  10  nsec  pulse  length.  Two  of  these  beams 
are  used  to  pump  a  dye  laser,  which  provides  output  centered  around  623  nm  (red).  The 
green  and  red  beams  from  two  of  the  lasers  are  formed  into  parallel  light  sheets  that  inter¬ 
sect  the  synthetic  streamwise  vortical  structure.  Each  of  the  green  and  red  light  sheets  is 
double-pulsed  by  time-separating  the  corresponding  pulses  from  the  two  pairs  of  NdiYAG 
lasers.  Two  stereo  PIV  systems  (LaVision)  are  used  to  record  particle  images  from  the 
green  and  red  light  sheets.  Optical  filters  allow  each  pair  of  cameras  to  see  only  one 
light-sheet  color.  The  two  SPIV  systems  are  slaved  together  and  operated  via  a  custom 
master-slave  software  arrangement  developed  for  this  application  by  LaVision.  Process¬ 
ing  of  the  two  separate  stereo  PIV  image  pairs  produces  the  velocity  field  u(x,t)  and,  more 
importantly,  the  full  velocity  gradient  field  Vu(x,t),  as  indicated  in  Fig.  5.34. 

In  the  current  study,  this  system  was  operated  in  the  single  stereo  PIV  mode.  The  dou- 
ble-frame/double-exposure  measurements  were  made  with  two  fast  shuttered  PIV  cameras 
(PCO  SensiCam  SVGA,  1280  x  1024  pixels,  6.7  mm  x  6.7  mm  pixel  size).  The  double 
image  from  both  cameras  was  stored  in  a  single  buffer.  There  are  thus  four  frames  in  each 
buffer,  with  the  order  of  each  frame  shown  in  Fig.  5.35.  The  3-component  vector  calcula¬ 
tion  was  performed  using  the  3D  stereo  cross-correlation  mode  in  the  Davis  6.2.0  (LaVi¬ 
sion)  software  package.  As  is  common  in  the  evaluation  of  PIV  recordings,  the  images 
were  divided  into  32  x  32-pixel  “interrogation  windows”.  The  cross-correlation  algorithm 
operates  on  the  intensities  inside  each  interrogation  window,  yielding  one  velocity  vector 
for  each  window.  With  the  interrogation  window  size  set  to  32  x  32  pixels,  using  a  camera 
with  a  CCD  chip  having  1280  x  1024  pixels,  the  image  is  divided  into  40  x  32  cells.  The 
processing  algorithm  used  a  50-percent  overlap  between  neighboring  interrogation  win¬ 
dows,  which  produced  a  final  image  comprised  of  an  array  of  79  x  63  vectors  {i.e.,  4977 
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uniformly-spaced  vectors).  Furthermore,  an  adaptive  multipass  algorithm  was  used  in  the 
Davis  software  which  allowed  for  successively  smaller  final  interrogation  window  sizes 
that  would  be  possible  without  adaptive  window  shifting.  This  approach  also  improves 
the  spatial  resolution  of  the  final  vector  field. 

Two-dimensional  velocity  vectors  were  calculated  for  each  of  Camera  1  and  Camera  2  by 
applying  the  stereo  cross-correlation  to  four-frame  buffers  of  the  type  shown  in  Fig.  5.35. 
The  three-dimensional  vector  at  each  point  was  then  calculated  from  the  two-dimensional 
vectors  obtained  from  Camera  1  and  Camera  2  by  using  the  local  viewing  angle  and  rela¬ 
tive  orientation  of  the  two  cameras.  This  was  achieved  by  a  transformation  defined  from 
mapping  functions  created  for  Camera  1  and  Camera  2  during  the  stereo  three-dimen¬ 
sional  calibration  process  {e.g.,  Soloff  et  al  1997).  The  resulting  system  of  equations  is 
overspecifed,  and  thus  allows  for  an  additional  check  to  determine  if  the  computed  vectors 
are  valid  to  improve  the  accuracy  of  the  measurements. 

The  two  532  nm  green  beams,  with  energies  in  excess  of  100  mJ,  were  combined  in  a  laser 
sheet  that  crossed  the  test  section  perpendicular  to  the  freestream  flow.  This  laser  sheet 
was  approximately  1  mm  thick  and  2  cm  wide.  The  laser  fiashlamps  and  Q-switches  oper¬ 
ate  at  10  Hz  to  produce  10  ns  pulse  widths.  The  time  separation  between  the  two  laser 
pulses  that  provided  the  highest  good-vector  percentage  was  found  to  be  90  qs  for  this 
experiment.  This  time  delay,  as  well  as  the  laser  Q-switches  and  fiashlamps,  were  con¬ 
trolled  by  a  PC-based  programmable  timing  unit  (PTU)  integrated  with  the  rest  of  the  PIV 
system  software  (LaVision  DaVis  6.2).  The  PTU  also  synchronizes  the  triggering  of  the 
lasers  and  cameras. 


5,8,3  PIV  Results 

All  the  PIV  results  shown  in  this  section  cover  a  field  of  view  of  12  mm  x  15  mm.  The 
measurements  thus  span  from  the  surface  of  the  plate  to  12  mm  above  the  plate  (the 
y-direction)  and  approximately  15  mm  in  the  transverse  direction  (x  direction).  For  con¬ 
venience,  the  center  of  the  wall  actuator  is  defined  as  (x,  y)  =  (0,  0).  The  measurements 
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were  made  20  cm  in  the  downstream  direction  (z  direction)  from  the  leading  edge  of  the 
plate. 

The  first  set  of  images,  shown  in  Fig.  5.36,  gives  the  transverse  and  streamwise  velocity 
field  for  the  laminar  boundary  layer  over  the  flat  plate  at  a  streamwise  velocity  of  17=  0.6 
m/s.  In  this  baseline  case,  the  streamwise  vortex  generator  was  removed,  so  that  there  was 
no  vortex  in  the  flow.  The  Blasius  solution  for  the  laminar  boundary  layer  thickness  on  a 
fiat  plate  indicates  an  expected  value  for  of  9.5  mm  at  these  conditions.  This  value  is 
comparable  to  the  results  shown  in  Fig.  5.37,  where  is  found  to  be  approximately  8 
mm.  The  use  of  a  flap  at  the  trailing  edge  of  the  plate,  set  to  a  small  angle  to  provide  pos¬ 
itive  circulation  around  the  plate  and  thereby  produce  a  slight  favorable  pressure  gradient, 
results  in  a  reduced  growth  rate  of  the  boundary  layer  and  accounts  for  the  smaller  mea¬ 
sured  value  of  Sggo/,.  Figure  5.37  also  verifies  that  the  boundary  layer  velocity  profile 
remains  essentially  the  same  at  the  various  transverse  locations  shown.  Typical  vorticity 
and  shear  strain  rate  fields  for  this  laminar  boundary  layer  are  shown  in  Figs.  5.38  and 
5.39,  respectively. 

Next,  the  velocity  field  for  the  boundary  layer  with  the  synthetic  streamwise  vortex,  cre¬ 
ated  by  the  small  wing  upstream  of  the  measurement  location,  is  shown  in  Fig.  5.40.  In 
this  second  baseline  case,  there  was  no  flow  issuing  from  the  wall  actuator.  The  center  of 
the  vortex  is  seen  to  be  located  approximately  4  mm  above  the  plate,  and  less  than  1  mm  to 
the  left  of  the  center  of  the  slot  actuator.  The  downward  motion  of  fluid  created  by  the 
vortex  on  its  left  side  significantly  decreases  the  thickness  of  the  boundary  layer  in  that 
region,  as  seen  in  Fig.  5.41  at  x  =  -5  mm.  The  opposite  is  true  on  the  right  side  of  the  vor¬ 
tex,  as  seen  in  Fig.  5.41  for  x  =  5  mm.  At  x  =  0  in  Fig.  5.41,  which  is  close  to  the  center  of 
the  vortex,  the  streamwise  component  of  the  velocity  shows  an  inflection  point,  as  would 
also  be  expected.  In  addition,  due  to  the  presence  of  the  vortex  an  accumulation  of  low 
speed  fluid  occurs  near  the  wall  at  the  right  side  of  the  vortex.  In  Fig.  5.42,  the  instanta¬ 
neous  and  ensemble-averaged  vorticity  fields  cOz  associated  with  the  velocity  fields  in  Fig. 
5.40  show  the  presence  of  the  streamwise  vortex  with  a  reasonably  well-defined  core. 
Similarly,  the  presence  of  the  vortex  in  this  viscous  region  generates  the  shear  strain  rate 
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field  shown  in  Fig.  5.43. 


The  PIV  results  presented  up  to  this  point  corresponds  to  the  wall  actuator  being  “turned 
off’,  so  that  there  is  no  blowing  (or  suction)  by  the  actuator  at  the  wall.  The  resulting  con¬ 
figuration  is  shown  schematically  in  Fig.  5.44a,  where  a  vortex  of  strength  F  is  located  at 
a  distance  h  =  A  mm  above  the  wall  and  at  a  distance  b  =  \  mm  from  the  center  of  the  slot 
actuator.  When  the  wall  actuator  is  “turned  on”,  so  that  there  is  a  constant  volume  flow 
rate  q  per  unit  actuator  length  introduced  at  the  wall,  then  the  resulting  configuration  is  as 
shown  schematically  in  Fig.  5.44h.  In  this  case,  when  the  vortex  has  the  sense  of  circula¬ 
tion  as  shown  and  is  located  to  the  left  of  the  actuator,  then  the  volumetric  pumping  at  the 
wall  should  displace  the  vortex  to  the  left  by  a  distance  /  from  its  original  location. 

Corresponding  PIV  results  for  the  velocity  field  under  the  influence  of  such  wall  actuation 
are  shown  in  Figs.  5.45-5.48.  As  anticipated,  the  vortex  is  seen  to  be  displaced  in  the  lat¬ 
eral  direction  to  the  left  as  a  result  of  the  low-speed  fluid  being  pumped  by  the  actuator. 
Furthermore,  this  low-speed  flow  is  accumulated  to  the  right  of  the  actuator,  as  can  be 
seen  in  Fig.  5.45,  due  to  the  influence  of  the  region  of  the  vortex  with  an  upward  motion. 
This  tends  to  modify  the  boundary  layer  as  shown  in  Fig.  5.46.  Also,  in  the  vorticity  in 
Fig.  5.47,  the  vortex  core  can  be  readily  seen  to  be  moved  to  the  left  in  relation  to  Fig. 
5.41,  directly  demonstrating  the  lateral  displacement  of  the  vortex.  Similarly,  a  reorienta¬ 
tion  of  the  shear  strain  rate  occurs  under  the  presence  of  wall  actuation,  as  seen  in  Fig. 
5.48.  Similar  PIV  results  were  obtained  at  several  different  values  of  the  volume  flow  rate 
q  introduced  by  the  wall  actuator. 

For  a  given  initial  position  of  the  vortex  relative  to  the  actuator,  the  lateral  displacement  of 
the  vortex  induced  by  the  actuation  was  seen  to  depend  on  the  vortex  strength  F  and  the 
volume  flow  rate  q.  The  PIV  results  indicate  a  strong  correlation  between  the  relative  lat¬ 
eral  displacement  l/h  of  the  vortex  and  the  “relative  pumping  strength”  qlT,  as  shown  in 
Fig.  5.49.  Furthermore,  this  plot  shows  that  there  is  a  relative  pumping  strength  limit  at 
which  a  plateau  is  reached,  beyond  which  further  increase  in  qlT  produces  little  or  no 
additional  relative  lateral  displacement  of  the  vortex.  One  possible  explanation  for  this  is 
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that,  as  was  noted  in  Fig.  2.3,  as  the  pumping  strength  q  inereases  so  does  the  loeal  exit 
Reynolds  number  of  the  aetuator  flow.  Eventually,  the  flow  produced  by  the  actuator  will 
reach  sufficiently  high  Reynolds  number  that  it  will  exhibit  a  jet-like  behavior,  rather  than 
the  low  Reynolds  number  radially-dominated  flow  indicated  in  Fig.  5.44.  This  would  act 
to  reduced  the  effectiveness  of  the  actuator  in  producing  lateral  displacement  of  the 
streamwise  vortex,  consistent  with  the  results  seen  in  these  experiments. 
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to  AC/DC  power  supply  to  oscilloscope 


syringe  pump 


Figure  5.1.  Experimental  configuration  used  for  testing  of  elementary  packed-chan¬ 
nel  actuators  described  in  §5.1.  Displacement  of  the  syringe  pump  permitted  volume 
flow  rate  measurements  for  various  applied  field  strengths  and  electrokinetic  pumping 
parameters.  The  pressure  transducer  permitted  measurements  of  the  pressure  achieved 
in  plugged  actuators. 
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Figure  5.2.  Basic  three-layer  design  of  MEKA-0  through  MEKA-5  electrokinetic 
microactuators  arrays  in  §5,  with  center  layer  containing  porous  polymer  matrix  struc¬ 
ture  in  which  electrokinetic  flow  occurs.  The  1  pm  pore  size  indicates  a  theoretical  fre¬ 
quency  response  limit  near  500  kHz.  Channel  diameters  ranged  from  mm-scales  for  the 
MEKA-0  and  MEKA-2  UAV-scale  microactuator  arrays,  to  250  pm  for  the  MEKA-5 
hydronautical-scale  array,  which  contained  25,600  such  individual  microactuator  chan¬ 
nels. 
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Figure  5.3.  Typical  center-layer  structure  showing  electrokinetic  driver  channels  filled 
with  porous  polymer  matrix  structure.  Channels  are  fabricated  to  match  the  unit-cell 
actuator  pattern,  and  completely  fdled  with  porous  polymer.  Electric  field  applied 
along  channel  length  creates  electrokinetic  flow  in  porous  matrix. 
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(a) 


(b) 


SHE5H5 


Figure  5.4  .  Typical  top-layer  structure  showing  (a)  actuator  nozzles  at  top  surface 
used  to  produce  point  volume  source  flow,  (b)  cross-sectional  view  showing  cavity, 
electrode,  and  channels  for  each  actuator,  and  (c)  electrodes  and  leadouts  to  unit-cell 
contacts. 
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Figure  5.5.  MEKA-0  3x3  electrokinetic  microactuator  array,  showing  (a) 
experimental  test  set-up,  (b)  array  in  electrolyte  reservoir,  (c)  array  fabricated  in 
glass,  and  (d,e)  front  and  back  surfaces  of  array  showing  metallization  and  porous 
polymer  in  channels. 


(a) 


(b) 


Figure  5.6.  Closeup  views  of  MEKA-0  eleetrokinetic  mieroaetuator  array  prior 
to  polymer  filling,  showing  (a,b)  top-surface  ring  electrodes  for  each  actuator,  (c) 
bottom-surface  common  electrode,  and  (d)  channel  edge  chipping  resulting  from 
mechanical  drilling  in  glass  wafer. 
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(c) 


(d) 


Figure  5.6.  (concluded) 


93 


to  AC/DC  power  supply 


Figure  5.7.  Schematic  of  experimental  apparatus  used  for  steady  flow  rate  measure¬ 
ments  with  MEKA-1  electrokinetic  microactuator  array.  Applied  electric  field  strength 
was  determined  from  voltage  of  top-surface  electrode  and  thickness  of  center  layer. 
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Time  (min) 


Figure  5.8.  Experimental  results  for  steady  electrokinetic  pumping  performance 
achieved  with  MEKA-1  array  for  two  different  applied  voltages.  Total  displaced  vol¬ 
ume  increased  linearly  with  time,  and  varied  with  field  strength  as  indicated  in  (4.23). 
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Figure  5.9.  Measured  frequency  response  of  an  electrokinetic  actuator  for  a  sinusoi 
dally  time-varying  applied  field  of  various  frequencies  co,  showing  results  from  a  low 
frequency  pressure  transducer  and  a  high-frequency  microphone.  Solid  line  gives  co“ 
roll-off  corresponding  to  loss-less  AC  performance. 
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(a) 


Figure  5.10.  Electrode  and  leadout  pattern  for  10  x  10  MEKA-1  array  ,  showing 
(a)  individual  top-surface  ring  electrode  and  leadout  for  each  individual  actuator 
channel,  and  (b)  common  bottom-surface  electrode. 
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Figure  5.11.  Comparison  of  actuator  channel  edge  quality  achieved  by  (a) 
mechanical  drilling  of  glass  substrate  on  MEKA-0  array,  and  (b)  CO2  laser  drilling 
of  glass  substrate  on  MEKA-1  array. 
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(b) 


Figure  5.12.  Photographs  of  1  pm  porous  polymer  matrix  structure  filled  in  10  x 
10  MEKA-1  array. 
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(a) 


(b) 


Figure  5.13.  Poor  channel  quality  obtained  by  CO2  laser  drilling  of  300  gm  diam¬ 
eter  holes  on  500  gm  center-to-center  spacing  in  2  mm  thick  glass  substrate  on 
MEKA-2  array,  showing  (a)  entry  side,  and  (b)  exit  side. 
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(a) 


(b) 
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Figure  5.14.  Actuator  channels  fabricated  by  CO2  laser  drilling  of  300  gm  diame¬ 
ter  holes  on  1000  gm  center-to-center  spacing  in  2  mm  thick  acrylic  substrate  on 
MEKA-3  array,  showing  (a)  complete  wafer-scale  view,  and  (b)  exit-side  view  of 
actuator  channels. 


101 


(a) 


Figure  5.15.  Close-up  views  of  MEKA-3  actuator  channels  on  (a)  entry  side,  and 
(b)  exit  side. 
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Figure  5.16.  Layout  of  electrokinetic  microactuator  channels  in  center  layer  of 
MEKA-5  full-scale  hydronautical  array,  showing  25,600  individual  microactuator 
channels  grouped  into  a  40  x  40  pattern  of  unit  cells,  each  containing  a  4  x  4  unit-cell 
structure  on  325  pm  center-to-center  spacing.  Array  is  shown  at  actual  size. 
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Figure  5. 
how  full 
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(a)  Top  Layer 


(b) 


Middle  Layer 


Polymer 


Figure  5.18.  Basic  elements  fabricated  in  25,600-element  MEKA-5  hydronautical- 
scale  array,  showing  (a)  top  layer  with  cavities,  electrodes  and  leadouts,  and  nozzles, 
and  (b)  center  layer  with  actuator  channel  and  porous  polymer  matrix.  Individual  sub¬ 
layers  are  shown  in  Figs.  5.19-5.22.  MEMS  fabrication  process  steps  are  shown  in  Eig. 
5.23. 
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(a) 


Figure  5.19.  (a)  Center  layer  of  25,600-element  MEKA-5  hydronautieal-scale 

array  fabricated  in  thin  flexible  mylar,  showing  tile-scale  view  250  gm  diameter 
actuator  channels  filled  with  porous  polymer  matrix  structure,  (b)  Close-up  view  of 
250  gm  diameter  actuator  channels  filled  with  porous  polymer  matrix  structure  in 
MEKA-5  hydronautieal-scale  array,  showing  individual  4x4  unit-cells. 
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Figure  5.20.  SEM  micron-scale  photograph  of  the  porous  polymer  matrix  structure 
resulting  from  two-component  polymer  fabrication  process  of  the  type  used  in  the 
MEKA-5  array. 
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Figure  5.21.  Layout  of  top-surface  electrodes  sublayer  mask  for  25,600-element 
MEKA-5  array,  showing  disk  electrodes  with  leadouts  to  contacts  near  edge  of  each 
unit  cell,  and  common  voltage  bus  lines  running  horizontally  and  vertically  between 
adjacent  unit  cells. 
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(a) 


"Nozzles"  Mask  •  •  •  • 


•  •  •  • 


•  •  •  • 


200  p.m 


•  •  •  • 


(b) 


"Cavities"  Mask 


200  p.m 


Figure  5.22.  Additional  masks  used  in  MEMS-fabrieation  of  top  layer  for 
MEKA-5  hydronautical-scale  array,  showing  (a)  nozzles  mask,  and  (b)  cavities 
mask. 
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Process  Step 

Process  Description 

1 

Prepare  starting  wafers 

•  0  p-type,  single  side  polished,  <100>,  1-20  [W-cm] 

•  Rinse:  DI-H2O,  2  min;  spin  dry 

2 

Sputter:  A1  release  layer 

•  Target  #3,  Sputter  20000  A  A1  at  5  A 

•  7  mT  Ar  with  a  pre-sputtet  time  of  2  min. 

3 

Mask  1 :  Polyimide  Nozzles  Mask 

•  Spin:  PI  2731  @  500  rpm/8  sec,  4000  rpm/30  sec 

•  Softbake:  75  min.  @  55° 

•  Expose:  on  EV  Aligner  for  75  sec 

•  Develop:  to  be  done  on  spinner 

•  Puddle  DE9040  for  ~70  secs 

•  Transition  to  RI9180  with  250  rpm/10  sec 

•  Rinse  for  ~30  sec  at  3000  rpm/30  sec 

•  Postbake:  30  min.  @  190°C 

•  Postbake:  30  min.  @  300°C 

4 

Mask  2:  Metallization  Mask 

•  Drybake:  15  min.  @110°C 

•  Spin:  1827  @4Krpm 

•  Softbake:  30  min.  @  90° 

•  Expose:  26.5  sec  on  EV 

•  Develop:  1.1  minute  in  MF3 19 

5 

Etch  Oxide  (~10  seconds  in  Buffered  HE) 

6 

Sputter  Cr/Au 

•  Load  into  chamber  immediately  after  Step  5 

•  Pump  down  to  5x10"^  Torr(1.5  hr) 

•  Target  #  1 ,  Sputter  400  A  Cr  at  650  W 

•  7  mT  Arwith  a  pre- sputter  time  of  2  min,  3T5" 

•  Target  #3,  Sputter  4000  A  Au  at  0.5  A 

•  7  mT  Ar  with  a  pre-sputter  time  of  2  min.,  14' 

7 

Lift-Off  of  Cr/Au 

•  Lift-off:  Soak  in  Acetone  overnight 

•  Sonicate  in  ultrasonic  bath  for  ~10  sec 

•  Rinse  with  IPA  and  DI  for  3  minutes;  spin  dry  and  inspect  for  flakes 

8 

Mask  3:  SU-8 

•  Drybake:  10  min.  @  110°C 

•  Spin:  HMDS/SU-8  25  @5  sec/15  sec  500  rpm/690  rpm 

•  Softbake:  5  min.  @  70°,  50  min.  @  90° 

•  Expose:  on  EV  Aligner 

•  Bake:  1  min.  @  50°,  13  min.  @  90° 

•  Develop:  SU-8  developer  in  beaker  for  5  min., 

•  Change  solution,  7  min.,  30  sec  spray 

•  Rinse:  DI-H2O,  2  min.;  spin  dry 

9 

Device  Release 

Table  5.1.  MEMS  process  sequence  for  fabrication  of  top  layer  in  MEKA-5  array. 
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(a)  Aluminum  is  deposited  on  a  525  pm  silicon  substrate  (DC  magnetron  sputter) 


Al  (2  pm) 

Substrate 


(b)  Polyimide  is  deposited  (spin  cast) 


Polyimide  (6  pm) 
Al 

Substrate 


(c)  Patterned  and  unmasked  polyimide  is  exposed  with  AZ  400K  to  define  tthe  50  pm  nozzles 


Polyimide 

Al 

Substrate 


(d)  Photoresist  (Shipley  1827)  coat  is  deposited  (spin  cast) 


Photoresist 

Polyimide 

Al 

Substrate 


(e)  Patterned  and  unmasked  photoresist  is  isotropically  wet  etched  with  MF-319 


Photoresist 

Polyimide 

Al 

Substrate 


Figure  5.23.  MEMS  process  sequence  (a-j)  developed  for  top  layer  of  MEKA-5 
array. 


Ill 


(f)  Cr  and  Au  are  deposited  (DC  magnetron  sputter) 


Cr/Au 

Polyimide 

Al 


Substrate 


(h) 


Su-8  is  deposited  (spin  cast) 


Su-8 

Cr/Au 

Poiyimide 

Ai 


Substrate 


(i)  Patterned  and  unmasked  Su-8  is  etched  with  Su-8  deveioper 

Su-8 

-  -  -  -  -  -  -  -  Cr/Au 

II  II  II  ^  I  Poiyimide 

-  Ai 

Substrate 


0) 


Aiuminum  etch  with  phosporic-nitric-acetic  acid  mixture 


Su-8 

Cr/Au 

Poiyimide 


Figure  5.23.  (concluded) 
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(a) 


(b) 


Figure  5.24.  (a)  Tile-scale  view  of  electrodes  sublayer  fabricated  via  the  process 

sequence  in  Fig.  5.23  for  the  25,600-element  MEKA-5  hydronautical-scale  array. 
(b)  Close-up  view  of  the  resulting  unit-cell  structure  in  the  electrodes  sublayer  for 
the  25,600-element  MEKA-5  hydronautical-scale  array. 
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(a) 


Figure  5.25.  Unit-cell-scale  view  of  the  electrodes  sublayer  for  the  25,600-ele¬ 
ment  MEKA-5  hydronautical-scale  array,  showing  (a)  a  single  unit  cell,  and  (b)  the 
bus  structure  separating  adjacent  unit  cells. 
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(a) 


Figure  5.26.  (a)  Tile-scale  view  of  cavities  sublayer  fabricated  in  Su-8  via  the 

process  sequence  in  Fig.  5.23  for  the  25,600-element  MEKA-5  hydronautical-scale 
array,  (b)  Unit-cell-scale  view  of  the  cavities  sublayer  showing  a  single  unit  cell. 


115 


synthetic  streamwise  vortex  honeycomb 


test  section 


Figure  5.27.  Schematic  showing  basic  layout  of  the  large-scale  wind  tunnel  experiments  used  to  demonstrate  synthetic  streamwise 
vortex  manipulation  by  volumetric  pumping  from  a  wall  actuator,  indicating  wind  tunnel,  flat  plate,  and  vortex  generation  by  thin  air¬ 
foil.  Stereo  PIV  measurements  were  used  to  measure  lateral  displacement  of  the  streamwise  vortex  by  various  levels  of  wall  pump¬ 
ing. 


Figure  5.28.  Flat  plate  used  for  the  present  experiments  on  streamwise  vortex  manip¬ 
ulation,  showing  rounded  leading  edge,  trailing  edge  eirculation  flap,  and  digital 
micropositioning  stage. 
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(a) 


(b) 


Figure  5.29.  Induced  draft  wind  tunnel  assembled  for  experiments  on  streamwise 
vortex  manipulation,  showing  (a)  the  seeding  section  and  the  flow  straightening 
section,  and  (b)  the  test  section  with  slots  to  accomodate  the  positioning  stage. 
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Figure  5.30 


Wind  tunnel  used  for  the  experiments  on  streamwise  vortex  manipulation,  showing  stereo  PIV  aceess  windows. 


(a) 


(b) 


Figure  5.31.  Laboratory  configuration  for  the  stereo  PIV  measurements,  showing 
(a)  the  master/slave  computers  for  dual-plane  stereo  PIV  capabililty,  and  (b)  the 
four  Nd;YAG  lasers  and  the  test  section  installed  on  the  optical  table. 
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Figure  5.32.  Major  components  of  the  dual-plane  stereo  PIV  system,  showing  (a) 
the  dye  laser,  (b)  the  four  PIV  cameras,  and  (c,d)  the  wind  tunnel  test  section  with 
sheet-forming  optics  and  dual-plane  stereo  PIV  cameras. 
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(d) 


Figure  5.32.  (concluded) 
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4  Nd:YAG  Lasers 


Figure  5.33.  Basic  components  and  configuration  for  the  two-color  dual-plane  stereo 
PIV  system. 
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Figure  5.34.  Two-color  dual-plane  stereo  PIV  concept  for  measurements  of  all  nine 
components  of  the  velocity  gradient  tensor  field  Vu(x,t),  consisting  of  all  three  compo¬ 
nents  of  the  vorticity  vector  field  and  all  six  components  of  the  strain  rate  tensor  field, 
via  differentiation  of  the  three-component  velocity  vector  fields  in  two  closely-spaced 
PIV  measurement  planes. 
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Laser  1  pulse 


Camera  1,  first  frame 


Laser  2  pulse 


Laser  1  pulse 


Laser  2  pulse 


Camera  1 ,  second  frame 


Camera  2,  first  frame 


Camera  2,  second  frame 


Figure  5.35.  Typical  PIV  images  taken  with  stereo  PIV  sytem  using  double  frame 
double  exposure  mode. 
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Figure  5.36.  Typical  results  from  PIV  measurements  of  (a)  the  instantaneous 
transverse  velocity  vector  field  and  (b)  the  ensemble  mean  transverse  velocity  field, 
with  color  contours  giving  the  streamwise  velocity  in  the  boundary  layer. 
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Figure  5.37.  Measured  mean  streamwise  velocity  profiles  along  the  wall-normal  direction  at  several  transverse  locations. 
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(b) 
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Figure  5.38.  Typical  results  from  PIV  measurements  of  (a)  the  instantaneous 
transverse  velocity  vector  field  and  (b)  the  ensemble  mean  transverse  velocity  field, 
with  color  contours  giving  the  streamwise  vorticity  field  in  the  boundary  layer. 
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Figure  5.39.  Typical  results  from  PIV  measurements  of  (a)  the  instantaneous 
transverse  velocity  vector  field  and  (b)  the  ensemble  mean  transverse  velocity  field, 
with  color  contours  giving  the  shear  strain  field  in  the  boundary  layer. 
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Figure  5.40.  Typical  results  from  PIV  measurements  of  (a)  the  instantaneous 
transverse  velocity  vector  field  and  (b)  the  ensemble  mean  transverse  velocity  field, 
with  color  contours  giving  the  streamwise  velocity  in  the  presence  of  a  synthetic 
streamwise  vortex  in  the  boundary  layer. 
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Figure  5.41.  Measured  mean  streamwise  velocity  profiles  along  the  wall-normal  direction  at  several  transverse  locations  for  the 
same  flow  conditions  as  in  Fig.  5.43. 
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Figure  5.42.  Typical  results  from  PIV  measurements  of  (a)  the  instantaneous 
transverse  velocity  vector  field  and  (b)  the  ensemble  mean  transverse  velocity  field, 
with  color  contours  giving  the  streamwise  vorticity  field  in  the  presence  of  a  syn¬ 
thetic  streamwise  vortex  in  the  boundary  layer. 
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Figure  5.43.  Typical  results  from  PIV  measurements  of  (a)  the  instantaneous 
transverse  velocity  vector  field  and  (b)  the  ensemble  mean  transverse  velocity  field, 
with  color  contours  giving  the  shear  strain  field  in  the  presence  of  a  synthetic 
streamwise  vortex  in  the  boundary  layer. 
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(a) 


Figure  5.44.  Schematic  indicating  the  synthetic  streamwise  vortex  location  with 
respect  the  center  of  the  wall  actuator  under  no  wall  actuation  (a)  and  wall  actuation 
(b). 
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Figure  5.45.  Typical  results  from  PIV  measurements  of  (a)  the  instantaneous 
transverse  velocity  vector  field  and  (b)  the  ensemble  mean  transverse  velocity  field, 
with  color  contours  giving  the  streamwise  velocity  in  the  presence  of  a  synthetic 
streamwise  vortex  being  displaced  by  wall  actuation. 
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Figure  5.46.  Measured  mean  streamwise  velocity  profiles  along  the  wall-normal  direction  at  several  transverse  locations  for  the 
same  flow  conditions  as  in  Fig.  5.48. 
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Figure  5.47.  Typieal  results  from  PIV  measurements  of  (a)  the  instantaneous 
transverse  velocity  vector  field  and  (b)  the  ensemble  mean  transverse  veloeity  field, 
with  color  contours  giving  the  streamwise  vortieity  field  in  the  presence  of  a  syn¬ 
thetic  streamwise  vortex  being  displaeed  by  wall  actuation. 
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Figure  5.48.  Typical  results  from  PIV  measurements  of  (a)  the  instantaneous 
transverse  velocity  vector  field  and  (b)  the  ensemble  mean  transverse  velocity  field, 
with  color  contours  giving  the  shear  strain  field  in  the  presence  of  a  synthetic 
streamwise  vortex  being  displaced  by  wall  actuation. 
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Figure  5.49.  Measured  relative  lateral  displaeement  Hh  of  the  streamwise  vortex 
as  a  funetion  of  the  relative  pumping  strength  q/T,  following  the  notation  in  Fig. 
5  Alb.  For  suffieiently  low  pumping  rates,  these  results  from  the  PIV  measure¬ 
ments  demonstrate  effeetive  lateral  displaeements  of  streamwise  vortiees  by  wall 
actuation. 
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CHAPTER  6 


WALL  SHEAR  STRESS  SENSOR  ARRAY 


The  system  developed  in  this  study  for  aetive  sublayer  eontrol  uses  the  signals  from  a  set 
of  mieroseale  wall  shear  stress  sensors,  eolloeated  at  the  unit-eell  level  with  the  microactu¬ 
ators  described  in  §5,  as  inputs  to  the  simple  look-up  table  control  logic  described  in  §3. 
This  chapter  first  reviews  the  current  state-of-the-art  in  microscale  wall  shear  stress  sen¬ 
sors,  and  then  summarizes  the  fundamental  heat  transfer  considerations  that  are  relevant  to 
the  design  of  such  microscale  sensors.  It  then  outlines  a  simple  sensor  device  based  both 
on  these  fundamental  considerations  and  on  the  results  of  previous  studies  reported  in  the 
literature.  It  also  develops  a  MEMS  process  based  on  the  current  state-of-the-art  from  pre¬ 
vious  studies  that  in  principle  would  permit  fabrication  of  such  collocated  sensors  in  the 
top  layer  of  the  basic  three-layer  design  and  integration  with  signal  processing  electronics 
of  the  type  noted  in  §3.  The  resulting  predicted  performance  of  this  sensor  array  meets  the 
requirements  in  Figs.  2.4-2. 6  for  active  sublayer  control  on  full-scale  vehicles  under  real¬ 
istic  operating  conditions. 


6.1  Review  of  Micro  Wall  Shear  Stress  Sensors 

Current  state-of-the-art  approaches  for  wall  shear  stress  measurements  with  microscale 
devices  are  essentially  all  based  on  either  “thermal  sensors”  or  “floating  element  sensors”. 
The  former  typically  yields  higher  accuracies,  however  an  equally  important  consideration 
for  measurements  under  real  vehicle  operating  conditions  is  that  such  thermal  sensors  can 
be  built  with  significantly  simpler  structures  that  require  no  moving  parts,  and  involve 
only  cavities  that  are  fully  sealed  and  thus  do  not  easily  become  contaminated  or  plugged. 
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Such  thermal  sensors  ean  also  be  readily  conformed  to  curved  surfaees.  However,  thermal 
sensors  rely  on  an  indireet  measurement  that  relates  the  shear  stress  to  the  heat  transfer 
rate  from  the  sensor  to  the  fluid.  As  a  consequenee,  one  of  the  central  design  issues 
involves  reducing  the  heat  losses  to  the  substrate,  since  this  ean  seriously  reduce  the  sensi¬ 
tivity  of  the  sensor.  Floating  element  sensors,  on  the  other  hand,  allow  for  a  direct  mea¬ 
surement  of  the  wall  shear  stress  that  does  not  rely  on  indireet  correlations  between  the 
shear  stress  and  some  other  physical  process,  such  as  heat  transfer.  Their  main  disadvan¬ 
tage  is  that  the  principle  on  whieh  they  operate  fundamentally  requires  moving  parts.  As  a 
result,  sueh  sensors  are  inherently  more  complex  to  fabricate  than  thermal  sensors,  and 
from  an  operational  point  of  view  they  are  more  fragile  and  contain  gaps  and  small  crev¬ 
ices  that  can  easily  beeome  contaminated  or  plugged  in  a  realistic  operating  environment. 
Such  sensors  also  cannot  be  as  readily  conformed  to  curved  surfaces. 

Several  approaches  for  developing  practical  wall  shear  stress  sensors  based  on  these  two 
approaches  can  be  found  in  the  literature  {e.g.,  Huang  et  al  1996,  1999;  Kimura  et  al  1999; 
Liu  etal  1999;  Gad-el-Hak  1989,  1994,  1996;  Lofdahl  &  Gad-el-Hak  1999;  Padmanabhan 
et  al  1995;  Shajii  et  al  1992;  Schmidt  et  al  1988;  Breuer  2000).  Of  particular  relevance  for 
the  present  study,  Schmidt  et  al  (1988)  fabricated  a  mieroscale  floating  element  sensor 
with  a  sensing  area  of  500  x  500  qm^  using  MEMS  teehniques,  and  Huang  et  al  (1996) 
developed  a  microscale  thermal  sensor  to  measure  surfaee  shear  stresses,  but  noted  that 
the  resulting  sensitivity  of  their  sensors  was  less  than  desirable. 

For  wall  shear  stress  measurements  of  the  streamwise  structures  in  the  near-wall  region  of 
turbulent  boundary  layers,  comparatively  large  arrays  of  such  sensors  would  be  needed  to 
identify  these  structures.  Arrays  of  thermal  sensors  would  appear  to  be  more  appropriate 
for  this  purpose  due  to  their  inherently  higher  robustness,  as  noted  above,  and  the  fact  that 
their  inherent  simplicity  more  readily  permits  fabrication  in  dense  arrays  collocated 
among  actuators  and  processing  eleetronics.  Related  to  this,  Kimura  et  al  (1999)  have 
successfully  used  MEMS  proeesses  to  fabricate  a  microscale  wall  shear  stress  sensor  array 
composed  of  100  thermal  sensors,  and  used  this  array  under  laboratory  conditions  to  mea¬ 
sure  instantaneous  spanwise  distributions  of  surfaee  shear  stress.  Their  sensors  were 
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spaced  300  pm  apart,  and  the  sensitive  area  -  defined  by  the  membrane  where  the  polysil- 
icon  hot  wire  was  loeated  -  measured  200  x  200  pm  .  Their  array  was  arranged  in  four 
rows  eonsisting  of  25  sensors  each,  however  these  rows  were  spaeed  over  2.5  mm  apart. 
The  relatively  large  size  of  their  sensors  and  the  large  separation  between  sensor  rows  lim¬ 
its  their  array  to  wall  shear  stress  measurements  at  relatively  low  Reynolds  numbers,  for 
whieh  the  spacing  between  sublayer  structures  beeomes  aceordingly  large  (as  noted  in  §2). 
These  authors  have  also  reported  that  the  frequency  response  of  their  sensors  was  signifi¬ 
cantly  improved  when  an  insulating  vacuum  cavity  was  used  under  the  membrane  on 
whieh  the  hot-film  element  is  located,  sinee  this  greatly  reduces  heat  losses  to  the  substrate 
material.  However  the  bandwidth  of  their  final  sensor  array  remains  unclear,  since  no 
direct  evaluation  of  its  time  response  was  reported. 


6.2  Heat  Transfer  Considerations 

Thermal  sensors  take  advantage  the  fact  that  the  heat  transfer  from  a  sufficiently  small 
heated  surface  depends  only  on  the  flow  eharacteristics  in  the  viscous  region  of  the  bound¬ 
ary  layer.  The  heat  generated  electrieally  in  the  device  is  conveetively  transferred  to  the 
fluid  and  conduetively  transferred  to  the  substrate,  so  that  in  general 

Q  =  Qf+Qs,  (6.1) 

where  represents  the  heat  transfer  rate  to  the  substrate  and  represents  the  heat  trans¬ 
fer  rate  to  the  fluid.  The  latter  ineludes  both  the  direct  heat  transfer  from  the  heated  ele¬ 
ment  to  the  fluid  and  the  indirect  heat  transfer  to  the  fluid  from  the  heated  portion  of  the 
substrate  that  is  exposed  to  the  fluid.  The  fundamental  assumption  relating  heat  transfer  to 
shear  stress  is  that  the  thermal  boundary  layer  that  develops  over  the  device  lies  entirely 
within  the  linear  region  of  the  velocity  profile.  If,  in  addition,  it  is  also  assumed  that  the 
magnitude  of  veloeity  fluctuations  increases  linearly  with  distance  from  the  wall,  then  it  is 
possible  to  relate  the  fluctuating  shear  stress  to  the  fluctuations  in  the  heat  transfer  rate  Qf 
to  the  fluid.  The  resulting  relation  between  the  wall  shear  and  the  heat  transfer  rate,  per 
unit  width  of  the  surface,  is  given  by 
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Nu  =  0.807(Pr(r)')'"'  ,  (6.2) 

where  Nu  is  the  Nusselt  number,  Pr  is  the  Prandtl  number  and  is  the  streamwise  extent 
of  the  heated  surfaee  normalized  by  wall  variables  (Lofdahl  &  Gad-el-Hak  1999).  The 
eombination  oi  PrP  is  the  Peelet  number  Pe  based  on  the  streamwise  length  of  the  film 
and  the  veloeity  at  the  edge  of  the  thermal  boundary  layer.  The  Nusselt  number  in  (6.2)  is 
defined  as 


Nu 


Qf 

wkAT’ 


(6.3) 


where  Qf  is  the  total  rate  of  heat  transfer  from  the  hot-film  element  to  the  fluid,  w  is  the 
width  of  the  film  (normal  to  the  flow  direetion),  k  is  the  eoeffieient  of  thermal  eonduetivity 
of  the  fluid,  and  AT  is  the  temperature  differenee  between  the  film  and  the  fluid. 


Sinee  part  of  the  heat  transferred  from  the  film  goes  direetly  into  the  fluid,  while  the  rest 
first  heats  the  membrane  and  the  substrate  whieh,  in  turn,  are  subsequently  eooled  by  the 
flow,  there  are  two  observations  that  are  eentral  to  the  design  of  sueh  sensors.  First,  the 
heat  transfer  to  the  membrane  on  which  the  hot-film  element  is  located  participates  in  the 
process  by  producing  an  increase  in  the  effective  value  of  P.  Second,  as  the  fraction  of  the 
total  heat  transferred  to  the  membrane  and  the  substrate  increases,  the  overall  sensitivity  of 
the  resulting  device  to  changes  in  the  wall  shear  decreases. 


In  order  to  address  these  points,  it  is  useful  to  rewrite  (6.1)  in  a  more  explicit  form  for  the 
case  where  the  film  is  heated  electrically.  When  the  hot-film  element  is  operated  in  con¬ 
stant-temperature  mode,  the  temperature  difference  AT  is  fixed  and  the  heat  transfer  rate  is 
proportional  to  the  square  of  the  voltage  across  the  film.  The  input  power  P  to  the  sensor 
can  then  be  written  as 

P  =  AT[h{x„)  +  k] ,  (6.4) 

where  the  power  dissipated  by  conduction  to  the  substrate  is  represented  by  a  generalized 
conduction  coefficient  k,  and  the  power  dissipated  by  convection  to  the  fluid  is  repre¬ 
sented  by  a  generalized  convection  coefficient  /i(t^).  The  heat  transfer  rate  to  the  fluid  is 
related  to  the  wall  shear  stress  according  to  the  relation  ~  QI  (Lofdahl  &  Gad-el- 
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Hak  1999).  Thus  in  (6.4)  can  be  written  as 

^('tw)  =  ,  (6.5) 

where  the  faetor  ^4  is  given  by  Liu  et  al  (1999)  as 

A  =  0.807^,(^^)  ,  (6.6) 

with  Ag  denoting  the  effective  area  of  the  thermal  element,  Cp  the  heat  eapacity  of  the 
fluid,  k  the  thermal  eonduetivity  of  the  fluid,  L  the  streamwise  length  of  the  resistive  heat¬ 
ing  element,  and  p  the  fluid  viseosity. 


From  (6.4)-(6.6),  for  any  fixed  input  power  P,  as  k  inereases  relative  to  A{x„y^^  the  sensi¬ 
tivity  to  the  wall  shear  stress  will  be  redueed.  Thus  effeetive  shear  stress  sensor  design 
must  minimize  the  heat  loss  to  the  substrate.  For  this  reason,  the  proposed  sensor  design 
improves  the  frequency  response  and  sensitivity  of  thermal  sensors  by  using  a  vaeuum 
eavity  underneath  the  thermal  sensor,  as  has  been  reported  by  Breuer  (2000),  Huang  et  al 
(1996),  and  Kimura  et  al  (1999).  Sinee  the  thin  film  element  is  then  separated  from  the 
substrate  by  a  vacuum  cavity  and  a  membrane,  the  heat  will  still  be  eondueted  from  the 
film  to  the  membrane  and  then  from  the  membrane  to  the  substrate.  Following  Huang  et 
al  (1999),  these  two  eonductive  heat  transfer  terms  ean  be  written  as 


and 


qs 


kJjT-T,) 

ds 


(6.7) 


_  2k4UTi-T,) 

qi - - - ,  (6.8) 

where  and  q^  represent,  respeetively,  the  eonduetive  heat  transfer  rate  from  the  sensor  to 
the  membrane  and  from  the  membrane  to  the  substrate,  while  A:,  and  k^  are  the  thermal  eon- 
duetivities  of,  respectively,  the  membrane  and  the  sensor  element,  and  is  the  half  length 
of  the  eavity.  The  faetor  of  2  in  (6.8)  is  a  result  of  the  bidireetionality  of  the  heat  transfer 
aeross  the  membrane.  An  initial  estimate  for  k  ean  be  obtained  by  evaluating  (6.8)  and 
assuming  an  ideal  membrane  with  no  g^-losses.  For  a  more  detailed  evaluation  of  k, 
detailed  heat  transfer  simulations  such  as  those  of  Huang  et  al  (1999)  would  be  required. 
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The  voltage  V  applied  across  the  film  element  can  then  be  obtained  from  (6.4)  as 


V"  =  AT  ■  R  ■  {h{x„)  +  k)  . 


(6.9) 


The  sensitivity  of  the  device  is  obtained  by  differentiating  (6.9)  with  respect  to  the  wall 
shear  stress  t:„  ,  giving 


1  dV  ^  1  dh  AT 
Vdx^  P 


(6.10) 


As  is  apparent  from  this  result,  the  sensitivity  is  increased  by  maximizing  the  overheat 
temperature  while  minimizing  the  dissipated  power.  This  can  be  done  in  practice  by  using 
small  thermal  masses  and  materials  with  a  high  temperature  coefficient  of  resistivity 
(TCR),  namely  a  high  temperature  difference  for  a  given  resistive  overheat  setting. 


6.3  Sensor  Array  Design 

The  thermal  losses  to  the  substrate  can  be  minimized  as  demonstrated  by  Huang  et  al 
(1996)  by  mounting  the  sensor  on  a  thin  membrane  over  a  vacuum-backed  cavity.  In  this 
case,  the  leakage  path  is  restricted  to  conduction  along  the  surface  of  the  membrane,  as 
indicated  by  Breuer  (2000).  This  is  demonstrated  by  assuming  a  sensor  with  a  characteris¬ 
tic  radius  i?,-  mounted  on  a  membrane  of  radius  R^.  The  temperature  distribution  can  then 
be  modeled  as  a  function  of  radial  distance  r  across  the  circular  membrane.  The  sensor 
element  is  approximated  as  maintaining  a  temperature  7)-  near  the  center  of  the  membrane, 
at  r  <  R^ ,  and  the  temperature  decreases  radially  to  the  ambient  value  T^  at  the  edge  of  the 
membrane,  at  r  =  R^.  The  resulting  one-dimensional  temperature  distribution  across  the 
membrane  is  then  simply 

This  is  shown  in  Fig.  6.3  for  a  typical  membrane  with  R^  =  50  pm. 

From  (6.11),  the  total  heat  flux  q  from  the  periphery  of  the  membrane,  at  r  =  7?^ ,  as  shown 
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by  Breuer  (2000),  is 


q  oc 


(6.12) 


where  t  is  the  membrane  thiekness.  To  minimize  these  losses  for  a  given  temperature  dif- 
ferenee  and  a  fixed  sensor  size,  the  faetor  tlh\iRJR^  in  (6.12)  needs  to  be  made  as  small 
as  possible.  On  the  other  hand,  this  is  limited  by  the  faet  that  the  membrane  also  serves  a 
struetural  function  as  one  of  the  vacuum  cavity  walls,  and  thus  needs  sufficient  strength  to 
support  the  pressure  difference.  This  pressure  difference  must  be  balanced,  as  indicated 
by  Breuer  (2000),  by  the  biaxial  stress  s  in  the  membrane,  given  by 


s^{s-So) 


(6.13) 


where  M  and  are  the  biaxial  elastic  modulus  and  residual  membrane  stress,  and  p  is  the 
external  pressure,  which  in  most  cases  would  be  one  atmosphere.  For  an  assumed  maxi¬ 
mum  upper  biaxial  stress  value  at  which  the  membrane  will  fail,  the  ratio  of  1  will  be 
constant.  A  calculation  was  carried  out  for  the  proposed  silicon  nitride  membrane,  giving 
Rq!  t  =  1169.  Most  of  the  values  for  constants  involved  in  this  calculation  were  obtained 
from  the  MEMS  Material  Database  at  the  MEMS  Clearinghouse  web  site,  where  material 
property  data  are  collected  from  publications  in  the  principal  MEMS  journals  and  IEEE 
conferences,  and  from  the  SUMMIT  Material  Properties  from  Sandia  National  Laborato¬ 
ries.  In  general  the  material  properties  from  these  two  sources  agreed  reasonably  well. 
Note  that  different  processes  used  to  deposit  the  membrane  material  yield  different  mate¬ 
rial  properties.  Eor  the  present  sensors,  the  silicon  nitride  membrane  was  designed  to  be 
deposited  using  LPCVD  (low  pressure  chemical  vapor  deposition)  which  yields  lower 
internal  stresses.  Although  there  remains  some  uncertainty  about  some  of  the  values  used 
here,  the  fundamental  sensor  design  and  the  resulting  performance  assessments  presented 
below  should  remain  valid. 


In  addition  to  effects  of  the  membrane  radius-to-thickness  ratio  Rg!  t,  a  further  consider¬ 
ation  is  the  maximum  deflection  which  the  membrane  will  experience  when  a  given  pres¬ 
sure  difference  exists  across  the  membrane  due  to  the  cavity  vacuum.  The  pressure 
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difference  across  the  membrane  due  to  the  cavity  vacuum  will  produce  a  deflection  of  the 
membrane,  as  shown  in  Fig.  6.4a.  Several  studies  have  analyzed  the  load-deflection 
behavior  of  such  a  membrane  {e.g.,  Tabata  et  al  1989).  The  approach  usually  taken  uses 
an  energy-minimization  technique  that  includes  the  contribution  of  intrinsic  tensile  stress 
and  provides  a  simple  closed- form  solution.  For  the  general  case  of  an  square  membrane, 
the  resulting  load-deflection  relationship  is  given  by  Tabata  et  al  (1989)  as 


P  = 


CiOth  ^  CiEth^ 


(6.14) 


where  P  is  the  pressure  difference  across  the  membrane,  E  is  Young’s  modulus  of  the 
membrane  material,  a  is  the  internal  film  stress,  t  is  the  membrane  thickness,  and  h  the 
resulting  membrane  deflection.  In  (6.14),  Cj  and  C2  are  energy  minimization  constants 
and  are  functions  of  the  membrane  shape.  For  the  present  membrane,  and  C2  are 
approximately  3.04  and  1.83,  respectively.  Figure  6.46  shows  the  resulting  pressure  P 
across  the  membrane  required  to  produce  a  given  deflection  h  for  membrane  thicknesses  t 
of  1  pm,  0.5  pm  and  0.2  pm.  The  results  show  that  for  a  membrane  thicknesses  below 
t «  0.5  pm ,  the  membrane  will  deflect  at  least  0.6  pm  when  a  pressure  difference  of  one 
atmosphere  is  applied  across  it.  As  a  design  constraint,  a  deflection  below  0.6  pm  is 
desired,  since  this  requires  the  vacuum  cavity  height  to  be  no  larger  than  2  pm,  and 
thereby  allows  comparatively  simple  MEMS  photolithographic  fabrication  methods  to  be 
used,  of  the  same  type  as  were  used  to  fabricate  the  top  layer  of  the  MEKA-5  hydronauti- 
cal  array  in  §5.  A  schematic  of  the  final  shear  stress  sensor  design  is  shown  in  Eig.  6.5. 


Erom  the  discussion  above  of  the  temperature  distribution  across  the  membrane,  the  mem¬ 
brane  stresses,  and  the  membrane  deflection,  as  well  as  the  maximum  sensor  spacing  in 
Eig.  6.1,  a  membrane  thickness  of  0.5  pm  and  radius  of  50  pm  has  been  chosen  for  ther¬ 
mal  sensor  array  design.  Sensors  of  this  size  could  be  readily  integrated  as  suggested  in 
Eig.  3.2  in  the  unit  cells  on  MEKA-5  hydronautical  array.  To  estimate  the  performance  of 
these  sensors,  the  resistance  of  the  heating  element  R  can  be  approximately  related  to  tem¬ 
perature  rby  a  linear  relation  R  =  i?o[l  +  a{T -  T^)] ,  as  noted  by  Huang  et  al  (1996). 
Here  R^  and  7),  are  the  average  resistance  and  temperature  at  a  reference  condition  (usually 
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ambient  values)  and  a  is  the  TCR  (temperature  coeflhcient  of  resistivity)  of  the  hot-film 
material.  An  important  parameter  governing  the  operation  of  such  hot-film  sensor  is  the 
overheat  ratio,  defined  as  =  {T-  T„)/ ,  where  T  is  the  temperature  of  the  heated  ele¬ 
ment.  In  operation,  it  is  more  practical  to  use  a  resistive  overheat  ratio,  defined  by 
=  (R-  Ro)/Ro,  where  R  is  the  resistance  of  the  heated  element.  The  relationship 
between  the  two  overheat  ratios  is  simply  =  aT^a-p  .  A  resistive  overheat  ratio  of 
«  0.2  is  typical  {e.g.,  Huang  et  al  1995;  Liu  et  al  1999).  In  general,  as  this  value  increases 
the  power  consumption  of  the  sensor  also  increases  and  the  frequency  response  is 
decreased.  On  the  other  hand,  if  is  too  low  then  this  will  decrease  the  steady-state  sen¬ 
sitivity  of  the  device.  Here  a  value  of  0.2  will  be  used  to  permit  comparisons  with  the 
results  from  these  previous  studies.  Equation  (6.9)  can  then  be  written  as  a  function  of  a«, 
a,  and  R^  to  give 

E"  =  (Ht‘"'+5)— (a«+ 1)R,  .  (6.15) 

a 

where  A  was  defined  in  (6.6)  and  B  is  given  by  in  (6.8)  divided  by  A f. 


6.4  Sensor  Performance 

The  considerations  above  suggest  graphing  the  ratio  of  input  power  to  temperature  differ¬ 
ence  against  as  shown  in  Fig.  6.6a.  The  slope  of  this  line  then  gives  the  value  of  A, 
and  the  intercept  with  they  axis  gives  the  value  of  B.  The  resulting  values  of  these  param¬ 
eters  for  the  sensor  design  developed  above  are  summarized  in  Table  6.1  and  compared  to 
those  obtained  by  Liu  et  al  (1999). 

In  general,  larger  values  of  A  correspond  to  sensors  with  higher  sensitivity.  Note  that  the 
theoretical  value  for  A  obtained  here  is  clearly  higher  than  the  values  obtained  experimen¬ 
tally  by  Liu  et  al  (1999)  when  they  operated  their  sensor  design  in  constant  current  (CC) 
and  constant  voltage  (CV)  modes,  but  is  about  a  factor  of  two  lower  than  their  result  for 
constant  temperature  (CT)  operation.  It  is  important  to  note  that  the  actual  value  of  A  for 
any  particular  sensor  depends  on  a  number  of  factors  beyond  the  design  of  the  sensor 
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itself,  such  as  the  performance  characteristics  of  amplifiers  and  filters  used  together  with 
the  sensor  in  the  circuit.  As  a  result,  the  present  calculated  theoretical  sensor  values  are 
only  approximations  for  the  actual  performance  that  would  be  achieved  experimentally.  If 
the  proposed  sensor  array  were  fabricated  and  tested,  it  is  likely  that  it  would  show  perfor¬ 
mance  parameters  A  and  B  closer  to  those  found  experimentally  by  Liu  et  al  (1999). 

As  regards  the  value  for  B,  which  represents  the  heat  lost  by  conduction  to  the  membrane 
and  to  the  substrate,  lower  values  of  B  provide  greater  sensitivity  of  the  sensor.  The 
present  theoretical  results  for  the  proposed  sensor  design  are  about  a  factor  of  three 
smaller  than  the  experimental  values  obtained  by  Liu  et  al  for  their  sensor  design,  consis¬ 
tent  with  the  fact  that  the  losses  from  conduction  to  the  membrane  were  not  included  in  the 
present  calculations.  As  noted  above,  it  is  furthermore  likely  that  additional  factors  that 
impact  the  overall  performance  of  the  sensor  in  a  complete  circuit  will  also  drive  the 
actual  value  for  B  closer  to  that  measured  by  Liu  et  al  for  their  sensor  design.  Neverthe¬ 
less,  the  results  in  Fig.  6.6a  suggest  that  the  proposed  sensor  design,  which  is  sized  to  per¬ 
mit  integration  in  the  present  MEKA-5  array,  would  have  performance  characteristics 
comparable  to  current  state-of-the-art  wall  shear  stress  sensors. 

Figure  6.66  shows  the  theoretical  predictions  for  the  measured  voltage  across  the  sensor 
element  for  various  values  of  the  wall  shear  stress  Direct  comparisons  of  this  result 
with  values  for  other  sensor  designs  reported  in  the  literature  are  not  possible  without  con¬ 
sideration  of  the  characteristics  of  the  circuit  in  which  the  sensor  is  operated,  as  noted 
above.  Figure  6.6c  shows  the  resulting  theoretical  sensitivity  of  the  proposed  sensor 
design  from  (6.10).  The  relatively  low  sensitivity  at  large  values  of  the  wall  shear  stress 
could  be  improved  by  using  a  differential  amplifier  in  the  associated  circuit.  It  is  likely 
that  the  overall  sensitivity  of  this  sensor,  when  accounting  for  the  circuit  in  which  it  is 
used,  will  be  similar  to  the  results  obtained  by  Huang  et  al  (1999),  since  the  present  values 
for  A  and  B  are  in  the  same  range  as  the  values  they  obtained. 
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6.5  Sensor  Fabrication  Process 


The  major  process  steps  by  which  the  proposed  microscale  wall  shear  stress  sensor  can  be 
fabricated  are  shown  in  Fig.  6.7,  and  have  been  largely  adapted  from  Huang  et  al  (1999), 
although  some  changes  were  made  to  the  geometry  and  type  of  electrodes  used.  One  key 
modification  is  made  in  Step  3  in  Fig.  6.7,  where  a  1.3  pm  layer  of  thick  silicon  dioxide  is 
grown  using  thermal  oxidation  at  1050  C  for  4  hours,  as  described  by  Wolf  &  Tauber 
(1986).  During  the  oxidation  process,  the  silicon/silicon  dioxide  interface  moves  further 
into  the  substrate.  At  the  end  of  this  oxidation  process,  44%  of  the  1 .3  pm  oxide  thickness 
will  be  contributed  by  oxidation  below  the  original  silicon  surface,  as  reported  by  Wolf 
and  Tauber.  Breuer  et  al  (2000)  describe  a  somewhat  different  approach  for  fabricating 
the  cavity  itself  that  provides  an  alternative  to  the  high  vacuum  nitride  deposition  step 
shown  in  Fig.  6.7.  Their  process  requires  a  “device”  wafer  that  is  patterned  with  the  cavity 
hole,  and  a  second  “capping”  wafer  that  is  coated  with  the  silicon  nitride  membrane.  The 
two  wafers  are  then  bonded  together  and  the  capping  wafer  is  etched  until  only  the  device 
wafer  with  the  nitride  film  remains.  The  metal  layers  in  Step  8  can  be  formed  using  e- 
beam  deposition  as  indicated  by  Breuer  (2000).  A  100  A  layer  of  titanium  can  be  used  as 
an  adhesion  layer  for  a  1000  A  layer  of  platinum.  Finally,  a  5000  A  layer  of  gold  would  be 
deposited  on  top  of  the  Pt/Ti.  All  three  metals  can  be  deposited  sequentially.  The  sensors 
can  then  be  defined  by  patterning  and  etching  (with  potassium  iodide)  the  top  layer  of  gold 
in  selected  areas  to  reveal  the  Pt/Ti  layer  below.  This  process  would  produce  a  sensor  of 
the  design  shown  in  Fig.  6.5. 
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Mode 

A 

B 

Current  design  (theory) 

0.25 

15 

Liu  et  al  CC  (experimental) 

0.08 

42.7 

Liu  et  al  CV  (experimental) 

0.10 

42.7 

Liu  et  al  CT  (experimental) 

0.56 

42.6 

Table  6.1.  Tabulated  values  for  A,  B  (defined  in  §6.4)  and  ealculated  from  Fig.  6.6a 
for  the  eurrent  design  and  eompared  to  those  obtained  by  Liu  et  al  (1999).  (CC  =  eon- 
stant  eurrent,  CV  =  eonstant  voltage,  CT  =  eonstant  temperature). 
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Full  Sensor  Length  (jam) 


Figure  6.1.  Microsensor  length  and  frequency  requirements  for  vehicle  types  in  table 
2.1.  These  sensor  requirements  will  allow  the  detection  of  the  incipient  bursting  of  the 
streamwise  vortical  structures  that  naturally  occur  in  a  turbulent  boundary  layer  as 
described  in  §2.  The  design  allows  for  collocation  with  the  MEKA  actuators  shown  in 
§5. 
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Figure  6.2.  Schematic  of  a  wall  shear  stress  thermal  sensor  of  length  L  embedded  on  a 
wall  and  capable  of  detecting  the  incipient  bursting  signature  of  the  coherent  structures 
in  a  turbulent  boundary  layer.  The  heat  transfer  losses  from  the  sensor  to  the  substrate 
are  shown.  The  heat  transfer  to  the  fluid  is  shown  by  Qf. 
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r  (|im) 


Figure  6.3.  Radial  temperature  distribution  are  shown  for  a  thermal  sensor  eomposed 
of  a  heated  thin  film,  at  a  temperature  =  350  K,  loeated  in  the  eenter  of  a  50  |um  cireu- 
lar  silicon  nitride  membrane. 
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Figure  6.4.  (a)  Cross  sectional  view  of  a  thermal  sensor  showing  a  silicon  nitride 

membrane  of  thickness  t,  with  a  sealed  vacuum  cavity  underneath,  being  deflected 
by  a  pressure  load,  and  (b)  maximum  membrane  deflection  at  its  center  point  for 
different  pressure  loads. 
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2  |jm  thick 
resistor 


Silicon 

substrate 


50  urn  radius 
membrane 


0.5  urn  thick 
membrane 


Figure  6.5.  Schematic  showing  a  thermal  sensor,  (a)  Top  view  shows  the  gold  elec¬ 
trodes  and  platinum  resistor  located  on  top  of  a  circular  membrane,  and  (b)  the  cross 
sectional  view  of  the  sensor  shows  also  the  silicone  membrane  with  a  vacuum  cavity 
underneath  used  to  increase  the  sensitivity  of  the  sensor  by  reducing  the  heat  transfer 
losses. 
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(Pa) 


Figure  6.6.  The  performance  of  the  designed  thermal  sensor  is  characterized  by 
(a)  the  ratio  of  input  power  to  temperature  difference  against  (b)  the  theoreti¬ 
cal  predictions  for  the  measured  voltage  across  the  sensor  element  for  various  val¬ 
ues  of  the  wall  shear  stress  x„  and  (c)  the  resulting  theoretical  sensitivity  of  the 
proposed  sensor  design  from  (6.10). 


158 


(a)  Silicon  nitride  is  deposited  by  LPCVD 


(b) 


Photolithographically  patterned  to  define  cavities  (100  pm  x  100  pm  ) 


Nitride 


Substrate 


Nitride 


Substrate 


( c)  Silicon  dioxide  is  grown  using  thermal  oxidation  (1.3  pm) 


Nitride 
Silicon  oxide 
Substrate 


(d)  500-nm  LPCVD  sacrificial  phosphosilicate  glass  (PSG) 


PSG 


Nitride 
Silicon  oxide 
Substrate 


Figure  6.7.  Major  fabrication  steps  required  for  the  production  of  the  designed  wall 
shear  stress  sensor.  These  steps  are  based  on  those  used  by  Huang  et  al  (1999)  with 
some  modifications  outlined  in  §6.5. 
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(e)  Patterned  and  unmasked  PSG  is  etched  with  hydrofluoric  acid 
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Silicon  oxide 
Substrate 


(f) 


Silicon  nitride  is  deposited  (1.2  |im) 


Nitride 
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(g)  Etch  hole  with  SF6  plasma.  Etch  PSG  and  nitride  with  HE 

etch  hole 


Nitride 


Substrate 


^  Seal  hole  and  deposite  electrodes 

Gold  leads  (1000  A),  Platinum  resistor  (1000  A) 


Nitride 


Substrate 


Figure  6.7.  (Caption  on  previous  page). 
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CHAPTER  7 


CONCLUSIONS 


Major  results  and  conclusions  from  the  present  study  of  electrokinetic  microactuator 
arrays  for  active  sublayer  control  in  turbulent  boundary  layers  can  be  summarized  as  fol¬ 
lows: 

1.  The  present  study  has  been  the  first  to  examine  the  electrokinetic  principle  as  the 
basis  for  a  new  class  of  microscale  actuator  arrays  for  active  sublayer  control  on  full 
scale  aeronautical  and  hydronautical  vehicles  under  realistic  operating  conditions; 
the  Helmholtz-Smoluchowski  scalings  that  govern  such  electrokinetic  actuator 
arrays  show  significant  performance  advantages  from  their  miniaturization  to  the 
microscale. 

2.  Specific  performance  requirements  for  microactuator  spacing,  flow  rate,  and  fre¬ 
quency  response  for  active  sublayer  control  have  been  determined  from  fundamental 
scaling  laws  for  the  streamwise  vortical  structures  in  the  sublayer  of  turbulent 
boundary  layers. 

3.  In  view  of  the  inherently  local  nature  of  the  sublayer  dynamics,  a  general  system 
architecture  for  microactuator  arrays  appropriate  for  active  sublayer  control  has  been 
developed  based  on  the  concept  of  relatively  small  and  independent  “unit  cells”, 
each  with  their  own  sensing,  processing,  and  actuation  capability,  that  greatly  simpli¬ 
fies  the  sensing  and  processing  requirements  needed  to  achieve  practical  sublayer 
control. 
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4.  A  key  innovation  in  this  system  arehiteeture  is  the  possibility  of  replaeing  previously 
elaborate  eontrol  proeessing  requirements  with  a  simple  look-up  table  approaeh 
implemented  at  the  unit-cell  level,  in  which  a  sensor  state  vector  provides  the  input 
to  a  simple  programmable  logic  array  (PLA)  that  determines  the  corresponding  actu¬ 
ator  state  vector;  such  an  approach  is  made  possible  by  the  comparatively  small 
number  of  sensors  and  actuators  within  each  unit  cell. 

5.  A  microscale  wall  shear  stress  sensor  has  been  proposed  that  meets  the  performance 
requirements  for  sublayer  control  on  full-scale  vehicles;  individual  sensor  states  are 
determined  from  the  sensor  output  by  an  averaging  and  thresholding  procedure  that 
essentially  eliminates  the  usual  concerns  regarding  drift  in  the  sensor  calibrations, 
and  effects  caused  by  changes  in  the  vehicle  speed  and  attitude. 

6.  A  fundamental  three-layer  design  has  been  developed  for  such  electrokinetic  micro¬ 
actuator  arrays,  in  which  electrokinetic  flow  is  induced  by  an  impulsively  applied 
electric  field  across  a  center  layer,  with  a  bottom  layer  containing  an  electrolyte  res¬ 
ervoir  and  a  common  electrode,  and  a  top  layer  that  contains  individual  electrodes 
and  lead-outs  for  each  microactuator  in  the  unit  cell. 

7.  A  recently-developed  porous  polymer  fabrication  technology  has  been  adapted  to 
permit  formation  of  fine-pore  matrix  structures  suitable  for  electrokinetic  pumping 
in  large,  dense  arrays  of  microchannels  located  in  the  center  layer;  porous  polymer 
matrix  structures  with  1  qm  average  pore  sizes  can  be  readily  created  within  the 
electrokinetic  pumping  channels. 

8.  Measured  volume  flow  rates  produced  by  such  porous  polymer-filled  electrokinetic 
microchannels  under  the  effect  of  steady  applied  electric  fields  have  been  shown  to 
meet  the  requirements  for  active  sublayer  control  on  full-scale  vehicles;  the  thin  lay¬ 
ers  into  which  these  electrokinetic  microactuators  can  be  fabricated  permits  the 
required  field  strengths  to  be  achieved  with  potential  differences  of  the  order  of  15- 
20  V. 
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9.  Simple  hydrodynamic  models  have  been  developed  for  the  frequency  response  of 
such  electrokinetic  microactuators  under  the  effect  of  sinusoidally  or  impulsively 
applied  fields;  these  show  that  electrokinetically-induced  flow  in  1  pm-scale  porous 
polymer  matrix  structures  achieve  essentially  quasi-steady  performance  at  frequen¬ 
cies  up  to  1  MHz. 

10.  Measurements  of  the  actual  frequency  response  achieved  in  such  electrokinetic  actu¬ 
ators  have  verified  essentially  loss-less  AC  performance  at  frequencies  as  high  as  10 
kHz;  the  observed  roll-off  at  higher  frequencies  may  be  due  to  inhomogeneities  in 
the  porous  polymer  structure  that  lead  to  inertial  damping  at  these  frequencies. 

11.  Microfabrication  techniques  have  been  developed  that  permit  mass  production  of 
large  numbers  of  individual  electrokinetic  microactuators  in  unit  cells  on  compara¬ 
tively  large-area  tiles;  CO2  laser  drilling  of  individual  actuator  channels  can  be 
accomplished  in  thin  flexible  plastic  substrates  that  are  compatible  with  the  solvents 
used  in  the  porous  polymer  constituents. 

12.  Several  generations  of  such  electrokinetic  microactuator  arrays  have  been  built  lead¬ 
ing  to  the  MEKA-5  full-scale  hydronautical  array,  composed  of  25,600  individual 
electrokinetic  microactuators  with  350  pm  center-to-center  spacings,  arranged  in  a 
40  X  40  pattern  of  unit  cells,  each  composed  of  a  4  x  4  matrix  of  actuators;  this  array 
was  successfully  fabricated  in  a  7  x  7  cm^  tile  in  250  pm  thick  mylar  substrate  mate¬ 
rial. 

13.  MEMS  design  and  fabrication  processes  were  used  to  produce  a  top  layer  for  the 
MEKA-5  hydronautical-scale  array,  composed  of  an  electrode  sublayer  patterned  in 
Cr/Au,  a  nozzles  sublayer  etched  in  polyimide,  and  a  cavities  sublayer  fabricated  in 
Su-8;  the  functionality  incorporated  in  this  top  layer  fabrication  demonstrates  the 
complete  elements  needed  for  microscale  electrokinetic  actuator  arrays  for  sublayer 
control. 
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14.  Stereo-PIV  measurements  successfully  demonstrated  lateral  displacement  of  syn¬ 
thetically-generated  streamwise  vortical  structures  by  volumetric  pumping  from  a 
wall  actuator  in  a  set  of  large-scale  wind  tunnel  tests;  results  characterized  the  rela¬ 
tive  lateral  displacement  achieved  as  a  function  of  the  dimensionless  pumping 
strength  qlT.  Despite  differences  between  the  large-scale  actuator  used  in  these  mea¬ 
surements  and  the  microscale  electrokinetic  actuators  that  are  the  subject  of  this 
study,  the  results  obtained  suggest  that  manipulation  of  naturally-occurring  sublayer 
vortices  in  turbulent  boundary  layers  by  volumetric  pumping  with  wall  actuators  is 
possible  at  comparable  dimensionless  pumping  strengths  qlT. 
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APPENDIX  A 


SLOT-TYPE  MICROACTUATOR  ARRAYS 


For  open-loop  control  based  on  the  oscillating  wall  approach,  volume  displacement  is 
through  slots  nominally  aligned  with  the  near-wall  flow.  Such  slot  microactuators,  shown 
in  Figure  A.l  are  between  adjacent  slot  pairs.  Adjacent  slots  pairs  produce  alternatively 
positive  or  negative  volume  displacements,  in  a  time-oscillating  fashion,  to  displace  fluid 
above  the  wall  along  the  cross-stream  direction.  The  volume  displacements  are  again  cre¬ 
ated  by  electrokinetic  pumping  within  the  much  larger  microactuator  channels  in  a  center 
layer.  In  this  case  the  slot  spacing  is  less  critical,  and  key  performance  requirements  are 
the  frequency  and  flowrate  needed  to  achieve  effective  spanwise  oscillation  of  the  fluid 
above  the  wall. 

For  the  open-loop  control  approach  based  on  slot  microactuators  it  is  apparent  that  the 
spacing  of  such  actuators  is  dependent  on  the  number,  N,  of  sublayer  vortices  typically 
located  between  adjacent  slot  pairs.  For  the  present  case  lets  choose  N  =  1  so  we  can  com¬ 
pare  the  flow  rate  requirements  for  the  slot  actuators  and  single  point  actuators. 

The  following  assumptions  are  made  while  calculating  the  flow  rate  for  a  slot  actuator: 

1. -  We  expect  to  move  the  midpoint  between  slots  a  distance  of  10  y'*’  back  and  forth 
between  cycles. 

2. -  The  Reynolds  number  should  be  chosen  small  to  prevent  inertia  effects. 

3. -  A  square  wave  voltage  will  be  used  from  now  on  as  the  driver  of  the  electrokinetic 
actuator. 

4. -  The  slots  will  be  treated  as  a  potential  source/sink  for  the  purpose  of  calculating  the 
flow  rate  and  exit  velocity. 
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The  flow  velocity  at  a  distance  r  from  a  potential  source  is  given  by 

9  Q  1 

Ur  =  2;^- 

2n  r 


(A.1.1) 


where  2Q,  the  source  strength,  is  the  flow  rate  per  unit  length  and  r  is  the  radial  distance 
from  the  source.  Notice  that  we  use  2Q  to  denote  the  total  flow  rate  per  unit  length  issuing 
from  the  source.  Figure  A.2  shows  an  ideal  view  of  the  slots  actuating. 

Thus,  the  mid  point  velocity  is  given  by 


•^midpoint 


^  Q  J _ -Q  1 

n A/2  n  A/2 


(A.  1.2) 


where  A  is  the  distance  between  slots. 
Simplifying  (A.  1.2)  we  obtain. 


•^midpoint 


_  4Q 


(A.  1.3) 


This  Uj^jjjpojnt  is  constant  during  the  oscillation  period  T/2. 

During  time  T/2  we  need  to  move  the  midpoint  a  distance  10  wall  units. 


lOL  =  u 


substituting  for  Uj^idpoint 


and  solving  for  Q  we  obtain. 


midpoint  r\  9 


Q  =  5kU^. 


(A.  1.4) 


(A.  1.5) 


(A.  1.6) 


We  have  chosen  T  =  100"'’,  although  further  checks  should  be  done  to  see  if  this  is  the 
appropriate  value. 


T  =  100 


'tw/p 


100 -, 
V 


(A.  1.7) 


where  fy  is  given  by 


/('Cw/p) 


(A.1.8) 
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Substituting  (A.  1.7)  into  (A.  1.6)  we  obtain 

Q  =  SttvN,  (A.  1.9) 

where  N  is  the  number  of  vortiees  between  slots  and  is  defined  as 


N  = 


A 

1001/ 


(A.1.10) 


Using  w,  the  width  of  the  slot,  as  the  eharacteristie  length  we  obtain  the  Reynolds  number 
as  follows 


Re  =  =  2 


(A.1.11) 


V  V 

Substituting  (A.1.10)  into  (A.1.11)  we  obtain  a  Reynolds  number  that  is  only  dependent 
on  the  number  of  vortiees  alloeated  between  slots. 

Re  =  SttN.  (A.1.12) 

Let  us  assume  that  the  gap  between  the  eleetrokinetie  beds  is  negligible  for  both  point  and 
slot  actuator.  Then,  a  sketch  of  points  and  slot  actuators  is  shown  in  Figure  A.  3.  The 
dimensions  have  been  chosen  to  simplify  the  results. 

Next,  we  will  show  that  the  flow  rate  that  two  point  actuators  need  to  deliver  in  order  to 
successfully  cancel  or  modify  the  sublayers  streaks  is  a  factor  of  exactly  75  times  smaller 
than  the  flow  rate  a  slot  actuator  will  need  to  oscillate  the  wall  at  the  correct  frequency. 


The  volume  displaced  by  a  point  actuator  is  given  by 

Vpoint  =  |7l;Rpoint^  (A.  1.1 3) 

where  Rp^wf  ^  10  /  is  the  radius  of  this  virtual  hemisphere.  Also,  the  flow  rate  per  stroke 
is  given  by 

Qpoint  ~  Vpoint^,  (A.1.14) 

where  f  is  the  frequency  at  which  this  virtual  hemispheres  will  need  to  be  created.  On  the 
other  hand,  the  flow  rate  per  stroke  displaced  by  a  slot  actuator  is  given  by  a  virtual  half 
cylinder  as 

Qslot  =  ^^‘4 
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where  ^5/0?  =  3-162 

^point- 


Figure  A.  3  shows  a  sketeh  where  there  are  two  point  aetuators  for  every  slot  aetuator. 
Thus,  in  order  to  eompare  the  flow  rates  between  slots  and  point  actuators  we  will  need  to 
look  at  the  ratio  of 


'slot 


3L 


slot 


2Qpoint  8  R 


3 


point 


(A.1.16) 


where  (A.1.13)  and  (A.1.14)  have  already  been  substituted  in  (A.1.16).  By  substituting  L 
=  2zl  =  200\^divA  '^poinr  10  (A.1.16)  the  following  ratio  is  obtained 


'slot 


2Q 


75 


point 


(A.1.17) 


This  large  difference  in  flow  rates  between  the  point  actuator  and  the  slot  actuator  is  plot¬ 
ted  in  Figure  A.4. 


We  still  need  to  prove  for  the  slot  actuator  that  we  need  to  displace  per  stroke  an  amount  of 
flow  equivalent  to  the  volume  of  half  a  cylinder  of  radius  31.62  1,.  and  length  fix  by  the 
length  of  the  slot.  Recall  that  the  flow  rate  per  unit  length  is  given  by  (A.  1 .9).  The  assump¬ 
tion  is  made  that  the  volume  displaced  per  stroke  per  unit  length  is  given  by 


TtR  slot  1 

2  AT’ 


(A.1.18) 


where  R5]q(  is  the  radius  of  the  virtual  volume  of  a  half  a  cylinder  displaced  by  the  slot 
actuator  per  stroke.  The  oscillation  time  T  was  given  by  (A.  1.7)  as 

T  =  100  —V  =  100 -.  (A.1.19) 

'Cw/p  V 

Substituting  (A.  1.7)  into  (A.1.17),  equating  the  result  with  (A.  1.9)  and  substituting  for 
Rslot  we  get 


Rsiot  =  31.621,. 


(A.  1.20) 
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Figure  A.l  Basic  Three-layer  design  of  slot-type  microactuator  arrays  for  “oscillating 
wall”  control  approach. 
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Figure  A.2  Slot  actuator.  The  left  actuator  is  issuing  fluid  while  the  right  actuator  is 
pumping  fluid  in. 
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4  point  actuators 


2  slot  actuators 


Figure  A.3  Drawing  of  a  slot  actuator  whose  electrokinetic  bed  is  twice  the  size  of  a 
point  actuator,  N  =  1 . 
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Equivalent  Flowrate  (i^L/min) 
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Figure  A.4  Comparison  of  flow  rate  requirements  for  slot  aetuators  and  point  actuators. 
(S  =  slots  actuators  and  P  =  point  actuators) 
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